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The  dissertation  focuses  on  the  physics-based  compact  noise  modeling  of 
nanoscale  complementary  metal-oxide-semiconductor  (CMOS)  devices.  The  models  are 
primarily  developed  to  be  implemented  in  circuit  simulators,  and  are  further  applied  to 
gain  insight  into  the  noise  behavior  of  advanced  metal-oxide-semiconductor  field-effect 
transistors  (MOSFETs).  For  a comprehensive  evaluation  of  noise  performance,  the 
models  describe  all  the  important  physical  noise  sources  and  effects  which  are  observed 
in  nanoscale  MOSFETs  with  ultrathin  gate  oxides.  The  models  account  for  the  gate 
leakage  current  and  noise,  the  drain  current  and  noise,  the  cross  correlation  effects 
between  the  gate  and  the  drain,  and  quantum  mechanical  (QM)  effects  in  the  inversion 
layer. 

The  gate  leakage  current  model  for  ultrathin  nitrided  oxide  MOSFETs  is 
presented,  and  is  based  on  inelastic  trap-assisted  tunneling  (ITAT)  combined  with  a semi- 
empirical  direct  tunneling  current  model.  In  order  to  explain  the  excess  noise  behavior 


Xll 


due  to  the  gate  leakage  current,  the  low  and  high  frequency  noise  models  are  formulated 
by  using  a barrier  height  fluctuation  model  and  Lorentzian-modulated  shot  noise  model, 
respectively.  Trap-related  processes  are  the  most  likely  cause  of  excess  current  noise 
because  slow  traps  in  the  oxide  result  in  low  frequency  dissipation  in  the  conductance  of 
oxides,  and  fast  traps  produce  the  Lorentzian-modulated  shot  noise  associated  with 
generation-recombination  (G-R)  process  at  higher  frequencies.  In  addition,  the  1 If  noise 
of  the  gate  leakage  current  can  be  used  to  extract  the  trap  properties  and  its  source  is 
identified  by  the  low  frequency  plateau  in  measured  conductance  loss.  By  using  the 
partition  noise  theory  and  the  BSIM4  partition  current  model,  the  cross  correlation  effects 
between  the  gate  and  the  drain  are  modeled  at  low  frequencies,  enabling  the  overall  noise 
performance  to  be  evaluated.  Furthermore,  by  incorporating  the  effects  of  a realistic  trap 
distribution,  the  parasitic  resistance,  the  doping  profile,  and  the  quadratic  mobility 
degradation  factor  on  the  noise  behavior,  a 1 If  drain  current  noise  model  for  ultrathin 
oxide  MOSFETs  is  formulated  on  the  basis  of  the  tunneling  assisted-thermally  activated 
trapping  process  of  carriers.  For  more  accurate  and  physical  prediction  of  noise  behavior 
in  ultrathin  («1.5nm)  MOS  devices,  the  QM  effects  on  the  noise  characteristics  are  taken 
into  account.  Especially,  the  compact  noise  model  of  gate  leakage  current  is  newly 
implemented  into  BSIM4,  and  the  modeling  is  validated  through  a comparison  of 
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simulation  results  and  measurements. 


CHAPTER  1 
INTRODUCTION 

1.1  Motivation 

With  the  steady  advances  in  semiconductor  technologies,  the  capability  of  the 
integrated  circuit  (IC)  has  accelerated  in  accordance  with  Moore’s  law  [FraOl,  Iwa99, 
Tau97].  According  to  the  2001  International  Technology  Roadmap  for  Semiconductors 
(ITRS-2001),  it  is  projected  for  complementary  metal-oxide-semiconductor  (CMOS) 
technology  that  64-Gb  dynamic  random  access  memory  (DRAM)  with  gate  lengths  down 
to  10  nm  is  produced  around  2016  [SiaOl].  The  continued  scaling  of  CMOS  technology 
has  driven  progress  in  very  large  scale  integration  (VLSI)  applications.  However,  as  the 
downsizing  of  the  components  of  VLSI  circuits,  such  as  the  MOS  field-effect  transistor 
(MOSFET),  continues,  it  is  expected  to  encounter  physical  and  technological  problems 
pertaining  to  scaling  of  CMOS  devices  into  the  nanometer  regime.  One  of  the  key  issues 
is  leakage  current  and  its  associated  noise  [SiaOl]. 

For  very  small  MOSFETs,  the  tunneling  currents  present  the  limits  to  useful 
device  operation  [FraOl].  The  most  constraining  limit  to  scaling  is  the  tunneling  current 
through  the  gate  insulator.  The  resulting  gate  leakage  currents  make  the  use  of  ultrathin 
oxides  (<2nm)  impractical  in  many  device  applications  where  high  performance  and  low 
standby  power  consumption  are  required.  Moreover,  the  thin  gate  oxide  is  associated  with 
the  loss  of  inversion  charge,  and  thus  transconductance,  due  to  inversion  layer 
quantization  and  polysilicon-gate  depletion  effects  [Tau97].  Another  possible  source  of 
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Table  1.1  Low  operating  power  (LOP)  logic  technology  requirements  ( ITRS  - 2001 
Edition ) 


Year  of 
Production 

Physical 
Gate  Length 
(nm) 

Physical 

EOT*, 

1 ox, physical 

(nm) 

Electrical 

EOT*, 

lox.electrical 

(nm) 

Subthreshold 

Leakage 

Current, 

W«*(pA///m) 

Power 

Dissipation 

per 

Device** 

(W/Device) 

2001 

90 

2.0 -2.4 

3.0 

100 

3.2xl0'n 

2002 

75 

1.8 -2.2 

2.8 

100 

2.9x10'" 

2003 

65 

1.6 -2.0 

2.6 

100 

2.1xl0'n  | 

2004 

53 

1.4-  1.8 

2.4 

300 

5.2x10'" 

2005 

45 

1.2-  1.6 

2.2 

300 

4.1x10'" 

2006 

37 

1.1-  1.5 

2.1 

300 

3.3x10'" 

2007 

32 

1.0-  1.4 

1.7 

700 

6.0x10'" 

2010 

22 

0.8-  1.2 

1.5 

1000 

5.3x10" 

2013 

16 

0.7-  1.1 

1.4 

3000 

l.OxlO'10 

2016 

11 

0.6-  1.0 

1.3 

10000 

2.0xl0'10 

* t , ,=t  +(e  /e  )re" , where xnQM  is  the  actual  charge  position  from  the  interface. 

ox,elecincal  ox, physical  v ox  si  ' n ° 1 

**  Static  power  dissipation  for  an  nMOS  device  with  dimension  W/Lgale-  3. 

leakage  current  is  the  direct  subthreshold  source-to-drain  tunneling  through  the  potential 
barrier  below  the  gate  [FraOl],  as  shown  in  Table  1.1.  For  low  power  MOS  transistors, 
the  subthreshold  source-to-drain  leakage  current  ISd,ieak  increases  quite  rapidly  with 
succeeding  years,  owing  to  the  reduction  of  the  threshold  voltage.  This  increase  in  Isd,ieak 
causes  the  static  power  dissipation  per  device  to  increase  somewhat  with  device  scaling, 
despite  the  reduction  in  the  power  supply  voltage  and  the  device  dimension  with  scaling. 
Since  the  gate  leakage  current  must  remain  less  than  hd.ieak,  a reasonable  restriction  is 
imposed  on  it.  Because  of  the  large  increases  in  Isd,ieak  with  device  scaling,  it  is  expected 
that  heavily  nitrided  oxy-nitrides  are  sufficient  to  meet  the  gate  leakage  requirements 


3 


until  the  end  of  the  roadmap.  Given  the  required  equivalent  oxide  thickness  (EOT)  in 
2005,  however,  the  gate  leakage  current  targets  cannot  be  met  using  oxy-nitride  because 
of  direct  tunneling  through  the  gate  oxide.  Furthermore,  there  are  serious  issues  of 
thickness  control  and  uniformity  for  ultrathin  oxy-nitride  gate  dielectrics,  and  reliability 
issues,  especially  with  large  gate  leakage  current.  From  the  standpoint  of  static  power 
dissipation  control,  it  is  also  desirable  to  reduce  the  gate  leakage  current. 

Although  remarkable  improvements  in  CMOS  technologies  have  been  made  in 
the  past  20  years,  noise  is  always  one  of  the  most  constraining  limits  to  scaling.  Due  to  its 
origin  in  the  basic  mechanisms  which  are  responsible  for  the  conduction  of  current 
through  devices,  noise  is  inevitable  even  in  scaled  MOSFETs.  Generally,  the  continuous 
minimization  of  MOS  deices  causes  a significant  increase  in  noise  amplitude,  leading  the 
feasibility  of  nanoscale  MOS  devices  to  be  questioned  [Bou98,  Cel99,  God97].  If  the 
noise  sources  are  uniformly  distributed  over  the  gate  area,  the  noise  spectral  density  is 
inversely  proportional  to  the  area  [Sim99].  Also,  the  device-to-device  fluctuations  in  the 
1 If  noise  performance  increase  in  the  small  geometry  devices.  For  ultrathin  oxide  MOS 
devices,  it  is  obvious  that  the  gate  leakage  tunneling  current  is  the  origin  of  low 
frequency  noise  sources,  and  adversely  affect  the  noise  characteristics  [Lee02,  Mom98]. 
In  order  to  reduce  the  gate  leakage  current,  the  present  gate  oxide  dielectric  will 
eventually  be  replaced  by  alternative  materials  with  different  trap  properties  and 
potential  barriers.  In  that  case,  the  low  frequency  noise  becomes  a more  significant  issue, 
since  the  noise  is  strongly  affected  by  the  physical  properties  of  the  gate  insulator,  such  as 
trap  density,  trap  distribution,  and  tunneling  effects.  In  addition,  the  low  frequency  noise 
in  ultrathin  MOS  devices  can  be  also  used  as  a sensitive  probe  for  investigating  localized 


4 


Table  1.2  Mixed-signal  device  technology  requirements  ( ITRS  - 2001  Edition  ) 


Year  of 
Production 

Physical 
Gate  Length 
(nm) 

NMOS  RF  Device 

NMOS  Analog  Device 

T0*  (nm) 

1//  noise* 
(/iV2^m2/Hz) 

T0*  (nm) 

1 //  noise 
(/tV2/tm2/Hz) 

2001 

65 

1.3-  1.6 

500 

2.5 -7.0 

1000  | 

2002 

53 

1.2-  1.5 

500 

to 

1 

— j 

o 

500 

2003 

45 

1.1  - 1.6 

300 

2.5 -5.0 

500 

2004 

37 

0.9-  1.4 

300 

2.5 -5.0 

500 

2005 

32 

0.8-  1.3 

300 

2.5 -5.0 

300 

2006 

28 

0.7-  1.2 

200 

2.5 -5.0 

300 

2007 

25 

0.6-  1.1 

200 

2.5 -5.0 

300  | 

2010 

18 

0.5 -0.8 

150 

1.3 -3.0 

150 

2013 

13 

0.4 -0.6 

100 

1.3 -3.0 

100 

2016 

9 

0.4 -0.5 

75 

1.3 -2.5 

100 

* 1// noise  spectral  density,  at  a frequency  of  1Hz,  normalized  to  an  active  area  of  1 yrm2.  Operation 

point  taken  at  100  mV  above  the  threshold  voltage 


phenomena,  such  as  degradation  and  failure  in  IC,  because  the  noise  and/or  fluctuations 
are  influenced  by  the  presence  of  localized  defects  and  irregularities  in  the  microstructure 
of  devices  [CioOO].  Furthermore,  in  the  field  of  mixed-signal  integrated  circuits 
containing  analog  and/or  radio  frequency  (RF)  circuitry  as  well  as  the  digital  circuitry, 
continuous  focus  on  1// noise  is  imperative  to  satisfy  the  increasing  demands  on  noise 
limitation.  The  1// noise  of  digital  MOS  devices  at  a frequency  of  1Hz  is  critical  for  the 
minimum  input  signal  in  mixed-signal  circuits  (lower  boundary  for  dynamic  range)  in 
low  frequency  circuits,  mixers,  and  voltage  controlled  oscillators  (VCOs)  [SiaOl].  Table 
1.2  shows  the  mixed-signal  device  technology  requirements  (ITRS-2001  Edition). 
Although  1 If  noise  in  MOS  devices  due  to  carrier  trapping-detrapping  in  traps  at  the  Si- 
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SiC>2  interface  is  more  prominent  than  in  bipolar  or  junction  field-effect  transistors 
(JFETs),  the  noise  requirements  continue  to  be  conflicting  with  succeeding  years. 

Noise  modeling  is  important  for  understanding  the  noise  behavior  and  optimizing 
the  low  noise  designs.  The  needs  for  noise  modeling  are  driven  by  the  continuous 
reduction  of  feature  size  in  scaling  MOSFETs,  which  not  only  causes  higher  demands  on 
model  accuracy  but  also  requires  models  for  currently  not  considered  material  properties. 
Two  models  widely  used  in  circuit  simulators  for  MOSFETs  noise  are  SPICE 
(Simulation  Program  for  Integrated  Circuits  Emphasis)  and  BSIM  (Berkeley  Short- 
channel  IGFET  Model)  noise  models.  The  SPICE  model  is  largely  empirical  and  is 
considered  to  be  oversimplified.  The  lack  of  a systematic  approach  for  determining  the 
empirical  parameters  imposes  another  serious  limitation  in  using  the  model.  The  BSIM 
1 If  noise  model  for  drain  current  is  a physics-based  model  applicable  to  all  MOSFET 
operation  regimes,  which  incorporates  both  carrier  number  fluctuation  and  surface 
mobility  fluctuation  [Hun90a,  Hun90b].  However,  this  model  is  at  best  10  years  old.  In 
addition,  the  direct  tunneling  gate  current  model  of  BSIM  not  only  underestimates  the 
current  magnitude  in  ultrathin  nitrided  oxide  devices,  but  also  there  is  no  noise  model  for 
the  gate  leakage  current.  Due  to  the  impact  of  gate  leakage  current  on  the  overall  noise 
performance,  a comprehensive  noise  model  is  required  covering  each  important  noise 
source  over  the  entire  frequency  range.  Furthermore,  since  the  small  size  and  high  doping 
concentration  create  quantized  energy  levels  within  the  channel,  noise  models  accounting 
for  quantum-mechanical  (QM)  effects  are  indispensable  in  nanoscale  MOSFETs.  It 
should  also  be  noted  that  a key  difficult  challenge  across  all  the  noise  modeling  area  is 
that  of  experimental  validation.  As  devices  shrink  and  become  vulnerable,  new  noise 
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measurement  techniques  that  provide  the  necessary  information  for  the  validation  of 
models  are  critical.  It  is  clear  from  the  above  that  in  order  to  consider  advanced  CMOS 
device  and  technology  changes  made  during  the  last  decade,  the  present  noise  models 
should  be  replaced  by  more  accurate  physics-based  noise  models  which  account  for  1) 
gate  leakage  current  effects,  2)  correlation  effects  between  the  gate  and  the  drain,  3)  short 
channel  effects,  4)  realistic  trap  properties,  5)  QM  effects  observable  in  carrier-trap 
interaction,  and  6)  experimental  validation. 

1.2  Overview  of  Noise  Sources 

The  design  of  low  noise  MOSFET  devices  and/or  circuits  can  be  achieved,  if  the 
fundamental  noise  models  of  devices  are  developed  and  understood  in  the  light  of 
physical  origin  of  noise  sources.  Four  important  kinds  of  noise  sources  in  MOSFETs  are 
thermal  noise,  shot  noise,  generation-recombination  (GR)  noise,  and  1//  noise,  as 
illustrated  in  Figure  1.1. 

1)  Thermal  noise,  Nyquist  noise,  or  Johnson  noise  of  MOSFETs  is  due  to  the 
random  thermal  motion  of  the  charge  carriers  in  the  channel.  It  produces  a fluctuating 
current  at  the  device  terminals.  Under  thermal  equilibrium,  the  thermal  noise  can  be 
obtained  either  from  drift  or  diffusion  carrier  transport  mechanisms.  The  MOSFET 
thermal  noise  expression  for  the  drain  current  Id  that  is  widely  used  in  the  literature  and  in 
circuit  simulators  is 


sl(f)  = rT-wTgd0 


(i.i) 


where  k is  Boltzmann’s  constant,  T is  the  absolute  temperature,  gdo  is  the  channel 
conductance  with  zero  drain-to-source  voltage,  with  2/3  < yT  < 1 in  the  linear  region,  and 
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yT  = 2/3  in  the  saturation  region  for  long  channel  device.  Thermal  noise  is  called  a white 
noise  at  low  frequencies,  because  its  power  spectral  density  is  flat  up  to  extremely  high 
frequencies  of  over  1012  Hz.  In  addition  to  channel  thermal  noise  in  MOS  devices,  two 
additional  thermal  noise  sources  are  the  gate  poly  resistance  and  the  bulk  resistance. 

2))  Shot  noise  is  generated  when  charge  carriers  encounter  a potential  barrier  and 
cross  barriers  independently  and  at  random.  The  power  spectral  density  for  the  gate 
leakage  current  Ig  in  ultrathin  oxide  MOSFETs  is  given  by 

Ss,g(f)  = FFano  -2qlg  (1.2) 


Gate  Leakage  Current  Noise 
due  to  trap-related  process 


Figure  1 . 1 Charge  transport  and  physical  origin  of  noise  sources  in  MOSFET 
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The  power  spectral  density  is  constant  up  to/«  \lTtrans  with  Ttrans  the  transit  time  across 
the  barrier.  The  analysis  of  shot  noise  is  performed  by  introducing  the  Fano  factor  FFano 
defined  as  FFar,°=Sig/(2qIg)  where  Ig  is  the  current  leakage  flowing  in  the  gate.  Shot  noise 
suppression  (0<F<1)  is  associated  with  a negative  correlation  between  current  pulses  due 
to  Coulomb  interaction  and  the  Pauli  exclusion  principle  [Ian97,  Ian98,  IanOO].  In 
contrast,  shot  noise  enhancement  is  associated  with  a positive  correlation  between  current 
pulses. 

3)  Generation-recombination  (GR)  noise  is  due  to  fluctuations  in  the  number  N of 
electrons  and  holes  in  the  conduction  and  valence  bands,  respectively.  These  fluctuations 
in  the  number  of  free  carriers  are  caused  by  the  trapping-detrapping  of  carriers  to  either 
bulk  or  surface  defect  centers.  Each  trap  center  is  characterized  by  its  energy  level,  trap 
density,  and  carrier  capture  coefficients.  The  power  spectral  density  of  the  fluctuations  in 
N is  given  by 


Sf(/)  = 4AAF 


1 + 4 n2f2z2 


(1.3) 


where  r is  the  relaxation  time  constant  of  the  trapping  process,  AN2  is  the  variance  of 
the  fluctuations  AN,  and  / is  the  frequency.  The  spectrum  appears  to  be  constant  at  low 
frequency,  then  drops  off  as  1 If  above  the  comer  frequency,  which  is  determined  by  the 
trap  characteristic  time.  If  the  GR  noise  stems  from  one  single  trap,  it  is  called  random 
telegraph  signal  (RTS)  noise.  As  shown  in  Figure  1.1,  GR  noise  in  MOSFETs  is  due  to 
random  trapping-detrapping  of  carriers  at  the  defect  or  impurity  centers  in  the  depletion 
region  between  the  conductive  channel  and  the  semiconductor  substrate.  The  GR  noise 
component  is  usually  overshadowed  by  the  higher  1// surface  noise  component.  However, 
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there  have  been  several  reported  instances  where  special  device  processing  has  led  to  the 
GR  noise  component  dominating  [Fan86], 

4)  1 If  noise  or  flicker  noise  owes  its  name  to  its  spectral  density  that  is  almost 
inversely  proportional  to  the  frequency  over  a wide  frequency  range  down  to  very  low 
frequencies  as  follows  : 


jAF 

S)  f (/)  = KF  • ——  (1.4) 

where  y is  the  frequency  exponent,  whose  value  usually  lies  in  the  interval  between  0.7 
and  1.3,  and  KF  and  AF  are  the  magnitude  coefficient  and  the  current  exponent, 
respectively.  As  depicted  in  Figure  1.1,  the  physical  origins  of  1// noise  of  drain  current 
in  MOSFETs  can  be  explained  by  two  major  models.  The  number  fluctuation  model, 
originally  proposed  by  McWhorter,  attributes  1 If  noise  to  the  random  trapping  and 
detrapping  processes  of  charge  carriers  in  the  oxide  traps  near  the  Si-SiCh  interface 
[McW57].  It  is  assumed  that  only  traps  within  a few  kT  of  the  Fermi-level  are  effective  in 
generating  noise  and  the  channel  can  exchange  charges  with  the  oxide  traps  through 
tunneling,  as  shown  in  Figure  1.2.  The  trapping  and  detrapping  in  these  traps  produce  GR 
spectra  with  distributed  relaxation  times.  Since  the  l//or  1// noise  spectrum  is  the  result 
of  a weighted  superposition  of  RTS  events,  the  studies  of  RTS’s  have  greatly  contributed 
to  our  understanding  of  the  nature  of  1 If  noise  in  MOSFETs.  The  mobility  fluctuation 
model,  based  on  Hooge’s  relationship,  considers  1// noise  as  a result  of  the  fluctuation  in 
bulk  mobility  due  to  carrier  scattering  at  lattice  vibrations.  Since  in  this  process  no  traps 
are  needed,  this  model  implies  that  there  is  1//  noise  even  if  no  traps  are  present.  In 
addition,  a unified  model  has  been  developed  in  order  to  explain  a broader  set  of  data 
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[Ghi91,  Hun90a,  Hun90b,  Jay89].  This  model  is  based  on  the  fact  that  fluctuation  in  the 
occupancy  of  the  oxide  traps  induces  correlated  fluctuations  in  the  carrier  number  and 
surface  mobility  [Hun90c], 

From  the  standpoint  of  noise  sources  in  ultrathin  oxide  MOSFETs,  the  equivalent 
input  noise  current  power  S/,eq(f ) and  voltage  power  Sy,eq  (f)  are  given  by,  respectively 

= Sft  (/)  + s!/  (/)  + *2¥C" ' ls[  (/)  + S)1/  (/)]  (1.5a) 

& m 

v,=4-k(/)+s;y(/)]  d.5b) 

g 

o m 


Figure  1 .2  Noise  spectrum  and  trap  distributions  in  MOS  structure 
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where  Cgs  is  the  gate-to-source  capacitance  and  gm  is  the  transconductance.  The 
equivalent  input  noise  current  at  low  frequencies  is  dominated  by  the  1 If  noise 
S/'//(/)and  shot  noise  S^(f) due  to  the  gate  leakage  current,  as  shown  in  Figure  1.3. 

Thus,  the  total  input  noise  voltage  can  be  expressed  as 

Sv,»(/)  = ^,(/)  + ^,„(/)|Z,r  + 2^/5v  .„(/)•  ^ Re  [c(Z ) ■ z;  ] (1.6) 

where  Zs  is  the  complex  impedance  and  Cif)  is  the  correlation  coefficient.  Due  to  the 
presence  of  the  gate  leakage  current,  the  total  input  noise  voltage  contains  the  correlated 
noise  component  caused  by  the  common  noise  sources.  These  effects  have  to  be  taken 
into  account  for  accurate  prediction  of  the  noise  performance  in  ultrathin  oxide 
MOSFETs. 


Figure  1.3  Equivalent  input  noise  current  components  in  ultrathin  oxide  MOSFETs 
showing  different  frequency  dependencies. 
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1.3  Organization 

In  this  dissertation,  a detailed  study  of  the  physics-based  compact  noise  models  in 
nanoscale  MOSFETs  will  be  presented,  as  outlined  in  Figure  1.4.  In  Chapter  2,  an 
analytical  model  of  the  gate  leakage  current  in  ultrathin  gate  nitrided  oxide  MOSFETs  is 
presented.  This  model  is  based  on  an  inelastic  trap-assisted  tunneling  (IT AT)  mechanism 
combined  with  a semi-empirical  gate  leakage  current  formulation.  The  tunneling-in  and 
tunneling-out  currents  are  calculated  by  modifying  the  expression  of  the  direct  tunneling 
current  model  of  BSIM  (Berkeley  Short-channel  IGFET  Model).  For  a microscopic 
interpretation  of  the  ITAT  process,  resonant  tunneling  (RT)  through  the  oxide  barrier 
containing  potential  wells  associated  with  the  localized  states  is  proposed.  We  employ  a 
quantum-mechanical  model  to  treat  electronic  transitions  within  the  trap  potential  well. 
The  ITAT  current  model  is  then  quantitatively  consistent  with  the  summation  of  the 
resonant  tunneling  current  components  of  resonant  energy  levels.  Moreover,  the  proposed 
ITAT  process  is  supported  by  the  low  frequency  excess  noise  components.  The  1// noise 
observed  in  the  gate  leakage  current  implies  the  existence  of  slow  processes  with  long 
relaxation  times  in  the  oxide  barrier.  In  order  to  verify  the  proposed  ITAT  current  model, 
an  accurate  method  for  determining  the  device  parameters  is  necessary.  The  oxide 
thickness  and  the  interface  trap  density  of  the  gate  oxide  in  the  20-30  A thickness  range 
are  evaluated  by  the  quasi-static  C-V  method,  dealing  especially  with  quantum- 
mechanical  and  polysilicon  effects. 

Chapter  3 presents  a physics-based  analytical  model  of  the  gate  leakage  current 
noise  in  ultrathin  gate  oxide  MOSFETs.  On  the  basis  of  the  gate  leakage  current 
mechanisms,  we  employ  a barrier  height  fluctuation  model  and  the  Lorentzian  modulated 
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shot  noise  of  the  gate  leakage  current  stemming  from  the  two  dimensional  electron  gas 
channel  to  explain  the  excess  noise  behavior.  The  excess  noise  can  be  interpreted  as  the 
sum  of  1 If  noise  and  the  Lorentzian  modulated  shot  noise.  Trap-related  processes  are  the 
most  likely  cause  of  excess  current  noise  because  slow  traps  in  oxide  can  result  in  low 
frequency  dissipation  in  the  conductance  of  oxides,  and  fast  traps  produce  the  Lorentzian- 
modulated  shot  noise  associated  with  generation-recombination  (G-R)  process  at  higher 
frequencies.  In  order  to  verify  the  proposed  noise  model,  the  simulation  results  are 
compared  with  experimental  data,  and  excellent  agreement  is  observed.  The  noise  model 
can  accurately  predict  the  excess  noise  behavior  of  MOSFETs  with  ultrathin  gate  oxide 
and  is  particularly  suitable  to  be  implemented  in  circuit  simulators. 

In  Chapter  4,  the  electrical  trap  properties  in  ultrathin  oxide  MOSFETs  are 
extracted  using  1 If  low  frequency  noise  components  of  gate  leakage  current.  The 
technique  is  based  on  low  frequency  noise  measurements  which  are  capable  of  probing 
trap  states  over  energy  and  space.  The  gate  voltage  is  swept  between  0.4  and  1.32  V to 
obtain  the  different  power  slopes  y of  the  1 If  noise  spectrum.  From  the  measured  noise 
data,  trap  parameters  such  as  the  slow  trap  densities  and  capture  cross  section  are 
accurately  determined  as  a function  of  trap  energy  and  trap  distance.  The  1 If  noise 
measurement  results  confirm  that  the  capture  cross  sections  depend  strongly  on  the 
surface  potential,  and  the  slow  trap  densities  estimated  near  the  band  edge  and  at  the  Si- 
SiC>2  interface  are  higher  than  those  near  the  midgap  and  at  deeper  distance.  Moreover, 
the  parallel  conductance  spectroscopy  of  slow  traps  is  reproduced  by  1 If  low  frequency 
noise  measurements,  and  is  verified  by  the  low  frequency  plateau  in  measured 
conductance  loss  Gp(f)/w.  For  the  purpose  of  investigating  the  dependency  of  trap 
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parameters  on  the  complex  admittance  of  oxide  traps,  the  Cole-Cole  plots  for  different 
trap  profiles,  a uniform  and  exponentially  decreasing  trap  distribution  in  space,  are 
calculated  using  the  theoretical  model.  The  noise  spectroscopy  technique  further  reveals 
the  dynamic  evolution  of  trap  states  and  gives  a clear  physical  understanding  of  the 
pertinent  low  frequency  noise  sources. 

In  Chapter  5,  an  analytical  model  and  analysis  of  the  low  frequency  noise 
performance  of  ultrathin  oxide  MOSFETs  are  presented.  The  model  includes  the  effects 
of  gate  tunneling  leakage  current,  low  frequency  noise  sources,  and  their  cross  correlation 
coefficient  between  the  input  and  the  output.  The  correlation  coefficient  is  not  negligible 
in  the  existence  of  direct  tunneling  current  through  the  gate  oxide,  and  has  an  impact  on 
the  overall  noise  performance  at  low  frequencies.  Based  on  the  partition  noise  theory  and 
the  BSIM4  gate  leakage  current  model  with  the  source-drain  partition,  the  correlation 
coefficient  is  obtained  from  a physical  and  rigorous  modeling  of  intrinsic  noise  sources. 
The  calculated  results  are  compared  with  correlation  noise  measurements  and  good 
agreement  is  observed.  An  analytical  formulation  of  low  frequency  noise  performance  is 
newly  derived  using  the  frequency-dependent  noise  parameters  P LF,  RLF,  and  &F . Noise 
performance  parameters,  such  as  minimum  noise  figure,  optimum  source  resistance  and 
reactance,  are  graphically  presented  to  accurately  predict  the  influences  of  the  gate 
leakage  current  on  the  low  frequency  noise  performance.  From  quantitative  results 
simulated  with  typical  device  parameters,  it  is  shown  that  the  proposed  model  can  be 
used  for  optimum  designs  of  ultrathin  oxide  MOS  devices  and  circuits  at  low  frequencies. 

Chapter  6 presents  a 1 If  drain  current  noise  model  for  deep-submicron  MOSFETs 
with  ultrathin  oxide.  Based  on  the  number  and  correlated  mobility  fluctuation 
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mechanisms,  the  model  is  derived  incorporating  a tunneling  assisted-thermally  activated 
process  and  a realistic  trap  distribution  inside  the  gate  oxide  layer.  The  effects  of  the 
device  structure  and  processing  technologies  on  the  noise  characteristics  are  taken  into 
consideration  through  a quadratic  mobility  degradation  factor,  a parasitic  resistance,  a 
doping  profile,  and  trap-related  parameters.  For  ultrathin  oxide  MOSFETs,  the  trapping 
efficiency  ratio  and  the  scattering  rate  are  expressed  in  terms  of  the  trap  distance  and  the 
inversion  carrier  density,  enabling  an  accurate  prediction  of  the  noise  behavior.  The 
model  parameters  are  extracted  and  determined  so  that  the  model  has  an  excellent 
agreement  with  the  measured  data  for  advanced  MOSFETs  operating  in  both  the  linear 
and  saturation  regions.  From  quantitative  results  simulated  with  the  extracted  data,  it  is 
shown  that  the  new  model  is  applicable  to  the  design  of  future  CMOS  devices  and  new 
device  processing  technologies,  and  is  suitable  to  be  implemented  in  circuit  simulators. 

In  Chapter  7,  a charge-based  quantum  correction  model  for  current  and  noise  in 
nanoscale  MOSFETs  is  presented.  The  model  is  analytically  formulated  by  incorporating 
quantum  mechanical  (QM)  effects  in  the  inversion  layer  into  the  current  and  noise 
models  of  BSIM.  For  efficient  computation,  the  quantum  correction  models  are 
calculated  for  the  lowest  quantized  energy  level,  which  is  a reasonable  approximation. 
The  QM  effects  on  the  inversion  charge  parameters,  such  as  the  surface  potential,  the 
threshold  voltage,  the  inversion  charge  density,  and  the  gate  capacitance  model,  are  taken 
into  account.  On  the  basis  of  the  quantum  correction  model  of  the  oxide  voltage  and 
oxide  thickness,  the  gate  leakage  tunneling  current  and  1 If  noise  model  are  evaluated  for 
ultrathin  gate  oxide  MOS  devices.  The  quantum  models  of  the  drain  current  and  noise  are 
obtained  by  accounting  for  the  QM  effects  through  the  inversion  charge  density  model. 


16 


Quantitatively,  it  is  confirmed  that  the  quantum  correction  of  the  oxide  thickness  is 
essential  to  the  accurate  prediction  of  the  gate  leakage  current  and  noise  model.  For  the 
1 If  drain  current  noise  due  to  the  number  fluctuation  components,  the  classical  and  QM 
models  behave  in  a qualitatively  similar  fashion  to  the  RTS  noise  in  both  the  weak  and 
strong  inversion,  showing  a noticeable  deviation  in  the  saturation  region.  It  is  also  shown 
that  QM  effects  in  the  inversion  layer  are  important  to  the  thermal  noise  modeling  and, 
for  the  accurate  noise  modeling  of  nanoscale  MOSFETs,  become  more  pronounced  for 
higher  doping  concentration  and  thinner  oxide  thickness. 

In  Chapter  8,  the  comprehensive  and  physics-based  compact  noise  models  for 
advanced  CMOS  devices  are  presented.  The  models  incorporate  all  the  important 
physical  effects  of  nanoscale  MOSFETs,  and  are  formulated  to  be  suitable  for  circuit 
simulators.  Especially,  the  implementation  of  the  gate  leakage  current  and  noise  models 
in  BSIM4  is  demonstrated. 

Finally,  Chapter  9 summarizes  the  primary  contributions  of  this  work  and  then 
concludes  with  suggestions  to  further  advance  noise  models  to  design  future  CMOS 


devices. 
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Figure  1 .4  Organization  and  overview  of  this  dissertation 


CHAPTER  2 

GATE  LEAKAGE  CURRENT  MODEL  AND  ANALYSIS 

2.1  Introduction 

Advancement  of  VLSI  technologies  has  resulted  in  continuous  scaling  down  of 
MOS  devices,  whose  feature  size  used  in  IC  chips  has  been  reduced  to  deep-submicron 
nowadays.  At  the  same  time,  the  continuous  reduction  of  the  gate  oxide  layer  thickness 
in  advanced  CMOS  devices  leads  to  excessive  gate  leakage  currents  and  reliability 
problems  [Ner95,  Sch94a,  Sch94b].  For  such  ultrathin  oxide  layers,  the  gate  tunneling 
current  is  one  of  the  limiting  factors  for  future  MOS  technologies  because  a significant 
tunneling  current  is  expected  even  at  normal  operating  conditions.  This  large  tunneling 
current  increases  stand-by  power  consumption  and  deteriorates  device  performance 
[GheOOa,  GheOOb].  Moreover,  the  gate  leakage  current  in  ultrathin  gate  oxide  MOS 
devices  should  be  considered  as  a low  frequency  noise  source  [Mom98]. 

Nitridation  treatment  of  thin  oxides  is  used  in  MOS  technology  to  improve  the 
dielectric  robustness  to  electrical  stress  and  to  develop  high  permittivity  gate  dielectric 
stacks.  However,  one  of  the  problems  associated  with  the  nitridation  process  is  that  it 
introduces  a large  number  of  electron  traps  into  the  insulating  layer  [Che88].  It  was 
shown  that  the  enhanced  current  in  nitrided  oxide  and  high  dielectric  stacks  are  mainly 
due  to  two-step  trap-assisted  tunneling  [Che88,  Fle92,  Fle93,  Hou99,  Yan96].  Since 
energy  loss  has  been  observed  during  tunneling  through  gate  oxides,  an  inelastic  process 
has  to  be  incorporated  in  these  trap-assisted  tunneling  models  [Tak99].  Endoh  et  al. 
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proposed  an  analytical  ITAT  model  for  stress-induced  leakage  current  (SILC)  in  6.8  nm 
thick  silicon  dioxide  films,  describing  the  quantitative  analysis  of  trap-related  properties 
[End99].  However,  a conclusive  ITAT  current  model  in  the  direct  tunneling  regime  at 
low  voltages  has  not  been  established  [GheOOb]. 

In  order  to  provide  a better  physical  understanding  of  ITAT  mechanisms,  a 
microscopic  interpretation  is  necessary.  Both  Coulombic  and  neutral  trapping  sites  in 
SiC>2  used  in  MOS  technologies  have  been  found  to  be  responsible  for  leakage  current 
and  noise  mechanisms  [Buc91,  Mom98].  Cheung  et  al.  showed  that  tunneling  is 
enhanced  by  a localized  point  defect  and  resonant  tunneling  may  be  the  missing  physical 
picture  [Che97].  If  we  treat  the  trapping  states  as  the  discrete  levels  of  the  quantum  wells 
associated  with  the  defects  [Buc91],  a theory  of  resonant  tunneling  through  localized 
trapping  sites  can  be  applied  to  the  ITAT  current  modeling  [Bro92,  Buc88,  Ric84],  In 
addition,  it  has  been  shown  that  the  inelastic  scattering  process  involved  in  the  ITAT  is 
producing  low  frequency  noise  [Ma92]. 

Accurate  and  efficient  methods  for  determining  the  gate  oxide  thickness  and  trap 
density  are  indispensable  for  device  modeling  and  process  monitoring.  The  gate  oxide 
thickness  is  reaching  a range  where  quantum-mechanical  (QM)  effects  should  not  be 
neglected  [CleOO,  Dmo97,  OmuOO].  QM  effects  mainly  affect  the  energetic  and  spatial 
carrier  distribution  in  inversion  layers,  which  induce  significant  modifications  in  the  CV 
characteristics  and  modeling.  Moreover,  the  polysilicon  gate  effect,  which  results  in  the 
creation  of  a space  charge  in  the  gate  region,  results  in  anomalous  CV  characteristics 
[QudOO].  All  these  effects  have  to  be  taken  into  account  and  evaluated  for  the  accurate 


extraction  of  device  parameters. 
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In  this  Chapter,  we  will  present  a detailed  study  of  a physics-based  compact 
model  for  the  gate  leakage  current  in  ultrathin  gate  nitrided  oxide  MOSFETs.  In  Section 
2.2,  an  ITAT  model  is  proposed  to  characterize  the  observed  conduction  enhancement. 


direct  tunneling  current  step.  In  Section  2.3,  resonant  tunneling  through  the  oxide  barrier 
containing  potential  wells  associated  with  the  electronic  defect  states  is  proposed  as  a 
microscopic  interpretation  of  the  ITAT  process.  Section  2.4  shows  how  the 
characterization  and  modeling  procedure  can  be  applied  to  devices  featuring  ultrathin 
nitrided  oxides  for  the  purpose  of  validating  the  leakage  tunneling  current  model  from  C- 
V,  I-V,  and  noise  measurements.  Finally,  Section  2.5  presents  some  conclusions. 


2.2.1  Direct  Tunneling  Gate  Leakage  Current 

When  the  oxide  layer  is  thicker  than  6 nm  and  the  applied  field  is  in  the  Fowler- 
Nordheim  (FN)  region,  the  leakage  current  can  be  explained  by  FN  electron  tunneling  in 
MOS  structures  [Dep95].  For  an  ultrathin  oxide  layer  and  an  applied  voltage  below  the 
potential  barrier  height,  the  classical  FN  formula  is  not  applicable.  From  an  independent 
electron  approximation  and  for  an  elastic  tunneling  process,  the  direct  tunneling  current 
density  from  the  Si  inversion  layer  to  the  poly-Si  gate  can  be  expressed  as  [Dep95] 


This  model  is  based  on  inelastic  trap-assisted  tunneling  combined  with  a semi-empirical 


2.2  Modeling  of  Gate  Leakage  Current 


(2.1) 


with 


-i2 


(2.2a) 


3/2 


= l-d -VFJ 


(2.2b) 


21 


(2.2c) 
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where  g is  the  electron  charge,  0*  is  the  energy  barrier  height,  h is  the  reduced  Planck 


constant,  tox  is  the  oxide  thickness,  mox  is  the  electron  effective  mass  of  the  conduction 
band  in  the  oxide,  and  Fox  is  the  electric  field  in  the  oxide  layer.  In  order  to  calculate  the 
gate  tunneling  currents  through  thin  oxides,  it  is  necessary  to  relate  the  oxide  voltage  Vox 
to  the  applied  gate  voltage  Vg,  since  Vox  depends  on  Vg  as  well  as  the  surface  potential  tf/s 
and  the  polysilicon  voltage  drop  as  follows 


Vox=Vg-VFB-Vs-Vpo,y  (2.3) 

An  accurate  computation  of  Vox  for  a given  Vg  requires  accounting  for  the  voltage 
drop  across  the  polysilicon  gate  due  to  the  formation  of  a depletion  layer  near  the  poly- 
SiC>2  interface.  For  a positive  voltage  applied  to  an  n+  polysilicon  gate  with  the  p- 
substrate  grounded,  the  polysilicon  will  deplete.  The  value  of  Vpoiy  is  given  by  [QudOO] 
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2 mVg  +l-yjuAmVt 
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where  Npoiy  is  the  polysilicon  doping  concentration  and  Vfb  is  the  flat-band  voltage. 

The  BSIM4  model  can  predict  all  the  significant  direct  tunneling  current 
components  which  are  electron  tunneling  from  the  conduction  band  in  the  substrate  to  the 
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gate,  electron  tunneling  from  the  valence  band  in  the  substrate  to  the  gate,  and  hole 
tunneling  from  the  valence  band  in  the  gate  to  the  substrate.  The  BSIM4  gate  oxide 
tunneling  current  model  is  given  by  [LeeOO,  YeoOO] 


DT.BS1M 


8 xh<t>b£0 


C (</>„, toxFox,VK  )-exp 
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3 hq\Fox\ 


(2.6) 


This  formula  is  very  similar  to  Equation  (2.1)  except  for  the  correction  function  C 
( Vg ,V0x,tox,<Pb)  as  expressed  by  [LeeOO,  YeoOO] 
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where  A is  a fitting  parameter  for  three  current  components,  (pbo  is  the  Si-SiCL  barrier 
height  ( 3.1  eV  for  electrons  and  4.5  eV  for  holes  ),  (pb  is  the  tunneling  barrier  height,  and 
the  function  N is  an  indicator  of  the  density  of  either  tunneling  carrier. 

2.2.2  Inelastic  Trap-Assisted  Tunneling  (ITAT)  Model 

We  formulate  an  improved  model  to  account  for  the  gate  leakage  current  of 
MOSFETs  with  ultrathin  gate  nitrided  oxides.  This  model  combines  the  inelastic  trap- 
assisted  tunneling  mechanism  with  the  semi-empirical  gate  leakage  current  model  of 
BSIM  described  above.  The  inelastic  trap-assisted  tunneling  component  describes 
electrons  tunneling  to  deep-lying  trap  states  and,  immediately,  released  to  positions 
deeper  than  the  original  position  and  subsequently  tunneling  to  the  gate  under  the 
influence  of  the  applied  electric  field  [Tak99].  In  this  process,  the  tunneling-in  current 
from  the  inversion  layer  to  the  traps  and  the  tunneling-out  current  from  the  traps  to  the 
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gate  are  calculated  by  modifying  the  formulation  of  direct  tunneling  in  the  BSIM  model. 
Figure  2.1  shows  the  schematic  band  diagram  of  two-step  inelastic  trap-assisted 
tunneling.  In  the  capture  tunneling  step,  the  electrons  undergo  a direct  tunneling  process 
into  the  trap  sites  with  a tunneling  probability  Tdtsi ■ They  reach  the  gate  via  an  additional 
direct  tunneling  step  with  a probability  W The  tunneling  probabilities  Tdtsi  and  Tortg 
are  related  to  the  tunneling-in  current  7Drm  from  the  inversion  layer  to  the  trap  states  and 
to  the  tunneling-out  current  Jdtoui  from  the  trap  states  to  the  gate.  Assuming  that  x is  the 
distance  from  the  Si-Si02  interface,  Ntrap{x,E)  is  the  sheet  trap  density  in  cm'2  at  a 
distance  * and  energy  E,  f(x,E)  is  the  electron  occupancy  of  the  traps,  Ag  is  the  gate 
tunneling  area,  and  cr,  is  the  capture  cross  section  of  the  traps,  the  current  densities  Jam 
and  Jdtoui  are  then  expressed  as 


where  <p,  is  the  barrier  height  of  the  oxide  trap  states,  and  g(0b^,Eox)  and  r((pb+x,Fox)  are 
the  uniform  generation  and  recombination  rate,  respectively. 

The  resulting  tunneling  current  Jjtat  of  this  trap-assisted  tunneling  process  is 
given  by  a detailed  balance  of  Jorn  and  Jdtoui ■ For  an  electron  to  contribute  to  the 
leakage  current,  it  must  first  tunnel  from  the  substrate  into  an  intermediate  trap,  then 
tunnel  out  of  the  trap  towards  the  gate,  events  which  are  characterized  by  the  tunneling 
rates  gs,  and  rtg.  Since  both  events  must  occur  in  sequence,  the  total  tunneling  rate  of  an 
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Figure  2.1  Schematic  band  diagram  of  two-step  inelastic  trap-assisted  tunneling  through 
the  gate  oxide 


electron  tunneling  through  the  entire  oxide  via  a specific  trap  site  at  a location  x within 
the  oxide  can  be  expressed  as  [End99] 


R(E  x F ) - 8Wb’*’Foxl)-r(<l>b-qFoxlx  + Eloss,tox-x,Fox2) 
g{(t>b , x,  FoxX ) + r(<t>b  - qFoxlx  + 

E loss  » ^ ox  X’  ^ oxl  ) 

and 

JmAE,x,FJ  = -S-o,N,r,p(x,qF„x-0l,)-RI.E,x,F„)  (2,10) 

A, 

The  local  electric  fields,  Fox\  and  Fox 2,  characterizing  the  separated  tunneling  regions  can 
be  obtained  by  considering  the  equivalent  MOS  capacitor  circuit  and  using  Gauss’s  law 


[End99] 
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If  a parabolic  dispersion  relation  is  used  for  the  transversal  energy  component  E± 
with  a transversal  electron  effective  mass  m±,  the  uniform  generation  and  recombination 
rates  can  be  represented  by  [Dep95] 

4 mn,  A [E/s  ^ 

8 (0b\  > x*  Kx\ ) = r^<Pb2,x">  Eox 2 ) = ^3  ' Tdt  {(pb\,2  ~ E + E fs , x,  Foxl2 , E± ) dELdE 

(2.12) 

where  h is  the  Planck  constant,  E/s  is  the  Fermi  energy  level  in  the  Si  substrate,  and  E is 
the  total  energy  of  the  tunneling  electron  measured  from  the  Si  conduction  band  edge  at 
the  Si-SiC>2  interface.  The  tunneling  transition  probability  Tdt  for  direct  tunneling  under 
the  WKB  approximation  is  given  by 
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Integrating  Equation  (2.12)  with  a substitution  of  a first  order  Taylor  approximation  of 
Equation  (2.13)  and  introducing  the  correction  function  of  BSIM,  we  obtain  the  analytical 
expression 
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where 
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Putting  the  above  equations  together,  the  total  IT AT  current  density  Juat  is  given  by 
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This  formula  is  a function  of  the  local  electric  fields  Foxl  and  Fox 2,  the  oxide  thickness  tox, 
the  barrier  height  of  the  oxide  trap  states  4>, , and  the  barrier  height  between  the  substrate 
and  oxide  (pb ■ Note  that  this  current  flow  through  the  oxide  layer  is  valid  for  Vox\  < <pb- 
2.3  Microscopic  Analysis  of  Current  Transport 
Resonant  tunneling  (RT)  through  the  oxide  barrier  containing  potential  wells 
associated  with  the  electronic  states  is  proposed  as  a microscopic  interpretation  of  the 
inelastic  trap-assisted-tunneling  process.  In  order  to  avoid  the  difficulty  of  treating 


27 


electrons  classically  when  they  are  confined  only  a few  nanometers  away  from  the 
trapping  sites,  we  can  employ  a quantum-mechanical  model  which  describes  the 
electronic  states  in  the  potential  well  [Buc91].  In  this  model,  the  electrons  are  bound  to 
localized  states,  and  electronic  transitions  from  any  localized  state  to  any  other  occur 
within  the  trap  from  interactions  of  the  electrons  with  phonons  and  photons. 

If  the  trap  state  has  an  energy  broadened  level  due  to  its  interaction  with  the 
emitter  states  lying  near  the  Fermi  energy,  then  resonance  tunneling  can  occur  through  a 
new  channel  where  the  tunneling  electron  goes  from  the  emitter  to  the  trap  state,  and 
subsequently  tunnels  through  the  reduced  barrier  between  the  trap  state  and  collector 
[Ban88,  BlaOO,  C0086].  This  resonance  tunneling  process  causes  an  enhancement  in  the 
tunneling  current  density  when  the  electron  energy  and  the  trap  energy  level  line  up. 
However,  tunneling  electrons  suffer  inelastic  events  in  the  potential  well.  Inelastic  events 
lead  to  a broadening  and  decrease  in  the  resonant  transmission.  The  decrease  of  the  peak 
transmission  with  increasing  inelastic  scattering  is  accompanied  by  an  increase  in  the  off- 
resonant  transmission.  Therefore,  the  theory  of  resonant  tunneling  through  isolated  trap 
states  should  be  modified  to  include  inelastic  scattering  processes  occurring  in  the 
trapping  state  [BlaOO,  Bro92,  But88]. 

In  the  neighborhood  of  a resonance,  the  total  tunneling  probability  for  an 
electron  to  traverse  through  a double  barrier  with  a quantum  well  is  [But88] 
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where  Enres  is  the  nth  resonance  energy  level,  j,,  is  a peak  value  at  nth  resonance 
energy,  T„  is  the  sum  of  all  the  partial  widths,  IV  is  the  total  elastic  width,  T,/  is  the  total 
inelastic  width,  YnDTs'  and  VnDT,g  are  the  partial  elastic  width  of  the  resonant  level,  v„  is 
the  attempt  frequency  of  the  nth  resonant  level,  and  p is  the  probability  for  a carrier 
approaching  the  trapping  state  to  be  scattered  into  the  lower  state.  The  widths  of  the 
resonance  are  a function  of  the  applied  voltage,  and  are  explicitly  affected  by  the  sheet 
trap  density.  The  formula  for  the  resonant  tunneling  current  density  is  then  [BlaOO] 
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where  Nres  is  the  number  of  resonant  energy  levels.  On  integrating,  the  current  density  is 
given  by 
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Within  the  resonant  tunneling  current  region,  (£/s-£’„re's)»r«  and  Enresy>Yn,  and  the 
current  density  is  then  approximated  by 
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This  formula  compares  directly  to  the  ITAT  current  density,  Equation  (16),  and  shows 
that  the  total  current  density  is  expressed  as  the  summation  of  the  resonant  tunneling 
components  of  resonant  energy  levels. 

2.4  Experimental  Verification 

We  illustrate  how  the  characterization  and  modeling  procedure  can  be  applied  to 
devices  featuring  ultrathin  oxides  for  the  purpose  of  validating  the  leakage  tunneling 
current  model  using  C-V,  I-V,  and  noise  measurements. 

2.4.1  Characterization  of  Ultrathin  Oxides 

MOS  capacitors  were  fabricated  using  a manufacturable  remote  plasma  nitrided 
oxide  (RPNO)  process  [Hat96,  NicOO].  About  1.7-4. 3 run  of  gate  oxide  was  grown  with 
both  furnace  and  rapid-thermal  based  processes.  Nitridation  was  performed  at  room- 
temperature  by  exposing  the  gate  oxide  to  a short,  high-density,  remote  helicon-based 
nitrogen  discharge.  This  was  then  followed  by  post-nitridation  annealing  in  an  inert  or 
low  partial-pressure  oxygen  ambient.  MOS  capacitors  were  fabricated  with  standard 
processing  techniques  on  both  p and  n-Si  substrate  with  n+  and  p+  poly-Si  gates,  using 
LOCOS  isolation  to  define  the  active  regions.  Following  N2-plasma  nitridation  and  post- 
nitridation  anneal  of  the  gate  oxide,  a 270  nm  undoped  poly-Si  gate  layer  was  deposited. 
For  n+  gates,  P was  implanted  (5.0xl015  cm'2,  45  keV).  Activation  was  performed  with  a 
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850  °C  / 30  min  anneal.  For  p+  gates,  the  poly-Si  was  implanted  with  different  doses  of  B 

(1.5xl015-4.5xl015  cm'2)  at  10  keV  energy  in  three  of  the  four  quadrants.  Rapid-thermal 
annealing  was  then  performed  at  1000  °C  in  Nt  for  different  durations. 

Figure  2.2  depicts  a schematic  energy  band  diagram  that  was  simulated  by 
running  the  program,  SCHRED_2.0.  This  illustrates  a n+  polysilicon-Si02-p-Si  MOS 
structure  in  inversion  showing  the  quantization  effects  of  the  substrate  electron  energy. 
Ec  is  the  Si  conduction  band  edge  and  EfS  is  the  Fermi  level  in  the  Si  substrate.  Ey  is  the 
energy  level  of  the  jth  subband  in  the  ith  valley,  and  79.1%  and  16.6%  are  the  carrier 
occupancy  for  each  energy  subband  in  each  valley  of  E-k  space. 

QM  and  polysilicon  effects  induce  significant  modifications  in  the  CV 
characteristics,  modeling,  and  parameter  extraction  procedure  [CleOO,  Dmo97,  OmuOO]. 
For  accurate  determination  of  the  gate  oxide  thickness  and  interface  trap  density,  all  these 
effects  have  to  be  taken  into  account.  Assuming  that  most  electrons  occupy  the  ground 
state  of  the  subband,  we  can  employ  the  following  approximate  expression  for  the  total 
capacitance  Ctotai- 
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Figure  2.2  Schematic  band  diagram  of  a n+  polysilicon-SiCh-p-Si  MOS  structure 
showing  the  quantization  effects  of  the  substrate  electron  energy 
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where  Cox  is  the  oxide  capacitance,  C.s,-  is  the  silicon  capacitance,  Cpoiy  is  the  capacitance 
of  the  polysilicon  depletion  region,  Vpoiy  is  the  voltage  drop  in  the  polysilicon  depletion 
region,  i//s  is  the  surface  potential,  Npoiy  is  the  doping  concentration  in  the  polysilicon 
gate,  Ndep  is  the  depletion  charge  density,  Co  is  a constant,  zav  is  the  centroid  position  of 
mobile  carriers  in  the  inversion  or  accumulation  layer,  and  m is  the  effective  electron 
mass.  For  convenience,  we  can  use  the  following  approximate  semi-classical  expression 


for  Qltys)  [OmuOO]. 
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where  Lp  is  the  Debye  length  and  F(qipJkT)  is  a semi-classical  function  of  the  surface 
potential.  In  Equation  (2.23a),  Co  is  introduced  to  make  a correction  of  effective  width  of 
the  inversion  layer.  From  self-consistent  simulation  results  for  Co  [OmuOO],  an  empirical 
and  useful  expression  can  be  obtained  for  the  (100)  silicon  surface  as 


C0  =0.41  + exp 
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where  Nm  is  a constant  ( =1.264xlOn  cm"2).  The  gate  voltage  applied  to  polysilicon- 
SiCh-p-Si  MOS  structure  is  described  by 

y,=y»+^+K-^-  <i27> 

^ox 

The  interface  trap  density  is  evaluated  from  the  quasi-static  C-V  technique. 
Because  the  interface  trap  capacitance  C„  is  in  parallel  with  the  silicon  capacitance,  the 
following  equation  is  used  for  extracting  the  interface  trap  density  Nit(  cm"2eV_1 ). 
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(a) 


(b) 

Figure  2.3  Comparison  of  C-V  and  interface  trap  density  evaluated  by  the  proposed 
model,  classical  model,  and  measured  data,  (a)  C-V  characteristics  and  (b) 
evaluation  of  interface  trap  density  for  n+poly-Si  nMOSFET  with  tsiON  = 3.21 
nm  ( tOXjeq  = 2.2  nm  ),  E siON  = 5.7eo,  Nsub  = 6.0xl017  cm"3,  Npoiy  = 1.85xl020 
cm"3,  and  Ag  = 140x70x4  gm2  at  T = 300  K 
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where  Clolai,m  is  the  measured  total  capacitance  and  C„  is  the  semi-quantum-mechanical 
semiconductor  surface  capacitance  including  the  inversion  layer  capacitance.  The 
interface  trap  density  is  derived  from  Equation  (2.28)  as  follows 


In  order  to  verify  the  proposed  model,  the  simulation  results  of  the  C-V 
characteristics  by  using  Equations  (2.22)-(2.23)  is  shown  in  Figure  2.3  (a)  with  the 
classical  simulation  results  and  the  measured  data  for  n+poly-Si  nMOSFET.  The 
equivalent  oxide  thickness  tox,eq  can  be  calculated  as  t0x,eq=(tJtsiON)  tsiON  [Vog98], 
where  the  nitrided  oxide  and  oxide  dielectric  constant,  Esion  and  eG*,  are  assumed  to  be  5.7 
Eo  and  3.9  8o,  respectively.  The  C-V  characteristics  simulated  by  the  proposed  model 
agree  well  with  the  measured  data.  Figure  2.3  (b)  shows  the  interface  trap  density 
evaluated  by  using  Equation  (2.29)  and  the  classical  model.  Near  the  band  edge,  the 
difference  between  two  models  increases  drastically.  Consequently,  when  the  quantum 
mechanical  influence  and  the  polysilicon  depletion  effect  are  anticipated  to  be  significant, 
the  classical  model  cannot  be  applied  to  extract  the  interface  trap  density  from  C-V 
measured  data. 

2.4.2  Current-Voltage  Characteristics 
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The  gate  leakage  I-V  characteristics  were  simulated  using  the  proposed  ITAT 
model  of  Equation  (2.16).  The  I-V  characteristics  are  dependent  on  the  unknown  trap 
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density  Nlrap(x,E).  It  is  assumed  in  our  calculation  that  Ntrap(x,E)  is  a constant  and  equals 
the  interface  trap  density  Nu,  irrespective  of  the  position  and  the  energy  level.  Figure  2.4 
(a)  shows  the  calculated  I-V  characteristics  of  the  IT AT  model  as  a function  of  the  gate 
voltage  with  the  energy  loss  during  ITAT  process  as  a fitting  parameter.  This  result  is 
compared  with  the  direct  tunneling  model  and  measured  data  for  n+poly-Si  nMOSFET. 
In  order  to  corroborate  the  simulation  results,  we  used  the  values  of  the  oxide  thickness 
and  trap  density  obtained  from  C-V  evaluation  including  the  QM  and  polysilicon 
depletion  effect.  The  simulation  was  carried  out  with  tsioN  = 3.21nm  ( tox,eq  = 2.2nm  ), 
(pb.SiON  = 2.6eV  [Vog98],  mox  = OAm0  [YeoOO],  cr,  = 5.0xl0'14cm2  [GheOOb],  N,rap  = 
2.0xl012  cm'2,  x,=  0.46  tsioN  , (f>t  = 2.5eV,  and  A = 0.6  [YeoOO].  From  the  comparison  of 
the  experimental  results  and  ITAT  theory,  it  was  found  that  the  energy  loss  ranges  from 

0. 1  to  0.6  eV  for  gate  voltages  0.3  - 1.5  V.  It  is  shown  that  the  results  calculated  by 
ITAT  have  better  agreement  with  the  experimental  results  than  those  from  the  direct 
tunneling  model.  This  is  due  to  the  fact  that  the  direct  tunneling  effect  gives  a small 
contribution  to  the  total  tunneling  current.  With  appropriate  changes  of  barrier  height, 
effective  mass,  and  voltage  polarity,  the  models  are  equally  applicable  to  p-channel 
MOSFETs.  The  simulation  was  performed  by  using,  fs,cw=3.21nm  (tox,eq= 2.2nm), 
(pb,siON= 2.5eV,  mox=0.32mo,  cr,=5.0xl0'14cm2,  /V/rap=2.0xl012cm'2,  *,=0.29 tSiON,  0,=2.4eV, 
and  A=0.4.  Again  good  agreement  between  the  measured  data  and  simulation  results  can 
be  observed  in  Figure  2.4  (b).  Figure  2.5  illustrates  the  effect  of  trap-related  parameters, 

1. e.,  trap  position  located  away  from  the  Si-SiC>2  interface  and  trap  energy  level  measured 
from  the  conduction  band  edge  of  the  oxide,  on  the  I-V  characteristics  for  n+  ploy-Si 
nMOSFET.  The  current  peaks  for  traps  located  in  the  middle  of  the  oxide  layer  and 
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Figure  2.4  Comparison  of  I-V  characteristics  simulated  by  the  proposed  model  and  DT 
model,  and  measured  data  with  (a)  tsioN  = 3.21nm,  Esion  = 5.7£o,  (pb.siON  = 
2.6eV,  mox=  0.4m0,  cr,=  5xl0'l4cm2,  N,rap=  2xl012cm'2,  A = 0.6  (b)  tSiON  = 
3.21nm,  E SiON=  5.7e0,  (pb.siON  = 2.5  eV,  mox=  0.32 m0,  cr,=  5xl0‘14cm2,  N,mp  = 
2xl012cm‘2,  A = 0.4 
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Figure  2.5  Illustration  of  the  effect  of  trap-related  parameters  on  I-V  characteristics  for 
n+poly-Si  nMOSFET  (a)  Influence  of  trap  energy  levels  and  trap  positions 
on  current  density  (b)  contour  plot  matching  with  experimental  data 
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Figure  2.6  Comparison  of  I-V  characteristics  calculated  by  the  measured  data  and  the 
proposed  resonant  tunneling  (RT)  model  with  v„  = 1014  s'1,  p = 0.01,  tSiON  = 
3.21nm,  x,=  0.27 tSiON,  mox=  0.4m0,  cr,=  5xl0'14cm2,  Nlrap=  2xl012cm'2 

exponentially  decreases  for  traps  away  from  the  center  of  the  oxide  and  at  deeper  energy 
levels,  as  shown  in  Figure  2.5  (a).  Figure  2.5  (b)  shows  the  contour  plot  of  current 
density  matching  with  experimental  results.  The  trap  position  x,  = 0.46  tsioN  and  energy 
level  0,=2.5eV  extracted  from  a fitting  procedure  are  consistent  with  the  contour  line  of 
trap. 

For  the  purpose  of  understanding  the  ITAT  process  at  the  microscopic  level,  we 
have  calculated  the  resonant  tunneling  (RT)  current  for  several  resonant  trap  energy 
levels,  as  shown  in  Figure  2.6.  For  our  calculations,  we  take  for  the  attempt  frequency  v„ 
=1014  s'1  [Ric98]  . The  RT  current  curves  for  resonant  energy  levels  are  characterized  by 
the  negative  differential  resistance  (NDR)  region  between  the  peak  and  valley.  The 
peak  current  occurs  when  the  energy  of  an  electron  tunneling  through  the  oxide  exactly 
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Figure  2.7  Gate  current  noise  spectral  density  for  n+poly-Si  nMOSFET  with  tstoN  = 
3.21nm  ( tox,eq  = 2.2nm  ),  Nsub=  6.0xl017cm"3,  Npoiy  = 1.85xl020cm'3,  and  A* 
= 140x70x4  /im2  at  T=300K 


corresponds  to  the  trap  energy  level.  Quantitatively,  it  is  found  that  the  ITAT  current  can 
be  reproduced  by  the  summation  of  the  resonant  tunneling  components,  based  on  a 
theoretical  model  which  takes  into  account  two  features,  the  trap  states  serving  as 
quantum  wells  and  trapped  electrons  undergoing  inelastic  events.  This  confirms  that  the 
ITAT  mechanism  is  probably  at  the  origin  of  the  gate  leakage  current  in  ultrathin  nitrided 
oxide  MOSFETs. 

2.4.3  Leakage  Current  Noise  Characteristics 

Figure  2.7  shows  a typical  gate  current  noise  spectrum  measured  at  low  voltages 
for  n+  poly-Si  nMOSFET  with  ultrathin  oxide  fabricated  by  remote  plasma  nitrided  oxide 
(RPNO)  process.  The  noise  spectrum  clearly  shows  a 1 If  frequency  dependence  with 
different  power  factors  y depending  on  the  gate  bias.  The  magnitudes  and  power  factors 
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of  the  low  frequency  current  noise  are  practically  a constant  between  0.4V  and  0.8V. 
However,  when  the  gate  voltage  is  above  0.8V,  the  magnitudes  and  power  factors  of  the 
low  frequency  noise  increase  with  an  increase  in  gate  voltage.  The  bias  dependence  of 
the  1 If  noise  indicates  the  existence  of  inelastic  processes  occurring  in  the  trap  states. 
Moreover,  this  is  in  accordance  with  the  dependence  of  energy  loss  on  the  gate  voltage. 
Physically,  the  1 If  noise  observed  at  low  frequencies  implies  the  existence  of  slow 
processes  with  long  relaxation  times  in  the  oxide  barrier.  Electron-phonon  scattering 
involved  in  this  process  gives  rise  to  low  frequency  dissipation  in  the  conductance  of  the 
oxide  traps.  This  in  turn  results  in  barrier  height  fluctuation,  and  cause  the  fluctuation  in 
the  tunneling  current. 

2.5  Conclusion 

A physical  model  and  analysis  applicable  to  the  low  voltage  gate  leakage  charge 
transport  in  ultrathin  gate  nitrided  oxide  MOSFETs  was  presented.  The  model  is  based 
on  inelastic  trap-assisted  tunneling  (ITAT)  mechanism  combined  with  a semi-empirical 
direct  tunneling  current  model.  The  model  includes  the  effect  of  charge  stored  in  the  trap 
states  and  the  dependence  of  the  energy  loss  on  the  oxide  field  during  tunneling  through 
the  trap  states.  It  was  found  that  the  results  simulated  by  the  ITAT  model  show  good 
agreement  with  the  experimental  results.  To  clarify  the  mechanisms  of  the  ITAT  current 
based  on  a collective  approach,  resonant  tunneling  (RT)  through  barriers  containing 
defect  state  potential  wells  was  formulated.  From  the  comparison  of  the  ITAT  model  and 
RT  formulation,  it  was  found  that  the  origin  of  the  gate  leakage  current  might  be 
attributed  to  the  ITAT  mechanisms  in  ultrathin  nitrided  oxide  MOSFETs.  Furthermore, 
the  ITAT  process  assumption  is  supported  by  the  low  frequency  excess  noise  spectra. 
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For  the  purpose  of  validating  the  leakage  tunneling  current  model,  the  oxide  thickness 
and  the  interface  trap  density  of  gate  oxide  were  determined  by  the  quasi-static  C-V 
method,  taking  quantum  mechanical  effects  and  polysilicon  depletion  into  account . The 
C-V  characteristics  simulated  by  our  model  agree  well  with  the  measured  data. 


CHAPTER  3 

NOISE  MODEL  OF  GATE  LEAKAGE  CURRENT 

3.1  Introduction 

One  of  the  important  problems  associated  with  the  downscaling  of  CMOS 
technologies  is  an  increase  in  the  importance  of  random  noise,  which  is  always  a practical 
limit  to  the  performance  of  a device  [Cel99,  Sim99].  As  a result,  the  signal-to-noise  ratio 
of  analog  circuits  is  reduced,  and  bit  errors  occur  in  data  transmission  systems. 
Especially,  the  scaling  of  device  dimensions,  resulting  in  the  use  of  thinner  oxide  layers, 
is  necessary  for  enhancing  circuit  speed  at  low  voltages,  providing  the  same  capacitance 
within  a smaller  area,  and  offering  lower  programming  voltages  for  non-volatile  memory. 

It  is  generally  accepted  now  that  the  miniaturization  of  devices  will  cause  an 
increase  of  low  frequency  noise  level  due  to  an  observed  scaling  law  [Cel99].  Since  low 
frequency  noise  is  strongly  technology  sensitive,  which  in  some  cases  can  also  be  used  as 
a predictive  or  diagnostic  tool  for  device  lifetime  and  reliability  issues,  the  introduction  of 
new  processing  steps  may  introduce  low  frequency  components  not  presently  observed  or 
described  by  existing  noise  theories  [Sim99].  Moreover,  the  low  frequency  noise  not  only 
determines  the  lower  end  of  the  dynamic  range  of  device  operation  in  the  base  band,  but 
also  has  potentially  a major  impact  on  the  phase  noise  of  nonlinear  circuits  in  the  GHz 
region  through  up-conversion  [KluOO].  In  addition  to  the  channel  current,  gate  leakage 
current  in  ultrathin  gate  oxide  MOS  devices  should  be  considered  as  a low  frequency 
noise  source  [Mom98].  It  is  known  that  the  appearance  of  noise  in  the  gate  current  of 
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MOS  structures  potentially  precedes  oxide  breakdown  [CruOOa,  CruOOb,  Ner95].  The 
gate  tunneling  current  noise  interpreted  at  a microscopic  scale  can  be  a powerful 
characterization  tool  in  evaluating  the  ultrathin  oxide  reliability  [Ner95].  However,  if  we 
can  understand  the  physical  origin  of  the  noise,  we  can  develop  strategies  for  minimizing 
it  or  its  influence.  Therefore,  an  accurate  modeling  and  characterization  of  noise  in  deep- 
submicron  MOSFETs  with  ultrathin  gate  oxide  is  needed  to  design  future  CMOS  devices, 
analog  circuits,  and  new  device  processing  technologies. 

In  this  Chapter,  we  will  present  a detailed  study  of  a physics-based  compact 
noise  model  for  deep-submicron  MOSFETs  with  ultrathin  gate  oxide.  In  Section  3.2,  on 
the  basis  of  the  gate  oxide  current  leakage  mechanisms,  we  employ  a barrier  height 
fluctuation  model  and  a Lorentzian-modulated  shot  noise  model  of  the  gate  leakage 
current  to  explain  the  excess  noise.  The  excess  noise  can  be  interpreted  as  the  sum  of  1 If 
noise  and  enhanced  shot  noise.  The  1 If  noise  component  at  low  frequencies  is  described 
by  thermal  noise  of  the  frequency-dependent  conductance  of  the  oxide  slow  traps  which 
modulates  the  current  tunneling  through  the  oxide  barrier.  The  positive  correlation 
between  tunneling  current  components,  associated  with  generation-recombination  (G-R) 
process  due  to  oxide  fast  traps,  produces  the  Lorentzian-modulated  shot  noise  at  higher 
frequencies,  leading  to  enhanced  shot  noise.  In  Section  3.3,  we  show  experimental 
results,  exhibiting  the  excess  noise  consisting  of  1 If  component  at  low  frequencies  and 
enhanced  shot  noise  at  higher  frequencies.  For  experimental  verification,  the  proposed 
noise  model  is  compared  with  experimental  results,  and  excellent  agreement  is  observed. 
Finally,  Section  3.4  presents  some  conclusions. 
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3.2  Formulation  of  Noise  Model 

We  present  an  analytical  noise  model  for  the  gate  leakage  current  through  a Si02 
barrier  in  the  inelastic  trap-assisted  tunneling  regime.  Trap-related  processes  are  the  most 
likely  cause  of  excess  current  noise  because  traps  in  the  oxide  can  result  in  low  frequency 
dissipation  in  the  conductance  of  oxides  and  produce  the  Lorentzian-modulated  shot 
noise  associated  with  generation-recombination  (G-R)  process  at  higher  frequencies 
[Ale96,  Che99,  Ghi87,  Kum77].  The  gate  current  noise  spectral  density  consists  of  two  or 
three  Lorentzians  in  addition  to  shot  noise  and  1 If  noise. 
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(3.1) 


where  the  first  term  corresponds  to  shot  noise  with  the  Fano  factor  F,  the  second  to 
flicker  noise  or  1 If  noise,  and  the  third  to  G-R  noise  of  discrete  time  constants  r,- . As 
shown  in  Figure  3.1,  one  possible  mechanism  for  1 If  noise  contribution  is  for  traps  to 
affect  locally  the  barrier  height  or  shape  due  to  thermal  noise  of  the  frequency-dependent 
conductance  of  the  oxide  slow  traps  [Ghi87,  IanOO,  Kle82].  This  fluctuation  in  turn 
modulates  the  tunneling  transmission,  and  causes  the  fluctuation  in  the  tunneling  current. 
The  mechanism  for  shot  noise  associated  with  G-R  process  can  be  explained  by  carrier 
transport  based  on  ITAT  [IanOO,  Lee02],  The  generation  rate,  i.e.,  the  transition  rate  from 
the  substrate  and  the  gate  to  the  unoccupied  trap  in  the  oxide,  and  recombination  rate,  i.e., 
the  transition  rate  from  the  occupied  trap  to  the  substrate  and  the  gate,  are  caused  by  the 
statistical  occupancy  fluctuation  of  fast  traps  in  the  oxide. 
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Figure  3.1  Flow  diagram  of  the  physics-based  processes  for  the  noise  model  calculations 
in  the  low  and  high  frequencies 
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The  fluctuations  in  the  generation  and  recombination  rates  in  turn  modulate  the  tunneling 
current  through  traps  in  the  oxide,  as  illustrated  in  Figure  3.1.  From  the  trap-assisted 
tunneling  associated  with  G-R  process,  therefore,  we  can  obtain  an  expression  consisting 
of  full  shot  noise  and  G-R  noise,  which  is  referred  to  as  the  Lorentzian-modulated  shot 
noise  [Vli81].  The  observed  enhancement  of  noise  at  higher  frequencies  is  probably  due 
to  G-R  noise  contributed  by  some  positive  correlation  [Ian97,  Ian98,  Rek99,  RekOO, 
Vas02]. 

3.2.1  Barrier  Height  Fluctuation  Model 

The  noise  equivalent  circuit  of  an  ultrathin  oxide  MOSFET  with  the  source  and 
drain  tied  together  is  described  in  Figure  3.2  [Eat72,  Oga97,  Vik95,  VogOO].  Gt,dc  is  the 
frequency-independent  dc  tunneling  conductance,  Rt  is  the  inversion  layer  resistance, 
Gp(f)  is  the  frequency-dependent  parallel  conductance  of  oxide  traps  in  n'eV'cm'2,  and 
Cp{  f)  is  the  frequency-dependent  parallel  capacitance  of  oxide  traps  in  Fcm'2.  The  1 If 
noise  might  be  originated  by  thermal  noise,  4kTGp(  f ),  due  to  the  frequency-dependent 
parallel  conductance  Gp(f),  producing  the  barrier  height  fluctuation.  In  terms  of  Gp(f)/w 
loss  observed  at  low  frequencies,  the  slow  traps  contributing  to  the  low  frequency  plateau 
in  Gp(f)/w  are  responsible  for  the  measured  1 If  noise  [Ghi87,  Kum77,  Kir89,  Nic67a, 
Nic67b,  Ure89,  Ure95], 

On  the  basis  of  the  Wentzel-Kramers-Brillouin  (WKB)  approximation,  the 
tunneling  transmission  coefficient  for  electrons  can  be  expressed  as  [Dep95] 


Wi)  = ex  p 


(3.2) 
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Gt,dc 


Figure  3.2  Noise  equivalent  circuit  of  an  ultrathin  oxide  MOSFET  with  the  source  and 
drain  tied  together,  corresponding  to  the  oxide  trap  occupancy  controlled  by 
semiconductor  electrons. 

where  the  classical  turning  point  xi  for  direct  tunneling  is  equal  to  tox.  Due  to  thermal 
noise  of  the  equivalent  conductance  modeling  oxide  trap  occupancy  fluctuations,  the 
barrier  height  0/,  fluctuates  as  follows 


</>b(x,t)  = <pb  +S<f)h(x,t)  (3.3) 


where  <50b  is  a small  fluctuation  in  <f>b{x,t)  which  is  a function  of  position  x and  time  t. 
The  fluctuation  in  the  tunneling  transmission  is  given  by 


Sr,(x,t) 


■S<t>b(x,t) 


(3.4) 
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Using  Equations  (3.2),  (3.3),  and  (3.4),  we  can  obtain  the  fluctuations  in  the  tunneling 
transmission  T,  as  follows  [Kle78,  Kle82] 


S<pb{x,t) 


:dX 


+ E,1-E  + ^El-q'^ 
m„r  t„. 


(3.5) 


Replacing  xi  = tox  for  direct  tunneling  and  making  a Fourier  transform,  we  find  the 
spectral  noise  density  of  the  transmission  to  be 


St,  if)  = T,2 


2 m„ 


ft 


SAx\x\f) 


m L „ V„, 


</>b  + Efs~E  + EL-q^x 

moX  Kx  J 


"V  m i/  „ 

0b+EJk-E  + ^E1-q-*-x 


-dxdx 


(3.6) 


* ox  y 


where  S^(;c  ,*  ,/)is  the  spectral  cross-correlation  density  in  thermal  noise  at  positions 
jc  and  ;c  , and  is  given  by  [Kle78,  Kle82,  Kum77,  Vik95] 


nglox 


S{xa-xb)~ 


1 


dxadxb 


(3.7a) 


with 


Gpif) 


R e[Z(/)]  = 


1 

27? 

2 7f 

~GP(J) 
2 Jtf 

2 

+[c„+c,(/)]! 

(3.7b) 


where  k is  Boltzmann’s  constant,  T is  the  temperature,  Cox  is  the  oxide  capacitance  in 
Fern'2,  6 is  the  Dirac  delta  function,  and  Re[Z(/)]  is  the  real  part  of  frequency-dependent 
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impedance  between  x = 0 and  x = tox  in  Figure  3.2.  Using  Equations  (3.6)  and  (3.7),  the 
spectral  noise  density  for  the  tunneling  transmission  is  found  to  be 


M-\tox  f 


Gp(f ) 
2 


Gp(f ) 
2 


- 2 

+k.+c,(/)]! 


ft 


If 


S(xa  -xb)-y-  dxadxh 


m, 


</>b+Efs-E  + ^E±-q-^x 


m„ 


-dxdx 


Til  V 

</>b  + Efs-E  + -^E±-q^x 


lox  y v v 


1 ox  y 


(3.8) 


On  integrating  Equation  (3.8),  we  get  the  spectral  noise  density  as  a function  of  the  oxide 
voltage  ( Vox  ^ <t>b  ) 


Sr  (/) 


Gp(f ) 
2# 


Gp(f) 

2/ f 


+[c„+c„(/)| 


b / 


(3.9) 


with 


16  moxkTq2 

w2\t2oxxi 


(3.10a) 


B(VJ  = ^X,(X,-qVj[sX,(X, -qVJ-3qX,}-2[2X!(2X,  -3 qVJ  + q^i]  (3.10b) 


Til  1/ 

X{=4>b  + Efs-E  + -±Et,  X2=q 

m„r  t„ 


(3.10c) 


where  Gp(  f )/2 nf  and  Cp{  f ) are  independent  of  frequency  at  low  frequencies.  The 
fluctuations  in  the  tunneling  current  are  generated  by  the  transmission  fluctuations. 
Assuming  that  a single  transmission  probability  applies  to  all  transitions  in  calculating  the 


50 


tunneling  current  [Gar83,  Vog98],  the  spectral  density  of  the  current  fluctuations  is  given 
by 


d In  7(7) ) 

dlnT. 


-i2 


Sr,(f)_.2 

'T  2 1 g 


ST,(f ) 


(3.11) 


The  oxide  traps  are  not  only  continuously  distributed  in  energy  but  are  also 
spatially  distributed  in  the  oxide.  In  the  case  of  direct  tunneling  of  carriers  in  the 
conduction  or  valence  band  into  traps  located  at  some  energy  and  at  some  distance  x 
away  from  the  interface,  the  parallel  conductance  Gp{  f ) and  capacitance  Cp(  f ) are 
expressed  as  [Bal78,  Ghi87] 


Nr(E,x) 

+ ^2f2^L(E,x) 


d/ FD 
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dEdx 


(3.12a) 
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(3.12b) 


where  /fD  is  the  Fermi-Dirac  function  and  NiiE^c)  is  the  energy  and  space  trap  density  in 
eV-1cnT3.  For  an  inelastic  tunneling  associated  with  energy  loss  due  to  multiphonon 
lattice  relaxation  [Kir89,  Lee02,  Tew94],  the  effective  slow  time  constant  tsiow(E,x ) in  a 
direct  tunneling  regime  can  be  expressed  under  the  Taylor’s  series  approximation  as 
follows 


TdoJE’X)  = Tslo*o(E)QXP 


kT 


+ A{E,x) 


= r0(£)-exp(/l,  •*) 


(3.13a) 


with 
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Tq(E)  = rslowo(EytXP\ 
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* y 


where  cr,v/""  is  the  capture  cross  section  of  slow  traps,  v,/,  is  the  thermal  velocity,  and  ns  is 
the  surface  electron  concentration.  Assuming  that  NiiEjc)  is  uniform  in  energy  and  in 
space  over  a distance  d into  the  oxide  and  tq(E)  is  a slowly  varying  function  of  energy, 
Equations  (3.12a)  and  (3.12b)  are  found  to  be  [Ghi87,  Kum77] 
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(3.14b) 


where  N,sIow(  Ef ) = Nj{Ejc)-d  is  the  oxide  slow  trap  density  per  unit  of  area  at  the  Fermi 
energy  level  in  eV'cm'2,  and  in  strong  inversion,  tq,,,  = l/(cr,i/ovV,/,ttv)  oc  \l{CoxVox) 


[Tew94], 
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3.2.2  Lorentzian-Modulated  Shot  Noise  Model 

Shot  noise  is  the  electrical  fluctuation  due  to  discreteness  of  the  charge  that 
provides  direct  information  on  the  correlation  of  different  current  pulses.  Full  shot  noise 
is  expected  if  the  electrons  tunneling  through  the  oxide  layer  show  Poisson  statistics.  The 
analysis  of  shot  noise  is  performed  by  introducing  the  Fano  factor  F defined  as  F = S/(f) 
/ ( 2ql  ),  Si  ( / ) being  the  spectral  density  of  current  fluctuations  as  a function  of 
frequency/and  I the  current  flowing  in  device.  Shot  noise  is  suppressed  ( 0 < F < 1 ) or 
enhanced  ( F > 1 ),  depending  on  correlation  introduced  in  the  behavior  of  carriers.  Shot 
noise  suppression  is  explained  by  Pauli  exclusion  and  Coulomb  interaction  acting  on  each 
trap  [Ian97,  IanOO,  Vli81],  Coulomb  repulsion  effectively  forbids  a second  electron  to 
occupy  the  other  state  and  inhibit  tunneling  of  electrons  in  a surrounding  region,  leading 
to  suppressed  shot  noise.  In  contrast,  as  an  electron  tunneling  into  the  trap  from  the 
substrate  channel  raises  the  potential  energy  of  trap  states,  more  states  are  available  for 
successive  tunneling  events  and  the  probability  to  tunnel  into  the  trap  increases,  implying 
that  electrons  entering  the  trap  are  positively  correlated  and  enhanced  shot  noise  is  to  be 
expected  [Ian98,  Rek99].  In  the  presence  of  generation-recombination  processes  with 
electron  traps,  it  was  found  that  the  potential  barrier  fluctuations  induced  by  the 
fluctuating  fixed  charge  by  long-range  Coulomb  interaction  are  at  the  origin  of  shot  noise 
enhancement  [Vas02]. 

We  describe  a model  for  noise  and  transport  based  on  trap-assisted  two-step 
tunneling  in  a MOS  structure,  as  shown  in  Figure  3.3.  Introducing  the  generation  rate, 
i.e.,  the  transition  rates  from  the  substrate  ( gsl ) and  from  the  gate  ( ggl ) to  the  unoccupied 
trap,  and  the  recombination  rate,  i.e.,  the  transition  rates  from  the  occupied  trap  to  the 
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substrate  ( r,s ) and  to  the  gate  ( r,g ),  the  generation  and  recombination  rate  can  be 
represented  by 


to . Foxi ) = qosrNsr  (x,  qFoxlx  - ^ )[l  - /,  (x,  ^F„,x  - d>„  )\g{<Db , x,  FnxX ) 

SgM,Sox  ~x’Fox2)  = qv!m'N!as' (x,qFoxXx- 0h )[l  - /, {x,qFoxXx  — <ph  )\g(</>,,tox  ~x,Fox2) 
r,M^ox -x,Foxu)  = qtrf* Njasl  (x, qFoxXx  - <pb )/,  (x, qFoxXx- <ph ) r(0, , tox -x,Fox2) 
rM,x,Fml)  = qcrfas'Nfas' (x, qFoxXx - <pb )ft (x, qFoxXx -<pb)  r(<pb , x, FoxX ) 

(3.15) 

where  g((pbrX,Foxi,2)  and  r((pbj:,Fox  1,2)  are  the  uniform  generation  and  recombination  rate, 
respectively,  and  Ntfas'{x,E)  is  the  oxide  fast  trap  density  at  a distance  x and  energy  E in 
eV'cm'2 . 


Figure  3.3  Model  for  Lorentzian-modulated  shot  noise  and  transport  based  on  inelastic 
trap-assisted  tunneling  in  a MOS  structure. 
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The  time  constant  Tfast  is  associated  with  the  rate  equation  for  carrier  number 
fluctuations  AN  with  Tfasl  1 and  Tfas,2  being  the  time  constants  related  to  the  variation  of  the 
generation-recombination  rate  due  to  the  carrier  number  fluctuations  through  the 
substrate  and  gate,  respectively.  It  is  assumed  that  N,fas'  and  f are  independent  of  the 
number  of  carriers  in  the  channel  per  unit  area,  N.  When  the  carrier  flux  is  one 
directional,  the  characteristic  time  constants  around  the  steady  state  values  are  given  by 


1 1 1 

— 1 

^ fast  (0b  F0x\,2  ) ^ fasti  to,  X'  Fox>  ) ^ fast  2 (01  > Kx  X*  Fox 2 ) 
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(3.16b) 


where  N is  the  average  number  of  N.  From  Equations  (2. 12),(3. 15),  and  (3.16),  we  can 
get  the  expression 
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where, 
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(3.18c) 

(3. 1 8d) 
(3.18e) 

(3.1 8f) 


These  equations  include  the  effect  of  charge  stored  in  the  trap  states  on  the  oxide  field 
during  tunneling  through  the  trap  states.  The  transition  rates  are  modulated  by  the 
statistical  fluctuation  of  trap  occupancy  in  the  oxide.  Therefore,  the  fluctuations  in  the 
recombination  rate  rlg  and  the  generation  rate  gst  around  the  steady  state  values  are  given 
by 

A rlg  ( N,t ) = — + R(N,t),  A g„  (N,t)  = — + G(N,t) 

T fast  2 T fasti 


(3.19) 
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where  AN  = N -N , and  R(N,t ) and  G(/V,r)are  the  driving  noise  forces,  respectively. 
It  is  shown  in  Equation  (3.19)  that  the  recombination  and  generation  rate  are  decomposed 
into  two  contributions  where  the  first  term  is  related  to  number  fluctuations  and  the 
second  term  is  due  to  the  full  shot  noise  associated  with  the  random  recombination  and 
generation  transition,  respectively.  Making  a Fourier  analysis  of  Arlg(N,t ) in  Equation 

(3.19)  and  assuming  the  correlator  for  Poissonian  current,  we  can  get  the  noise  spectrum 
for  the  gate  leakage  current 


Sslh,”(f)  = S?  + 
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with 


S°  =^^f=92^(/)^;(/)  = 5G(/)-5*(/)  (3.21a) 

S*(f)  = If  • /f‘  = q1  S,„  (/)  • Si „ (/)  = SAQ (/)  • (/)  (3.21b) 

S** (/)  = /f  -If  = q2  R e[SR(/)-5^(/)]  = Re[5e(/)-5;e(/)]  (3.21c) 


where  Ref  ] denotes  the  real  part  of  the  expression  [ ],  the  asterisk  refers  to 
complex  conjugate  quantities,  the  superscripts  of  the  noise  spectrum  S(f)  are  the 
fluctuating  quantities  considered,  and  IgQ  and  IgAQ  are  the  current  components  due  to  the 
total  inversion  charge,  Q = qN,  and  the  total  inversion  charge  fluctuation,  A Q = qAN, 
respectively.  Introducing  the  charge  fluctuations  A Q around  a steady  state  and  using  the 


Langevin  equation  for  A<2,  we  have 
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~AQ  + 
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(3.22) 


Making  a Fourier  transform  and  arranging  Equation  (3.22)  with  the  help  of  some 
algebra,  we  obtain 
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Assuming  that  Sq(/)  = Sg(f)  = 2ql  g and  substituting  Equations  (3.23a)  and 

(3.23b)  into  Equation  (3.20),  we  can  obtain  the  spectral  density  of  current  fluctuations  as 
follows 
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(3.24) 


This  formula  can  be  either  larger  or  smaller  than  the  full  shot  noise  2 ql  , 

depending  on  the  sign  of  characteristic  time  constants  Tfml(<pbJ,xl2,F0Xl2)  and 

T fas,2  (0,  » x2  ’ Foxi ) • Note  that  an  enhancement  in  the  transition  rate  occurs  when  the  time 

constant  becomes  a negative  value,  resulting  from  the  positive  correlation  between 
current  pulses.  Physically,  the  negative  time  constant  means  that  the  change  of  the  carrier 
flux  due  to  the  number  fluctuations  of  carriers  decreases  [Rek99,  RekOO]. 
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3.3.  Experimental  Verification 

Figure  3.4  shows  a typical  gate  current  noise  spectrum  measured  at  low  voltages 
for  n+  poly-Si  nMOSFET  with  ultrathin  nitrided  oxide.  The  noise  spectrum  clearly  shows 
a 1 If  frequency  dependence  with  the  power  factor  y = 1.0  ~ 2.0.  Gate  current  noise 
increased  with  the  increase  in  gate  voltage  in  the  wide  range  of  frequency,  while  slopes  of 
the  frequency  dependence  deviate  from  unity  with  increasing  voltage.  This  implies  that 
more  slow  trap  states  are  available  for  trap-assisted  tunneling  with  increasing  voltage. 
The  simulation  results  of  the  gate  current  noise  spectral  density  are  also  shown  as  a 
function  of  frequency  for  different  gate  voltages.  An  excellent  agreement  between  the 
measurement  and  simulation  at  both  low  and  high  frequencies  is  observed.  This  is 
attributed  to  the  fact  that  the  noise  model  is  based  on  the  trap-related  mechanism, 
resulting  in  low  frequency  dissipation  in  the  conductance  of  oxide  traps  and  producing 
the  Lorentzian-modulated  shot  noise  associated  with  G-R  processes  at  higher  frequencies. 

Figure  3.5  shows  the  comparison  of  experimental  data  and  simulation  results  of 
the  gate  leakage  current  noise  at  low  frequencies.  The  gate  current  noise  for  n+  poly-Si 
nMOSFET  was  measured  at  room  temperature  as  a function  of  dc  gate  current  with 
frequency  as  parameters,  as  shown  in  Figure  3.5  (a).  It  is  observed  that  there  are 
deviations  from  Ig  scaling  at  lower  current.  This  implies  that  the  low  frequency  noise 
behavior  at  low  voltages  is  dominated  by  a shot  noise  component.  At  higher  current 
level,  however,  the  simulation  results,  proportional  to  Ig2,  are  in  good  agreement  with  the 
measured  data  taken  at  low  frequencies  ( 10  Hz  and  50  Hz  ).  This  quadratic  gate  current 
dependence  is  expected  by  the  barrier  height  fluctuation  model  expressed  by  Equations 
(3.9)-(3.1 1).  In  order  to  investigate  the  behavior  of  1//  y noise  in  more  detail,  the 
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normalized  gate  current  noise  as  a function  of  gate  voltage  at  low  frequencies  is  shown  in 
Figure  3.5  (b).  The  normalized  noise  power  at  high  voltages  is  constant  as  a function  of 
gate  voltage,  indicating  that  both  the  low  frequency  noise  and  current  are  due  to  the  same 
ITAT  process  through  the  oxide  barrier. 

In  Figure  3.6  (a),  the  measured  Fano  factor  for  n+  poly-Si  nMOSFET  is  plotted 
as  a function  of  gate  voltage  for/=  10  kHz.  It  can  be  noticed  that  as  the  gate  voltage 
decreases,  the  measured  Fano  factor  of  nMOSFET  at  high  frequencies  reaches  a value 
close  to  10  . For  increasing  voltages,  the  Fano  factor  asymptotically  approaches  unity. 
The  Fano  factor  was  simulated  by  the  proposed  model  of  Equation  (3.24)  with/=  10 
kHz,  tsiON  = 3.21  nm,  EsiON  = 5.7eo,  <f>b,siON  = 2.6eV,  mox-  0.4m0,  N/as>  = 2xl012  eV'cm'2, 
x,  = 0.50 tsiON,  and  <p,=  2.33eV  at  T - 300K.  It  is  found  that  the  proposed  noise  model  has 
a good  agreement  with  the  experimental  results.  In  order  to  explain  the  enhanced  shot 
noise  with  decreasing  voltages,  inverse  time  constants  as  a function  of  gate  voltage  were 
calculated  by  using  the  expressions  of  Equations  (3.17)  and  (3.18),  as  shown  in  Figure 
3.6  (b).  An  enhancement  in  shot  noise  occurs  as  inverse  time  constants,  l/r/ut(1  and 

l/r/U5, 2,  approach  the  same  value  with  decreasing  voltages.  By  contrast,  full  shot  noise 

behavior  is  recovered  as  the  difference  between  inverse  time  constants  increases.  This  is 
due  to  the  positive  feedback  between  trap  states  and  tunneling.  When  an  electron 
tunneling  from  the  inversion  layer  to  the  traps  is  captured,  the  potential  barrier  increases, 
allowing  a decrease  of  carrier  flux  ( a negative  time  constant,  Tfas,\ ).  On  the  other  hand, 
when  a captured  electron  is  emitted  from  the  traps  to  the  gate,  the  potential  barrier 
decreases,  leading  to  an  increase  of  carrier  flux  ( a positive  time  constant,  Tfasa  ) and  thus 
an  enhanced  shot  noise  [Vas02],  Figure  3.7  illustrates  the  influence  of  trap-related 
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Frequency,  f(Hz) 


Figure  3.4  Comparison  of  measured  data  and  simulation  results  of  the  gate  current  noise 
spectral  density  as  a function  of  frequency  for  n+  poly-Si  nMOSFET  with 
tsiON  = 3.21nm,  Nsub  = 6.0xl017cm'3,  Npoiy  = 1.85xl020  cm'3,  eSiON  = 5.7e0, 
(pb,siON=  2.6eV,  mox  = 0.4m0,  N,rap  = 2xl012  cm'2,  and  N,slow  = Nfast  = 2xl012 
eV'1cm'2atr=300K. 
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Figure  3.5  Comparison  of  measured  data  and  simulation  results  of  (a)  gate  current  noise 
as  a function  of  gate  current  and  (b)  normalized  gate  current  noise  as  a 
function  of  gate  voltages  at  low  frequencies  for  n+  poly-Si  nMOSFET  with 
tsiON  = 3.21nm,  EsiON  = 5.7£o,  (pb.siON  = 2.6eV,  mox=  0.4 m0,  N,slow= 2xl012  eV'1 
cm-2,  x,  = 0.50tsiON , and  0,=  2.33eV  at  F=300K 
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Figure  3.6.  Gate  bias  dependence  of  (a)  Fano  factor  at  /=  10  kHz  and  (b)  inverse  time 
constants  simulated  by  the  proposed  model  with  tSio\  = 3.2  lnm,  Nsub  = 
6.0x10  cm  ',  Npo/v=  1.85xl0"°cm 3,  EsiON  = 5.7eo,  (f>b,siON  = 2.6eV,  mox  = 
0.4m0,  and  N,fasl  = 2x  1012  eV  'cnf2,  x,  = 0.50 tSi0N  , and  0,  = 2.33eV  at  T = 
300K. 
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Figure  3.7  Illustration  of  the  influence  of  trap  energy  levels  and  trap  positions  on  Fano 
factor  at  /=  10  kHz  and  = 0.5  V. 
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parameters  on  the  Fano  factor  F at  /=  10  kHz  and  Vg  = 0.5  V.  The  Fano  factor  peaks  for 
trap  located  in  the  middle  of  the  oxide  layer,  and  decreases  for  traps  at  deeper  energy 
levels.  This  is  consistent  with  the  trap-dependent  current  characteristics  [Lee02], 
implying  that  the  shot  noise  enhancement  is  due  to  the  trap-assisted  tunneling  associated 
with  the  generation-recombination  process. 

3.4  Conclusion 

A physical  and  analytical  model  describing  the  gate  leakage  current  noise  in 
ultrathin  gate  oxide  MOSFETs  was  presented.  In  order  to  characterize  the  noise  behavior 
exhibited  by  the  thin  nitrided  oxide  devices,  a physics-based  noise  model  of  gate  leakage 
current  was  proposed,  which  is  based  on  the  trap-related  tunneling  processes.  The  low 
and  high  frequency  noise  models  are  based  on  the  barrier  height  fluctuation  model  and 
Lorentzian-modulated  shot  noise  model  of  the  gate  leakage  current  to  explain  the  excess 
noise  behavior,  respectively.  Slow  traps  in  oxide  result  in  low  frequency  dissipation  in 
the  conductance  of  oxides,  and  fast  traps  produce  the  Lorentzian  modulated  shot  noise 
associated  with  generation-recombination  process  at  higher  frequencies.  It  was  verified 
from  the  measured  data  that  the  newly  developed  noise  model  can  accurately  predict  the 
excess  noise  behavior  of  MOS  devices  with  ultrathin  gate  oxide.  In  particular,  a new 
formula  for  the  gate  leakage  current  noise  is  suitable  to  be  implemented  in  circuit 


simulators. 


CHAPTER  4 

DEFECT  SPECTROSCOPY  USING  1// NOISE  OF  GATE  LEAKAGE  CURRENT 

4.1  Introduction 

With  the  introduction  of  ultrathin  oxides  in  MOS  technologies,  the  gate  leakage 
current  and  its  associated  noise  becomes  significant  even  at  low  frequencies  [Lee02, 
Mom98,  Ner95,  Sch94b].  This  has  not  only  major  implications  for  advanced  circuit 
design  but  also  impacts  the  experimental  study  of  defects  in  these  advanced  structures 
[Mom98,  Nic67b].  It  has  long  been  known  that  the  low  frequency  noise  behavior  of  MOS 
structure  is  directly  related  to  trap  properties  [Cel85,  Chr68,  Hun90a,  Nic67b].  Due  to  the 
effects  of  Si-Si02  interface  trap  states  on  the  reliability  and  quality  of  MOS  structures,  its 
electrical  and  physical  characterization  has  been  a focus  of  intensive  study  [Nic67a, 
Nic82],  However,  for  ultrathin  (<3nm)  gate  oxides,  an  accurate  and  reliable 
characterization  method  of  the  trap  properties  has  not  yet  emerged  [KonOl]. 

The  most  frequently  used  techniques  to  electrically  characterize  Si-Si02  interface 
are  the  conductance  method  [DuvOla,  KonOl,  Nic67a,  Nic82],  the  C-V  method  [CleOO, 
GreOl,  OmuOO],  and  the  charge  pumping  method  [MasOl].  However,  the  conductance 
method  loses  its  resolution  due  to  the  large  tunneling  current  in  ultrathin  gate  oxide  MOS 
devices  and  requires  elaborate  analysis  procedures.  In  addition,  the  C-V  method  needs 
quantum-mechanical  and  poly-silicon  depletion  correction  and  charge  pumping  method 
requires  an  additional  correction  for  excessive  leakage  current  that  can  lead  to  large 
errors  in  the  extracted  trap  parameters. 
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In  this  Chapter,  we  will  present  our  attempt  to  probe  the  trap  properties  in  MOS 
device  with  ultrathin  gate  oxides  by  measuring  the  1 If  low  frequency  noise  of  the  gate 
leakage  current.  Since  the  1 If  noise  depends  strongly  on  the  charge  fluctuation  of  the  trap 
states  in  MOSFETs,  information  about  the  trap  properties  can  be  extracted  from  1 If  low 
frequency  noise  measurements.  There  have  been  several  papers  on  the  applications  of  the 
1 If  noise  technique  to  determine  the  trap  densities  [Cel88,  Jay89,  Man99,  Vil99]. 
However,  all  the  low  frequency  noise  measurements  were  performed  on  the  channel 
current  in  MOSFETs  with  relatively  thick  oxides.  Besides  the  channel  current,  the  gate 
leakage  current  should  be  considered  as  a low  frequency  noise  source  [Lee02,  Mon98]. 
As  we  will  demonstrate,  the  gate  leakage  current  noise  occurring  at  microscopic  level  can 
be  a powerful  and  sensitive  tool  in  evaluating  the  trap  properties  in  ultrathin  gate  oxide 
MOSFETs. 

For  an  accurate  and  sensitive  extraction  of  trap  properties,  a new  technique  is 
proposed.  The  basic  theory  of  1 If  noise  in  terms  of  complex  admittance  of  the  trap  states 
will  be  presented  in  Section  4.2.  In  order  to  investigate  the  effects  of  the  spatial  trap 
distribution  on  the  complex  admittance  and  1 If  low  frequency  noise,  analytical 
expressions  are  derived  for  uniform  and  exponentially  decreasing  trap  distributions  in 
space.  After  presenting  the  experimental  details  in  the  next  section,  we  will  show  how  the 
characterization  and  extraction  procedure  can  be  applied  to  MOS  devices  featuring 
ultrathin  oxide  using  the  1 If  low  frequency  noise  of  the  gate  leakage  current.  Trap 
densities  and  capture  cross  sections  are  accurately  determined  as  a function  of  trap  energy 
level  and  trap  distance,  and  the  noise  characteristics  are  revealed  in  terms  of  the  complex 
admittances  of  the  trap  states.  Finally,  Section  4.5  presents  some  conclusions. 
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4.2  Ilf7  Noise  and  Admittance  Model 


An  analytical  noise  model  for  the  gate  leakage  current  through  a SiC>2  barrier  is 
presented  in  the  direct  tunneling  regime.  It  has  been  generally  accepted  that  trap-related 
processes  are  the  major  cause  of  1 If  fluctuations  in  MOSFETs.  The  1 If  noise  observed 
at  low  frequencies  might  be  caused  by  the  fluctuation  of  the  barrier  height  or  shape  due  to 
the  fact  that  the  frequency-dependent  conductance  of  the  oxide  traps  produces  thermal 
noise  resulting  in  a surface  potential  fluctuation,  as  shown  in  Figure  3.2.  This  barrier 
height  fluctuation  in  turn  modulates  the  tunneling  transmission,  and  causes  the 
fluctuations  in  the  tunneling  current. 

For  the  extraction  of  the  trap  properties,  the  noise  technique  may  be  used  by 
combining  the  complex  MOS  admittance  of  the  oxide  trap  states  [Mal71,  Saw79] 


Rearranging  Equations  (3.9)  and  (3.11)  with  the  help  of  some  algebra,  we  obtain  a 
calculable  expression  in  terms  of  the  frequency-dependent  parallel  conductance  and 
capacitance 


Yp(f)  = Gp(f)  + j2nfCp(f ) 


(4.1) 


(4.2) 


Provided  the  noise  data  over  the  necessary  frequency  range  in  Equation  (4.2)  are 
available,  the  trap  information  can  be  obtained  from  a Cole-Cole  plot  of  the  complex 
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admittance  in  which  the  frequency  loci  of  Yp{f)l27tf  are  drawn  on  the  complex 
Cp  ~Gp(f)/27tf  plane.  In  this  plot,  the  diameter  of  the  arc  is  equal  to  the  trap  density 

qN,,  and  other  trap  parameters  are  estimated  by  the  analysis  of  the  detailed  geometry  of 
the  arc  [Mal71,  Saw79].  Similarly,  using  Equations  (3.9)  and  (3.11),  and  defining 
Gp(f)/27tf  = N;low(Ef)-G(f ) and  C,(/)  = N?ow(Ef)-C(f),  we  get  an  analytical 
expression  for  the  interface  trap  density 


where 


^(Vox,f)  = 


V ^ y 
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(4.4c) 


For  an  exponentially  decreasing  trap  distribution  in  space,  we  consider  the  form 


[Bor99] 


Nt(E,x ) = Nt  - exp 


f ^ 


V oc) 


(4.5) 


where  a is  the  characteristic  decay  depth  of  the  trap  density  into  the  oxide.  Substituting 
Equation  (4.5)  into  Equations  (3.12a)  and  (3.12b)  and  assuming  that  in  most  cases  a 
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>10A,  the  integral  can  be  approximated  analytically.  The  resulting  parallel  conductance 
and  capacitance  are  found  to  be 


Cp(f)  = 


q2Nj‘°w(Ef ) [ tan-1  [2flfr0>„ ] 
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(4.6b) 


4.3.  Experimental  Details 

In  order  to  extract  the  trap  properties  of  the  MOS  structure,  knowledge  of  the 
current-voltage  and  low  frequency  noise  characteristics  is  required.  Prior  to  the  noise 
measurements,  all  devices  were  examined  for  gate  leakage  current  by  measurements  of 
their  dc  current-voltage  characteristics.  The  low  frequency  noise  measurements  were 
performed  for  frequencies  ranging  between  1 Hz  and  100  kHz  at  room  temperature.  The 
gate  voltage  was  swept  from  0.4  V to  1.32  V.  Figure  4.1  shows  a typical  gate  current 
noise  spectrum  measured  at  low  voltages  for  a n+  poly-Si  nMOSFET  with  ultrathin 
nitrided  oxide.  The  noise  spectrum  clearly  shows  a 1 If  frequency  dependence  with  the 
power  factor  y - 0.9  ~ 1.7. 

For  the  purpose  of  validating  the  proposed  method,  the  conductance 
measurements  are  performed  using  a standard  LCR  meter.  The  small  signal  conductance 
was  measured  at  different  gate  biases  in  the  inversion  regime,  sweeping  the  signal 
frequencies  from  100  Hz  to  1 MHz.  For  a MOSFET  with  the  source  and  drain  tied 
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Figure  4.1  Gate  current  noise  spectral  density  for  several  current  values  with  tox,eq  = 
2.2nm,  Nsub  = 6.0xl0!7cm'3,  Npoiy  = 1.85xl020  cm'3,  and  Ag  = 140x70x4  //m2 
at  T = 300K 


together,  measurements  were  made  in  the  parallel  mode,  as  illustrated  in  Figure  3.2 
[Eat72,  Lee02].  From  the  noise  equivalent  circuits  including  the  tunneling  dc 
conductance  GUdc  and  the  inversion  resistance  /?„  the  equivalent  parallel  conductance 
Gp(  f )/2nf  is  related  to  the  measured  conductance  Gm(  f ) and  capacitance  Cm(  f ) as 
follows  [Eat72,  VogOO] 


Gp(f) 
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where  Gc(  f ) and  Cf(  / ) are  the  conductance  and  capacitance  corrected  for  inversion 
resistance  /?„  respectively.  In  Equation  (4.7a),  the  dc  tunneling  conductance  G,jc  ( Vg ) is 
subtracted  from  Gc(  f ) to  obtain  the  ac  conductance.  G,jc{  Vg  ) is  obtained  by 
differentiating  the  static  gate  leakage  current  curve  taken  as  a function  of  gate  voltages. 

4.4  Extraction  Procedure  and  Results 

The  shape  of  the  low  frequency  noise  spectrum  is  strongly  dependent  on  the  trap 
distribution  in  the  oxide.  Actually,  although  traps  are  distributed  over  space  jc  into  the 
oxide  as  well  as  in  energy  E over  the  band  gap,  a complete  characterization  of  traps  in  the 
( x,  E ) window  is  extremely  difficult.  It  is  widely  accepted  that  trap  states  are  classified 
by  two  distinct  groups.  The  first  is  ‘fast  states’  with  capture  cross  sections  of  about  10'16 
cm  , and  the  second  is  ‘slow  states’  characterized  by  capture  cross  sections  less  than 
about  10"16  cm2  [Kir89].  The  slow  states  play  a fundamental  role  in  l//low  frequency 
noise  and  are  due  to  defect  states  extended  into  the  oxide  layer.  From  the  tunneling  point 
of  view,  the  effective  capture  cross  section  of  a slow  trap  located  at  a distance  * from  the 
interface  is  reduced  by  [Kir89,  Lee03] 


aeff  (x,  E)  = a;,ow(E)  ■ exp[-  X,  (Vox ) • x] 


(4.8a) 
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In  the  case  of  the  low  frequency  noise  technique,  the  tunneling  depth  * corresponds  to  a 
particular  measurement  signal  frequency  and  is  expressed  as  [Lee03] 


x — — — • In 
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2 7f-T0(E) 


(4.9) 


In  order  to  probe  the  trap  energy  level,  the  band  gap  is  scanned  by  changing  the  surface 
potential  i//s,  which  is  done  by  varying  the  gate  voltage.  Therefore,  we  can  access  the 
( x,  E ) window  of  trap  states  in  the  oxide,  by  using  frequency  and  bias-dependent  noise 
measurements.  Figure  4.2  (a)  shows  N,slow(  f ) curves  as  a function  of  frequency,  directly 
extracted  from  the  measured  noise  power  spectral  density  using  Equations  (4.3)  and  (4.4). 
This  trap  spectrum  clearly  shows  the  dynamic  behavior  of  oxide  trap  states  in  terms  of 
frequencies.  N,shw(f ) curves  increasing  with  frequencies  reveal  a l/f  (y  < 1.0).  The  trap 
densities  are  plotted  as  a function  of  trap  energy  E and  trap  distance  x,  as  illustrated  in 
Figure  4.2  (b).  The  plot  clearly  visualizes  that  the  slow  trap  densities  estimated  near  the 
band  edge  and  at  the  Si-SiC>2  interface  are  higher  than  those  near  the  midgap  and  deeper 
distance. 

Moreover,  the  gate  leakage  current  noise  technique  enables  the  determination  of 
the  equivalent  parallel  conductance  of  oxide  traps  at  low  frequencies.  This  can  be  done 
by  using  Equation  (4.2)  and  measured  noise  data.  Figure  4.3  (a)  shows  the  conductance 
loss  Gp{f)!2nf  extracted  from  the  noise  measurements  for  n+  poly-Si  nMOSFET. 
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Figure  4.2  Trap  density  directly  extracted  from  the  measured  low  frequency  noise 
spectra  (a)  as  a function  of  frequency  and  (b)  trap  energy  and  trap  distance. 
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It  was  assumed  that  the  trap  is  uniformly  distributed  in  space  for  Vg  = 0.8  V and  the  trap 
is  exponentially  decreasing  in  space  for  Vg  = 1.2  V.  The  low  frequency  plateau  in  the 
conductance  loss  is  responsible  for  the  1 If  noise  and  the  curvature  determines  the 
capture  cross  sections  of  the  oxide  traps  [DuvOla,  Kir89].  The  corresponding  capture 
cross  sections  are  obtained  from  the  extracted  conductance  loss  data  and  Equation  (4.8), 
and  are  plotted  in  Figure  4.3  (b).  It  is  interesting  to  note  that  the  capture  cross  sections 
depend  strongly  on  the  surface  potential  and  is  nearly  an  exponential  function  of  the 
surface  potential  with  a slope  of  q / kT.  The  strong  energy  dependence  of  the  capture 
cross  section  of  ‘slow  trap  states’  observed  at  low  frequencies  is  in  good  agreement  with 
the  previously  published  data  on  conductance  spectroscopy  techniques  [DuvOla,  Eat72, 
VogOO]. 

The  parallel  conductance  Gp(  f )/2nf  as  a function  of  frequency  is  measured  for 
different  gate  voltages.  Figure  4.4  (a)  shows  the  parallel  conductance  obtained  by  using 
Equation  (4.7)  and  Figure  3.2.  The  low  frequency  plateau  extracted  from  the  noise  results 
is  confirmed  by  the  conductance  measurements.  The  peak  of  Gp{  f )/2nf  reduces  with 
decreasing  gate  voltage  and  shifts  to  lower  frequency.  This  is  mainly  due  to  the  fact  that 
the  parallel  conductance  for  lower  gate  voltage  is  produced  by  smaller  trap  densities  and 
slower  time  constants  of  ‘fast  trap  states’.  As  shown  in  Figure  4.4  (b),  the  parallel 
conductance  curves  for  the  uniform  trap  distribution  with  the  tunneling  depth  0.5 t„x  and 
0.08  tox  into  the  oxide  are  plotted  for  comparison,  respectively.  The  existence  of  two 
distinct  traps  is  verified  by  the  model  fitting  with  the  measured  data.  The  fast  trap  states 
are  at  the  origin  of  enhanced  or  suppressed  shot  noise  associated  with  the  generation- 
recombination  process  [Lee02]. 
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For  the  purpose  of  investigating  the  dependency  of  trap  parameters  on  the 
complex  admittance  of  oxide  traps,  the  Cole-Cole  plots  for  different  trap  profiles  are 
calculated  using  the  theoretical  model.  Figure  4.5  (a)  and  Figure  4.5  (b)  illustrate  the 
shape  of  the  Cole-Cole  plots  of  the  trap  states  admittance  for  a given  frequency  range, 
assuming  that  the  trap  density  is  uniform  and  exponential  decreasing  in  space, 
respectively.  It  can  be  seen  that  the  arc  shape  depends  strongly  on  the  trap  distribution, 
and  that  the  computed  loci  for  the  higher  gate  voltage  have  a larger  arc  diameter  and 
lower  peak  frequency,  giving  a direct  determination  of  the  trap  parameters.  From  these 
results,  it  was  shown  that  a Cole-Cole  plot  obtained  by  the  measured  noise  data  could 
yield  qualitative  information  about  the  trap  properties. 

The  \/f  noise  of  MOS  structures  is  given  in  terms  of  the  complex  admittance  of 
trap  states,  because  the  noise  behavior  is  directly  related  to  the  trap  state  properties. 
Therefore,  in  order  to  investigate  the  dependence  of  trap  properties  on  the  low  frequency 
noise  in  more  detail,  the  Cole-Cole  plot  on  the  2-dimensional  Cp(f)-Gp(f)/ plane 

can  be  extended  into  3-dimensions.  Figure  4.6  (a)  shows  the  frequency  loci  of  the  noise 
in  terms  of  parallel  conductance  and  capacitance  in  3 dimensional  space,  assuming  that 
the  trap  density  is  uniform  in  space.  The  low  frequency  noise  reduces  with  increasing 
frequency,  giving  a quantitative  determination  of  trap  parameters.  Curves  of  Figure  4.6 
(b)  show  the  noise  power  spectrum  as  a function  of  conductance  for  the  parameter  of 
frequency.  For  Vg  = 0.6V,  the  noise  spectrum  becomes  independent  of  frequency  at  low 
frequencies,  occurring  at  frequencies  much  lower  than  the  conductance  peak  frequency. 
By  contrast,  for  Vg  = 1.32V,  the  noise  spectrum  has  its  highest  value  around  the 
conductance  peak  frequency,  tending  towards  a 1 If  (y  >1.0  ) type  frequency  dependence. 


Capture  Cross  Section,  a „ ( cm2 ) Conductance  Loss,  GpA/v  ( F/cm2 ) 
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Figure  4.3  Trap-related  parameters  extracted  from  the  low  frequency  noise 
measurements  for  n+  poly-Si  nMOSFET.  (a)  Conductance  loss  Gp{  f )/w, 
assuming  that  the  trap  is  uniformly  distributed  in  space  for  Vg  = 0.8  V and 
the  trap  is  exponentially  decreasing  in  space  for  Vg  = 1.2  V and  (b) 
capture  cross  section  of  oxide  traps.  ( tox,eq=  2.2  nm,  Nsub  = 6xl017cm'3, 
mox  = 0.4  mo,  and  T - 300K  ) 
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Figure  4.4  Conductance  loss  characteristics,  (a)  Measured  conductance  loss  G,,(f)/w  as 
a function  of  frequency  for  various  gate  voltages  and  (b)  theoretical  behavior 
of  Gp{  f )/w  for  the  slow  trap  ( d = 0.5  tox ) and  the  fast  trap  ( d = 0.08  tox  ), 
assuming  that  the  trap  is  uniformly  distributed  in  space. 
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Figure  4.5  Cole-Cole  plots  of  the  trap  states  admittance  calculated  from  the  theoretical 
model  ( frequency  step  A/  = 1 Hz  ),  assuming  that  (a)  the  trap  density  is 
uniform  in  space  and  (b)  the  trap  is  exponentially  decreasing  in  space 
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(a) 


Figure  4.6  Gate  current  noise  characteristics  ( frequency  step  A/  = 1 Hz  ),  (a)  as  a 
function  of  the  trap  states  admittance  for  a parameter  of  frequency,  and  (b) 
as  a function  of  the  conductance  loss  G,Jw,  assuming  that  the  trap  density  is 
uniform  in  space. 
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As  shown  in  Figure  4.6  (a)  and  (b),  it  is  apparent  that  these  plots  provide  a better 
physical  understanding  of  the  low  frequency  noise  source  and  give  information  about  the 
trap  states  responsible  for  1 If  noise  in  MOS  devices  with  ultrathin  gate  oxide. 

4.5  Conclusion 

The  trap  state  properties  of  MOS  devices  with  ultrathin  oxides  were 
characterized  by  1//  low  frequency  noise  of  gate  leakage  current.  The  trap  parameters 
such  as  slow  trap  densities  and  capture  cross  section  are  accurately  determined,  and  the 
low  frequency  loci  of  1 If  low  frequency  noise  are  clearly  shown  by  the  complex 
admittance,  revealing  the  dynamic  evolution  of  trap  states.  It  was  found  that  the  capture 
cross  sections  depend  strongly  on  the  surface  potential  and  decrease  exponentially  with 
the  surface  potential  towards  the  midgap.  The  1 If  noise  technique  allows  the 
determination  of  oxide  trap  densities  spectroscopy  over  the  accessible  frequency  range, 
and  the  values  near  the  band  edge  and  the  interface  are  higher  than  those  near  the  midgap 
and  deeper  distance.  It  was  also  demonstrated  that  the  parallel  conductance  spectroscopy 
can  be  reproduced  by  1//  noise  measurements  and  the  Cole-Cole  plots  can  be  used  to 
identify  the  different  trap  profiling  in  space.  Furthermore,  3-dimensional  plots  give  better 
insight  into  the  low  frequency  noise  behavior.  The  results  confirm  that  1//  noise 
technique  is  directly  applicable  to  extract  the  trap  properties  in  ultrathin  oxide  MOS 


devices. 


CHAPTER  5 

COMPREHENSIVE  NOISE  PERFORMANCE  AT  LOW  FREQUENCIES 

5.1  Introduction 

The  gate  tunneling  leakage  current  resulting  from  downsizing  CMOS 
technologies  is  increasingly  important  in  terms  of  device  operation  and  reliability 
[LeeOl,  Mom96,  Mom98,  Sch94].  As  the  gate  oxide  thickness  is  scaled  down  to  the 
direct  tunneling  regime,  the  tunneling  leakage  current  increases  drastically  and 
deteriorates  the  device  performance.  Moreover,  the  current  source  causes  a low  frequency 
noise  component  which  not  only  determines  the  lower  end  of  the  dynamic  range  of 
device  operation  in  the  base  band,  but  also  has  potentially  a major  impact  on  the  phase 
noise  of  nonlinear  circuits  in  the  GHz  region  through  up-conversion  [KluOO].  Since  the 
impact  of  the  gate  leakage  current  can  play  a significant  role  in  the  overall  device 
characteristics,  an  exhaustive  evaluation  of  device  performance  is  required  over  the  entire 
frequency  range. 

From  the  noise  modeling  point  of  view,  the  gate  leakage  current  introduces  a shot 
noise  current  source  at  high  frequencies  [Dan95],  in  addition  to  the  existing  channel 
thermal  noise  source  at  the  output  and  induced  thermal  noise  source  at  the  input.  The  high 
frequency  noise  sources  are  partially  correlated,  because  the  MOS  transistor  can  be 
treated  as  an  RC  distributed  network  meaning  that  the  high  frequency  gate  admittance 
contains  a conductive  component  [Cap88,  Zie86].  At  low  frequencies,  since  the  gate 
conductance  in  ultrathin  gate  oxide  MOSFETs  can  be  substantial  due  to  the  oxide  traps,  it 
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is  also  observed  that  a partial  correlation  between  the  drain  and  the  gate  current  noise 
sources  exists.  Therefore,  in  order  to  accurately  predict  the  noise  performance  at  low 
frequencies,  a physical  and  rigorous  modeling  of  the  correlation  effects  is  required. 

The  correlation  noise  measurements  at  low  frequencies  have  been  previously 
presented  for  heterojunction  bipolar  transistors  (HBT’s)  [Bru99,  BruOO,  Del97,  Jar97, 
MouOl],  and  the  influence  of  the  gate  leakage  current  on  the  high  frequency  noise 
performance  was  investigated  for  MESFET’s  and  MODFETs  [Puc74,  Van92].  To  date, 
no  detailed  study  of  low  frequency  correlation  and  noise  performance  with  an  effect  of 
gate  leakage  current  has  been  given  for  ultrathin  oxide  MOSFETs. 

In  this  Chapter,  a comprehensive  noise  performance  at  low  frequencies  is 
presented,  including  a physical  model  of  correlation  coefficient  and  its  measurements,  as 
well  as  the  minimum  noise  figure  and  the  optimum  source  impedance  in  the  presence  of  a 
gate  leakage  current.  The  correlation  coefficient  is  derived  on  the  basis  of  the  partition 
noise  theory  [Rig94,  Zie84,  Zie87]  and  the  BSIM4  gate  leakage  current  model  with  the 
source-drain  partition  [CaoOO],  and  is  compared  with  correlation  noise  measurements  in 
terms  of  the  current  ratio  of  the  gate  to  the  drain.  The  analytical  expression  of  low 
frequency  noise  performance  is  newly  formulated  using  the  frequency-dependent  noise 
parameters  P^,  RLF,  and  CLF,  and  the  simulation  results  are  graphically  presented  to 
investigate  how  the  overall  noise  performance  depends  on  the  gate  leakage  current  and 
the  correlated  noise  sources  at  low  frequencies.  It  is  confirmed  that  accurate  modeling 
and  understanding  of  physical  noise  sources  at  low  frequencies  are  indispensable  for 
optimum  designs  of  ultrathin  oxide  MOS  devices  and  circuits. 


83 


5.2  Noise  Theory 

5.2.1  Noise  Sources 

The  high  frequency  noise  behavior  of  MOSFETs  basically  arises  from  the 
channel  thermal  noise  source  at  the  output,  S]d  (/) , and  the  induced  thermal  noise  source 

at  the  input,  S ] (/) , as  shown  in  Figure  5.1.  These  two  sources  are  correlated  due  to  the 

capacitive  coupling,  and  the  corresponding  correlation  coefficient  CT(  f ) is  complex 
[Cap98,  Dan95,  Zie86].  In  the  low  frequency  range,  the  two  dominant  noise  sources  at 
the  output  are  the  1 If  noise  and  the  generation-recombination  (G-R)  noise  modeled  as 
follows 
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where  k is  the  Boltzmann  constant,  T the  absolute  temperature,  gmg  =dId/dV^  the 
transconductance,  (p'dlf  I fr")  the  1 If  equivalent  input  noise  resistance,  and  pfR  and  r. 

characterize  a G-R  noise  resistance  and  a time  constant  for  the  ith  component, 
respectively.  The  physical  origins  of  1 If  noise  can  be  explained  by  the  correlated  number 
and  mobility  fluctuation  model.  The  correlated  model  is  based  on  the  fact  that 
fluctuations  in  the  occupancy  of  the  oxide  traps  induce  correlated  fluctuations  in  the 
carrier  number  and  surface  mobility  [Hun90b].  G-R  noise  in  MOSFETs  is  caused  by 
random  trapping-detrapping  of  carriers  at  the  defect  centers  in  the  depletion  region.  The 
G-R  noise  is  usually  overshadowed  by  the  higher  1 If  surface  noise  component. 

In  the  case  of  MOSFETs  with  ultrathin  gate  oxide,  however,  the  low  frequency 
noise  performance  should  be  modified  due  to  the  presence  of  the  gate  leakage  current  Ig, 
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Figure  5.1  Admittance  representation  of  the  physical  noise  sources  of  ultrathin  oxide 
MOSFET  with  a gate  leakage  current  noise  source 


as  schematically  represented  in  Figure  5.1.  The  low  frequency  noise  source  at  the  input, 
SfF(f),  physically  follows  the  same  mechanism  as  observed  in  the  output,  and  is 

represented  as 
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where  (p'Jf  / fr* ) is  the  1 If  gate  noise  resistance  and  g g =dl  g /dVgs  is  the  gate 

conductance.  Note  that  the  low  frequency  gate  noise  is  proportional  to  \/f,  while  the 
high  frequency  gate  noise  shows  the  quadratic  dependence  on  frequency  [Cap88,  Zie86]. 
The  low  frequency  noise  sources  between  the  input  and  the  output  are  partly  correlated, 
and  the  cross  correlation  coefficient  &F  may  not  be  negligible.  The  Lorentzian- 
modulated  shot  noise  source  Sf  (/)  is  generated  by  a trap-assisted  tunneling  (TAT) 
through  the  gate  oxide,  and  can  be  analyzed  by  the  Fano  factor  Ff“"°  defined  as 
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5.2.2  Noise  Models  and  Correlation  Function 

In  a direct  tunneling  gate  oxide  MOSFET,  there  is  a gate  leakage  current 
component  Ig  besides  the  drain  current  component  Ij.  The  ratio  of  gate  leakage  current  to 
drain  current,  Ig/Ij,  decreases  approximately  in  proportion  to  the  square  of  channel 
length,  Lr  and  typically  lies  between  1 and  10‘6,  depending  on  channel  length  [Mom96]. 
As  shown  in  Figure  5.2,  the  source  current  Is  is  partitioned  into  Id  flowing  to  the  drain 
with  a partition  factor  A,  and  Ig  flowing  to  the  gate  with  a partition  factor  (1-  A).  The 
average  value  of  the  partition  factor  A represents  the  probability  that  electrons  arrive  at 
the  drain.  Assuming  that  electrons  arriving  at  the  drain  or  at  the  gate  depend  only  on  the 
starting  point  of  the  electrons  in  the  channel,  the  drain  current  noise  SFF  (/) , the  gate 

leakage  current  noise  SFF  (/)  and  the  cross  noise  spectrum  SFF,  (/)  at  low  frequencies 
are  given  by  [Rig94,  Zie84,  Zie87] 
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where  Sff  (/)  is  the  source  current  noise  before  partitioning  into  the  drain  and  the  gate. 

The  cross  correlation  coefficient  at  low  frequencies  between  the  drain  and  the  gate 
current  noise  is  defined  as 
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Figure  5.2  Noise  current  components  in  ultrathin  oxide  MOSFET.  Is  is  partitioned  into 
Id  at  the  partition  factor  A,  and  Ig  at  the  partition  factor  (1-  A). 
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Using  Equations  (5.4a)-(5.4c),  (5.5),  and  (5.6),  and  rearranging  in  terms  of  the  ratio  of  the 
gate  leakage  current  noise  to  the  drain  current  noise,  S" (/)/ S" (/) , we  obtain  an 
expression  for  the  correlation  function  as  follows 


(5.7) 


For  the  evaluation  of  the  correlation  function  at  low  frequencies,  it  is  necessary  to 
obtain  a physical  model  of  S" (/)/ S" (/) , taking  the  noise  generation  process  into 

account.  Under  the  assumption  that  at  low  frequencies  both  the  drain  and  gate  noise  are 
caused  by  a random  trapping-detrapping  of  carriers  in  localized  trap  states  near  the 
Si-SiC>2  interface,  we  can  use  the  Klaassen’s  formula  for  the  l//r  noise  current  generator 
at  y and  y along  the  channel  [Kla7 1 ] 

AS"  (y,y,f)  = Hd°(y)-—— l- -S(y-y)  for  / = 0,1  (5.8) 

fr'"  8(v,(y)) 

where  5(y)  is  a delta  function,  H{y)  the  local  noise  source  related  to  a trap  efficiency, 
g(Vj)  the  channel  conductance  per  unit  length  at  y,  Vo(y ) the  channel  voltage  at  y,  and 
Vi(y)  the  perturbation  voltage  induced  by  the  gate  leakage  current  ( Vo(y)  » UiOO  ).  If 
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the  noise  source  is  distributed  uniformly  over  the  range  0 < y < yi,  the  noise  spectrum 
ratio  may  be  represented  by 


S"(f) 

S?if) 


^^w'(y)l^MS{y~y)dydy 


(5.9) 


where  L is  the  effective  channel  length.  Assuming  Hd  (y)  = H*  (y)  and  carrying  out  the 
integration,  bearing  in  mind  that  the  drain  and  gate  current  are  respectively  [Rao71], 


= g(v*W) 


dVjAy) 

dy 


and 


i'=8(vt(y)) 


dV\(y) 

dy 


yields  a simple  expression  as  follows 
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In  this  formula,  the  channel  voltage  Vo(y)  as  a function  of  y is  given  by  [Tau98] 
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where  m is  the  body-effect  coefficient  ranging  from  1.1  to  1.4,  VYthe  threshhold  voltage, 
Vgi  the  gate-to-source  voltage,  and  V*  the  drain-to-source  voltage.  According  to  BSIM4 
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gate  leakage  current  model  including  the  source-drain  partition  [CaoOO],  the  perturbation 
voltage  Vi(y)  may  be  rewritten  as 


V,O0  = 
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where  q is  the  electron  charge,  the  tunneling  barrier  height,  h the  Planck  constant,  tox 
the  gate  oxide  thickness,  mox  the  electron  effective  mass  of  the  conduction  band  in  the 
oxide,  Pigcd  a fitting  parameter  with  a default  value  of  1,  ^ the  mobility,  and  C„x  the  gate 
oxide  capacitance  per  unit  area.  Jg o is  the  direct  tunneling  current  density  with  Vds=  0 and 
is  given  by  [CaoOO] 
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with  a correction  function 
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where  Fox  is  the  electric  field  in  the  oxide  layer,  Ai  a fitting  parameter,  0/,o  the  Si-SiC>2 
barrier  height  ( 3.1  eV  for  electrons  and  4.5  eV  for  holes  ),  and  N an  indicator  of  the 
density  of  tunneling  carriers. 

5.2.3  Noise  Performance  at  Low  Frequencies 

For  the  calculation  of  noise  performance  at  low  frequencies,  the  dimensionless 
noise  coefficients  P and  R are  defined  as  [Puc74]. 


where  |f21|  = gm^and  |Fll|«2#CgJ.  At  low  frequencies,  R is  not  due  to  the  capacitive 

coupling,  but  due  to  the  trap-assisted  tunneling  through  the  gate  oxide  layer.  The  gate-to- 
source  capacitance  Cgs  in  the  linear  region  and  in  the  saturation  region  is  calculated 
according  to  the  following  formula,  respectively 
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where  W is  the  effective  channel  width  and  Vgd  is  the  gate-to-drain  voltage.  The  cross 
correlation  coefficient  C(f)  between  the  drain  and  the  gate  can  be  written  as 
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C(f)  = 


V^(/)-\(/) 


(5.19) 


With  knowledge  of  the  basic  noise  sources,  it  is  possible  to  calculate  the  low 
frequency  noise  performance  using  PLF,  RLF,  and  &F . We  shall  consider  fundamental 
noise  parameters  as  a function  of  the  noise  spectrum  ratio  at  low  frequencies.  The 
minimum  noise  figure  Fff  can  be  obtained  in  terms  of  the  optimum  source  resistance 
RFFopt,  the  noise  resistance  rnLF , the  noise  conductance  gFF , and  the  correlation 
impedance  ZFF  as  follows  [Puc74] 
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where  Rg  and  Rf  are  the  parasitic  gate  and  source  resistance,  respectively,  and  fc  is  the 
intrinsic  cut-off  frequency  (=  gmg  / 2 7iCgs ) . For  Rg  = Rf=  0,  the  minimum  noise  figure  is 
simplified  to  the  form 
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The  low  frequency  noise  figure  can  be  written  in  the  form 
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where  Z5  = Rs  + jXs  is  the  source  impedance,  and  Zlfopt  and  Y^pl  = G^pl  + jB^ipl  are  the 
optimum  external  source  impedance  and  the  optimum  source  admittance  at  low 
frequencies,  respectively.  From  Equations  (5.17),  (5.21),  and  (5.24),  expressing  Y^pl 
for  Rg  = R/=  0 in  terms  of  the  noise  ratio,  we  can  get 


g„  s,r</) 


> mg 


8*  Sff(f) 


C-  =-2 


1-C 


LF 


s>:  (/) 
5,  S"(f) 


(5.25a) 


(5.25b) 


93 


In  addition  to  the  influence  of  the  correlation  on  the  low  frequency  noise  performance, 
the  shot  noise  due  to  the  gate  leakage  current,  S',  (/)  = 2 qlgFFano , has  an  impact  on  the 

noise  behavior  by  modifying  the  noise  parameters,  and  Fff , as  follows  [Dan95] 
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Using  Equations  (5.3),  (5.25),  and  (5.26a),  the 
leakage  current  taken  into  account,  Z*F 
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5.3  Correlation  Noise  Measurements 

Prior  to  noise  correlation  measurements,  all  devices  were  examined  for  the  drain 
and  gate  leakage  current  by  measurements  of  1-V  characteristics.  Figure  5.3  (a)  and  (b) 
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illustrate  Id,g  ( Vas ) in  the  strong  inversion  at  vg*  = 1 -0  V and  Id,g  (Vgs)  in  the  linear  region 
at  Vds  = 0.3  V,  respectively.  The  gate  leakage  current  in  the  direct  tunneling  regime  is  not 
negligible,  and  the  ratio  of  gate  current  to  drain  current  ranges  from  ~ 10  4 to  ~ 10"2  V 
for  nMOSFETs  with  W = 10  yam,  L = 10  /urn,  and  tox  = 2.2  nm.  A dc  analysis  of  all 
devices  was  made  including  two  conductances  ( g g = dl g /dVgf  , gd  = dl d ldV(h ) and  two 

transconductances  ( gmg  = dId/dVgs,  gmd  = dlg  ldVds).  These  parameters  are  used  to 

obtain  the  noise  equivalent  circuit  and  to  extract  the  cross  noise  power  spectrum.  As 
shown  in  Figure  5.3  (c)  and  (d),  the  transconductance  gmd  has  a negative  value  in  the 
linear  region  and  the  lowering  of  the  gate  conductance  gg  points  out  the  peak  of  the 
transconductance  gmg. 

Noise  correlation  measurements  were  performed  using  a HP3582A  double 
channel  spectrum  analyzer  in  the  frequency  range  1 Hz  to  25  kHz.  Brookdeal  5004  low- 
noise  voltage  amplifier  were  used  for  both  the  drain  and  gate  current  noise  measurement. 
The  voltage  noise  at  the  gate  and  the  drain  of  MOS  transistors  in  a common-source 
configuration  are  simultaneously  obtained,  allowing  cross  correlation  noise  acquisition. 
The  drain  current  noise  5,  (/)  and  gate  current  noise  S,  (/)  are  directly  measured  at  the 

drain  and  gate  of  the  device,  and  the  cross  noise  spectrum  S/(/  (/)  is  determined  from 

the  cross  correlation  coefficient  and  noise  measurements.  Figure  5.4  (a)  shows  a typical 
drain,  gate,  and  their  cross  noise  spectrum  measured  at  Id=  9.0x1 0'5  A and  lg  = 6.5x1 0'7 
A.  The  noise  spectra  clearly  show  a 1 If  frequency  dependence  with  the  spectral  power 
factor  y close  to  unity.  As  shown  in  Figure  5.4  (b),  we  observe  that  a partial  correlation 
between  the  drain  and  the  gate  current  noise  sources  appears  and  varies  with  the  current 
ratio.  In  the  lOHz-lkHz  frequency  range,  the  existence  of  correlation  effects  is 


Conductance,  g (A/V).  g (A/V)  Current.  ld(A).  I (A) 
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Drain  Voltage,  Vd#(V)  Gate  Voltage.  Vfl#(V) 

(a)  (b) 


0 0.5  1 1.5 

Gate  Voltage.  V8>(V) 


(c) 


(d) 


Figure  5.3  Measured  drain  and  gate  current  as  a function  of  (a)  drain  voltage  Vds  at  Vgs- 
1.0V,  (b)  gate  voltage  at  VdS  - 0.3  V,  and  measured  conductance  and 
transconductance  as  a function  of  (c)  drain  voltage  Vjs  at  V^.5=  1.0  V,  (d)  gate 
voltage  Vgs  at  V<is  = 0.3  V for  nMOSFETs  with  W = 10  //m,  L=  10  /urn,  and 
t„x  = 2.2  nm 
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Frequency,  f (Hz) 
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(b) 

Figure  5.4  Measured  correlation  noise  characteristics  (a)  Measured  drain,  gate,  and  their 
cross  noise  spectrum  and  (b)  measured  correlation  coefficient  between  drain 
and  gate  as  a function  of  frequency  / for  nMOSFETs  with  W = 10  /urn,  L = 
10  fim,  and  tox=  2.2  nm 
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probably  caused  by  common  noise  sources  which  are  slow  trap  states  available  for  trap- 
assisted  tunneling. 

5.4.  Simulation  Results 
5.4.1  Correlation  Function  Characteristics 

The  variations  of  noise  ratio  S , and  the  correlation  coefficient 

CLF  (/)  with  current  ratio  IJId  are  plotted  in  Figure  5.5  (a)  and  (b),  respectively. 

These  values  are  compared  with  those  calculated  from  Equations  (5.7)  and  (5.1 1)  - (5.16) 
for  nMOSFETs  biased  at  Vgs  = 0.8  - 2.0  V and  Vds  = 0.3  V with  W = 10  ^m,  L = 10  fxm, 
t„x  = 2.2  nm,  /i  = 300  cm2/V-s,  m = 1.2,  VT=  0.2  V,  <pb=?>.\  eV,  and  mox  = 0.4  mo,  and 
Jd  -yg  = 0.5.  Good  agreement  is  obtained  for  the  noise  ratio  at  / = 100  Hz  and  for  the 
correlation  coefficient  at/=  10  Hz  and/=  100  Hz.  In  order  to  investigate  the  correlation 
behavior  in  more  details,  the  dependencies  of  the  cross  correlation  coefficient  CLF  at  low 
frequencies  on  the  current  ratio  Ap,  the  noise  ratio  Sj*  (/)/ S,LF  (/) , the  noise  source 

range  yt/L,  and  frequency / are  shown  in  Figure  5.6.  As  illustrated  in  Figure  5.6  (a),  the 
value  of  correlation  coefficient  is  reduced  as  both  the  current  and  noise  ratio  decrease, 
implying  that  the  correlation  coefficient  can  be  negligibly  smaller  in  a device  with  thicker 
gate  oxide  layer.  Figure  5.6  (b)  indicates  that  the  correlation  coefficient  depends  strongly 
on  the  range  over  which  noise  is  generated  and  traces  the  noise  source  along  the  channel. 
This  is  consistent  with  the  fact  that  the  correlation  increases  with  the  stress  which 
produces  the  percolation  leakage  path  near  the  drain  side  [CheOl,  JonOl].  According  to 
Equations  (5.7)  and  (5.1 1),  the  frequency  dependence  on  the  correlation  coefficient  is 
shown  when  there  is  the  difference  between  the  drain  and  the  gate  frequency  power 


exponent. 
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(a) 


(b) 

Figure  5.5  Comparison  of  measured  data  and  simulation  results  of  (a)  noise  ratio 
between  gate  and  drain  (/=  100  Hz  ) and  (b)  correlation  coefficient  as  a 
function  of  current  ratio  Ig/Ij  for  nMOSFETs  biased  at  Vgs=  0.8  - 2.0  V and 
Vds  = 0.3  V with  W - 10  yum,  L = 10  /urn,  tox  = 2.2  nm,  [x  = 300  cm2/V-s,  m = 
1.2,  VY=  0.2  V,  <f>b=  3.1  eV,  mox  = 0.4  mo,  and  yci-yf,=  0.5 
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i ^ 


Noise  Source  Range,  ^/L  0 0 Frequency,  Log10(f)  (Hz) 

(b) 


Figure  5.6  Characteristics  of  the  cross  correlation  function  CLF.  (a)  dependency  of  CLF 
on  the  current  ratio  and  the  noise  ratio  and  (b)  dependency  of  C'F  on  the 
frequency  and  the  noise  source  range  for  nMOSFETs  biased  at  = 1.0  V, 
Vds  = 0.3  V with  W = 10  n m,  L = 10  /urn,  tox  = 2.2  nm,  [i  = 300  cm  /V-s,  m = 
1.2,  VY=  0.2  V,  (pb=  3.1  eV,  and  mox  = 0.4  m0,  Ap=3.9xl0'4,  yi=0.01  L,  and 
yd-yg=0.5 
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(a) 


(b) 


Figure  5.7  Optimum  source  impedance  (a)  as  a function  of  frequency  and  (b)  as  a 
parameter  of  frequency  for  various  current  ratios  in  nMOSFETs  with  W = 10 
fim,  L = 10  fim,  tox-  2.2  nm,  ld-  0.1mA,  gm=  0.1  mA/V,  FFann=l,  yd=  1.0, 
yg=  1.0,  SldLF=  l.OxlO'16  A2/Hz  ( @/=  1Hz  ),  m=1.2,  y,=0.01L,  and  Slg/S,d 
* SjSmg 
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Figure  5.8  Characteristics  of  the  minimum  noise  figure  F^  . (a)  dependency  of  F^  on 
the  current  ratio  and  frequency  and  (b)  dependency  of  F^fn  on  the  optimum 
source  resistance  and  reactance  for  different  current  ratio  in  nMOSFETs  with 
W=  10  fim,  L - 10 //m,  tox=  2.2  nm,  Id  = 0.1mA,  gm=  0.1  mATV,  FFano=l,  yd 
= 1.0,  yg  = 1.0,  S,dLF=  l.OxlO16  A2/Hz  ( @ /=  1Hz  ),  m=  1.2,  y^O.OlL,  and 

S,g/SId  ~ g gig  mg 
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5.4.2  Noise  Performance  Characteristics 

Figure  5.7  (a)  and  (b)  show  how  the  optimum  source  impedance  varies  with  the 
current  ratio  and  frequency.  Notice  that  the  optimum  source  impedance  decreases  with 

increasing  current  ratio  and  the /_l  dependence  of  the  optimum  source  resistance  changes 

—1/2 

to  / dependence  as  the  gate  current  becomes  larger.  Figure  5.8  (a)  visualizes  the 


dependencies  of  the  minimum  noise  figure  Ft 


LF 

min 


leakage 


on  the  gate  leakage  current  and 


frequency.  Observe  that  the  minimum  noise  figure  increases  with  increasing  gate  leakage 
current  and  clearly  exhibits  l//r  dependence  at  low  frequencies.  As  illustrated  in  Figure 
5.8  (b),  the  variations  with  the  gate  leakage  current  and  frequency  emphasize  the 
important  role  played  by  these  parameters  in  the  optimum  low  frequency  noise 
performance.  Therefore,  the  strong  dependence  of  the  minimum  noise  figure  on  the 
current  ratio  and  frequency  reflects  the  importance  of  minimizing  the  gate  leakage  current 
and  noise  of  ultrathin  oxide  MOSFETs  at  low  frequencies. 


5.5  Conclusion 

An  analytical  model  describing  the  low  frequency  noise  performance  in  ultrathin 
gate  oxide  MOSFETs  was  presented.  For  a comprehensive  evaluation  of  noise 
performance,  the  model  was  formulated  on  the  basis  of  a unified  approach  involving  the 
gate  leakage  current  and  the  cross  correlation  coefficient  between  the  drain  and  the  gate. 
For  the  physical  modeling  of  the  correlation  coefficient  at  low  frequencies,  intrinsic  noise 
sources  were  taken  into  account,  with  the  help  of  the  partition  noise  theory  and  the 
BSIM4  partition  current  model.  It  was  found  that  the  correlation  behavior  is  characterized 
by  the  current  ratio  lgllj,  the  noise  source  range,  and  the  frequency  power  factor.  The 
correlation  characteristics  were  verified  by  correlation  noise  measurements  which  can 
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reveal  the  noise  sources  distributed  along  the  channel.  In  addition,  in  order  to  investigate 
the  noise  performance  characteristics  at  low  frequencies,  the  frequency-dependent  noise 
parameters  PLF,  RLF,  and  CLF  were  introduced  and  the  minimum  noise  figure,  optimum 
source  resistance  and  reactance  were  graphically  presented.  It  was  shown  that  the  details 
of  low  frequency  noise  performance  are  modified  due  to  the  presence  of  the  gate  leakage 
current  and  the  cross  correlation  between  the  input  and  the  output.  The  results  confirmed 
that  the  proposed  model  allows  accurate  prediction  of  the  noise  performance  and 
optimum  designs  of  ultrathin  oxide  MOS  devices  at  low  frequencies.  Furthermore,  the 
model  is  analytically  presented  in  circuit-design  oriented  form  and  therefore  is  suitable  to 
be  implemented  in  circuit  simulators. 


CHAPTER  6 

1 If  NOISE  MODEL  FOR  DRAIN  CURRENT 

6.1  Introduction 

It  is  widely  accept  that  the  shrinkage  of  MOS  devices  will  cause  an  increase  of  the 
low  frequency  noise  level  due  to  the  scaling  law  [Cel99].  Specifically,  the  scaling  of 
device  dimensions  results  in  the  thinning  of  the  oxide  layer,  leading  to  excessive  gate 
leakage  current  [Sch94].  In  addition  to  the  channel  current,  the  gate  leakage  current  in 
ultrathin  gate  oxide  MOS  devices  should  be  considered  as  a low  frequency  noise  source 
[Mom98], 

The  existing  low  frequency  noise  models  are  based  on  the  trapping-detrapping  of 
carriers  by  the  tunneling  assisted  process  [Cel85,  Chr68,  Hun90a,  Hun90b]  or  thermally 
activated  process  [Fan91,  Sur88].  It  was  observed  from  random  telegraph  signals  (RTS) 
in  small  MOSFET  that  switching  rates  have  been  thermally  activated,  often  with  different 
activation  energies  for  capture  and  emission  by  lattice  relaxation  [Kir86,  Ral84]. 
Moreover,  an  inelastic  process  associated  with  energy  loss  during  trap-assisted  tunneling 
was  proposed  to  model  stress-induced  leakage  current  (SILC)  in  ultrathin  oxides  [Lee02, 
Tak99].  Generally,  the  low  frequency  noise  behavior  is  characterized  by  the  trap 
distribution  near  the  Si-SiC>2  interface.  Since  traps  are  distributed  over  space  inside  the 
oxide  as  well  as  in  energy  over  the  bandgap,  the  energy  and  space  distribution  of  trap 
states  should  be  rigorously  modeled  [Cel87],  The  dissimilar  shapes  of  the  energy 
distribution  near  the  conduction  and  valence  band  edge  lead  to  differences  in  the  gate 
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voltage  and  temperature  dependence  of  the  noise  observed  for  n and  pMOSFETs 
[Cha94]. 

As  the  channel  length  and  the  oxide  thickness  are  shorter  and  thinner,  the  effects 
of  the  two-dimensional  electric  field  and  parasitic  resistance  on  the  mobility  and  I-V 
characteristics  become  extremely  important  and  have  been  shown  to  degrade  the 
transconductance  of  devices  [Ong87,  Wu85],  These  effects  may  produce  abnormal  bias 
dependence  of  the  low  frequency  noise  behavior  [KolOl,  KumOO].  For  ultrathin  oxide 
MOS  devices,  the  impact  of  the  quadratic  mobility  degradation  factor  on  the  1 If  noise 
was  investigated  [Mas98].  However,  the  lack  of  comprehensive  noise  models  leads  to 
inaccurate  prediction  of  noise  performance  in  advanced  MOS  devices.  In  addition, 
although  the  process  technology  can  have  significant  influence  on  the  noise 
characteristics,  a process-dependent  noise  model  has  not  yet  emerged  [Hun90a,  Hun90b]. 
Therefore,  in  order  to  design  future  CMOS  devices,  analog  circuits,  and  new  device 
processing  technologies,  an  accurate  modeling  and  characterization  of  noise  in  deep- 
submicron  MOSFETs  with  ultrathin  gate  oxide  is  strongly  needed,  taking  into  account 
effects  of  the  device  miniaturization. 

In  this  Chapter,  we  will  present  a 1 If  drain  current  noise  model  for  deep- 
submicron  MOSFETs  with  ultrathin  oxide,  operating  in  both  the  linear  and  saturation 
regions.  This  model  incorporates  a more  realistic  trap  distribution,  and  the  effects  of  the 
two-dimensional  electric  fields  and  the  parasitic  resistance.  Physically,  the  model  is  based 
on  a tunneling  process  of  carriers  between  the  Si-Si02  interface  and  oxide  trap  states, 
followed  by  a thermally  activated  process  due  to  lattice  relaxation.  The  dependence  of  the 
drain  current  noise  on  the  doping  profile  is  described  through  the  bulk  charge  parameter. 
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In  Section  6.2,  an  improved  analytical  mobility  and  current  model  for  advanced  MOS 
devices  including  the  parasitic  source-drain  series  resistance  is  developed  with  the 
dependence  of  the  effective  field  on  the  doping  profile.  Section  6.3  deals  with  an 
effective  time  constant  model  based  on  the  tunneling  assisted-thermally  activated 
(TATA)  processes.  Also,  we  formulate  a new  model  for  trap  distribution  over  energy  and 
space.  In  this  formulation,  the  trap  density  inside  the  bandgap  increases  as  energy 
becomes  close  to  the  conduction  and  valence  band  edges,  and  the  trap  density  versus  trap 
distance  is  an  exponentially  decreasing  distribution.  On  the  basis  of  TATA  process,  a 
new  noise  model  is  formulated  by  taking  the  device  structure  and  processes  into 
consideration  through  a realistic  trap  distribution,  a bulk  charge  parameter,  and  parasitic 
resistance.  Moreover,  this  noise  model  accounts  for  the  gate-bias  dependence  of  the 
frequency  exponent  of  the  noise  power  spectrum.  In  Section  6.4,  in  order  to  demonstrate 
the  applicability  of  the  developed  model,  we  illustrate  how  the  device  characterization 
and  parameter  extraction  procedure  can  be  applied  to  MOS  devices  featuring  short 
channel  length  and  ultrathin  oxide  thickness.  In  Section  6.5,  the  developed  drain  current 
and  noise  models  are  compared  and  verified  with  experimental  results  for  nMOSFETs. 
Finally,  some  conclusions  are  presented  in  Section  6.6. 

6.2  Mobility  and  I-V  Model 

The  electron  mobility  in  the  surface  inversion  layer  along  the  channel  direction  y 
depends  on  normal  field  E^iy)  and  lateral  field  Eiiy)  and  can  be  expressed  as  [Han82, 
Wu85] 


MnO 

1 + &\EN(y)\  + /32\EL(y)\ 


(6.1a) 
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with 


Qb (y)  = -p£si(lNsub[V(y)  + i//s-v J,  ys  = 2vth  In 

l ni 


(6.1c) 


where  funo  is  the  maximum  electron  mobility  in  the  inversion  layer,  (3\  is  an  empirical 
constant,  /?2=/Wv.s/  with  the  scattering-limited  velocity  vsi,  Cox  is  the  gate  oxide 
capacitance  per  unit  area,  esi  is  the  silicon  dielectric  permittivity,  Vgs  is  the  internal  gate- 
source  voltage,  VFB  is  the  flat  band  voltage,  \j/s  is  the  surface  inversion  potential,  V(y)  is 
the  potential  distribution  along  the  channel,  £ is  a bulk  charge  parameter  through  which 
the  dependence  of  the  effective  field  on  the  doping  profile  is  described  [Kru88],  Qb(y)  is 
the  silicon  bulk  charge,  Nsub  is  the  substrate  doping  density,  n,  is  the  intrinsic  carrier 
density,  v,h  = kTlq  is  the  thermal  voltage,  and  q is  the  electron  charge.  The  drain  current  of 
nMOSFETs  operating  in  strong  inversion  is  written  by  [Wu85] 


(6.2a) 


with 


Qn(y)  = Cox[Vgs-VT(Vds)-V(y)] 


(6.2b) 


(6.2c) 


where  Qn(y ) is  the  electron  charge  density  per  unit  area,  W is  the  effective  channel  width, 
and  Vds  is  the  internal  drain-source  voltage.  Putting  Equation  (6.1c)  into  Equation  (6.2c) 
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and  expanding  the  threshold  voltage  around  V&’  = 0,  we  can  obtain 


VT  (Yds  ) _ ^T0  k builds 

, (6'3) 
Q 4 (V's-vJC,,  * 

Assuming  that  the  parasitic  source  and  drain  resistance  are  characterized  by  a 
resistor  R,  the  relationship  between  the  external  voltage  (Vgs,Vds)  and  the  internal  voltage 
(Vgi  ,Vds ) can  be  written  as  [Cho92] 

Vv=V„-R-It,  v;=V*-2  RId  (6.4) 

Using  Equations  (6.1)-(6.4)  and  integrating  along  the  channel  length  L,  the  drain  current 
Id.iin  for  short-channel  MOSFETs  biased  in  the  linear  region  is  then  approximated  as 
[Cho92,  Wu85] 


r 

Vfb  + Vs  + ^ Nsub  (y/ s — vth ) 


I 


d,lin 


w „ tygs  VtO  Kul^dsVds 

LMneff  “ 1 + Vds/ECL 


(6.5a) 


with 


Yn,eff  ~ 
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1 + ax  (v^  V^0  kbulkVds )+  a2  tygs  vT  khulk  Vds  )2 


(6.5b) 


2v  l w 

Ec  > kbulk  kbulk , ctx  = (X0  + ——2R/in0Co 

"n.eff  ^ L 


(6.5c) 


where  /. in,eff  is  the  effective  mobility,  a\  and  are  the  mobility  degradation  parameters, 
vSat  is  the  carrier  saturation  velocity,  and  Ec  is  the  critical  saturation  field.  Note  that 
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Equation  (6.5b)  is  an  analytical  effective  mobility  model  which  can  be  applied  to 
ultrathin  gate  oxide  nMOSFETs  [Mas98,  McL95,  Ong87],  In  the  saturation  mode  of 
operation,  the  lateral  electric  field  reaches  the  critical  field  Ec  and  carriers  travel  with  the 
saturated  velocity  vsa,  near  the  drain  section  of  the  length  Lsat.  Using  Equations  (6.4)  and 
(6.5a),  and  considering  channel-length  shortening  effects,  the  drain  current  in  the 
saturation  region  is  obtained  as  [Cho92] 


C fyg*  ^r°  kbulkVdsa,  ha, 

d.sal  ~ r^n.eff^'ox  / __  \ 


(6.6a) 


with 


(6.6c) 


vdsat=vgs-vFB-¥s- 


bulk  y dsat 


(6.6c) 


(6.6d) 


where  Vjsat  is  the  drain  saturation  voltage  and  Xj  is  the  drain  junction  depth. 
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6.3  Formulation  of  Noise  Model 
6.3.1  Effective  Time  Constant  Model 

It  was  found  that  the  capture  or  emission  time  constant  is  energy-activated  and 
could  be  either  enhanced  or  reduced  by  the  applied  gate  voltage.  The  traditional 
McWhorter’s  approach  based  on  the  tunneling  process  may  have  problems  on  this 
account  [And90,  McW57,  Tew94].  For  ultrathin  oxide  MOSFETs,  we  suggest  a model 
based  on  a tunneling  process  of  carriers,  followed  by  the  thermally  activated  process  due 
to  lattice  relaxation  multiphonon  capture-emission  in  oxide  trap  states,  as  shown  in 
Figure  6.1  [Lee02,  Tew94].  This  is  called  the  tunneling  assisted-thermally  activated 
(TATA)  model.  To  model  the  switching  response,  we  assume  that  both  the  capture  and 
emission  time  constants  may  be  represented  by  the  sum  of  two  time  constants,  one  being 
tunneling  assisted  and  the  other  being  thermally  activated.  Specifically,  we  write 


1 


(6.7) 


where  Tc,e  is  a capture  or  emission  time  constant.  The  thermally  activated  time  constant 
Tth  is  given  by 


T,h  = exP 


kT 


= Ttho  exp(^rt  F) 


(6.8) 


where  t,m  is  the  thermally  activated  attempt  time,  E is  the  activation  energy,  and  X,h  is  an 
exponential  factor  of  thermal  activation  process.  From  a Taylor’s  series  approximation 
for  the  tunneling  time  constant  Ttun  in  a direct  tunneling  regime,  r,un  can  be  simplified  to 
[And90,  Tew94] 


Ill 


•X 


(6.9a) 


with 


ox 


ox 


(6.9b) 


where  Tlun0  is  the  tunneling  attempt  time  («  10'10  s),  A,un  is  the  exponential  factor  of 
tunneling  process,  tox  is  the  gate  oxide  thickness,  mox  is  the  effective  mass  of  tunneling 
earners  ( 0.40mo  for  electrons  and  0.32mo  for  holes  ) , (f>b  is  the  tunneling  barrier  height 
( 3.1  eV  for  electrons  and  4.5  eV  for  holes ),  h is  the  reduced  Planck  constant,  Fox  is  the 

electric  field  in  the  oxide  layer,  Vox  is  the  oxide  voltage,  and  x is  the  distance  of  oxide  trap 
located  from  the  Si-SiC>2  interface.  Using  Equations  (6.7),  (6.8),  and  (6.9),  the  effective 
time  constant  can  be  rearranged  in  strong  inversion  as 


where  tc o>e0  is  a capture  or  emission  attempt  time  constant,  Xlh,c  and  Xlh,e  are  exponential 
factors  of  the  capture  and  emission  thermal  activation  process,  respectively,  and  Xtunx 
and  Xlunx  are  exponential  factors  of  the  capture  and  emission  tunneling  process, 
respectively.  The  rate  of  the  capture-emission  process  relaxing  to  an  equilibrium  state  can 
be  defined  in  terms  of  its  effective  time  constant  Tejf  which  is  given  by 


r 


Tc,e (E,x)  = Tc0e0 (E,x)- exp  Ath(cx)E  + qAlun(ce)  1- 


•x 


(6.10) 


r 


effO  ' eXP  ^-th.eff  E + 1 


•X 


(6.11a) 
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Figure  6. 1 Energy  band  diagram  of  the  tunneling  assisted  - thermally  activated  (TATA) 
process 


with 


Tc0(E,x)-Te0(E,X ) 
Tc0(E,x)  + Te0(E,x)’ 


^ th,eff  — ^ih.c 


+ Ah.e  < Aun.eff  ~ \ tun.c  + K un 


(6.11b) 


where  reffo  is  an  effective  attempt  time  constant,  and  A ,h,ejf  and  \lun,eff  are  exponential 
factors  of  the  thermal  activation  and  tunneling  process,  respectively. 


6.3.2  Trap  Distribution  Model 

A nonuniform  trap  distribution  over  a distance  in  the  oxide  near  the  Si-SiC>2 
interface  can  also  cause  a plateau  at  low  frequencies  that  is  similar  to  generation- 
recombination  noise  [Fan86].  Besides  nonuniformity  in  the  trap  distribution  in  the  oxide, 
band-bending  due  to  a gate  bias,  a work  function,  or  oxide  charges  pulls  down  the  traps, 
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causing  the  carriers  to  encounter  more  traps  [Cel87,  Sur86].  Therefore,  the  trap  charge 
distribution  is  expressed  over  energy  E and  distance  x in  the  oxide  as  follows 


( V„  - VT()  } 

Nj(E,x)  = Nt0  - exp 

f,|E-E,|  + 

• X 

l *ox  y 

(6.12) 


where  Et  is  the  intrinsic  Fermi  level,  and  £2,  and  £3  are  exponential  factors 
characterizing  the  distribution  of  traps  in  the  oxide. 

It  is  worthwhile  to  consider  the  practical  aspect  of  trap  distribution  over  energy  in 
a MOS  structure.  In  spite  of  the  fact  that  Equation  (6.12)  has  the  capability  of  probing  the 
trap  distribution,  this  expression  may  be  insufficient  to  obtain  a reasonable  fit  to  the 
experimental  results.  We  propose  a new  model  for  trap  charge  distribution  over  energy  at 
the  Si-Si02  interface.  In  the  new  model,  the  conduction  and  valence  band  tail  has  the 
form  of  a sum  of  two  exponentials,  respectively 


nTc(E)  = nTc0 
nTv(E)  = nTv0 


(!-£)'  exp 


(1  ~ ' exp 
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e~e} 


C 

*re  0 y 

+ be ' exP 

1 

t 

e-eA 

r E-Ev 

+ L -exp 

1 vo  y 


kT 


vi  y 


(6.13a) 


(6.13b) 


where  Ec  and  Ev  are  the  conduction  and  valence  energy  band,  respectively,  and  Ntco,  £c, 
Tco,  and  Tci  are  the  trap  parameters  for  the  conduction  band  tail,  and  NTv0,  £v,  TM,  and  Tvl 
are  the  trap  parameters  for  the  valence  band  tail.  However,  this  form  results  in  a 
complicated  formula  for  the  noise  spectrum.  Simplifying  Equations  (6.13a)  and  (6.13b), 
and  including  the  spatial  distribution  of  traps,  we  can  obtain  a more  realistic  model  for 
the  trap  charge  distribution 
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(l-£)exp 


E-E  „ V„  -VT, 


N^-M(E,x)  = NT0 


kT 


- + <f2^ -x  + ^x 


+(1-0 


exp 


E-E.,  „ V.  -VTl 


kT. 


■ + £2  ~ —x  + £3x 


+ 2^c 


(6.14) 


Figure  6.2  shows  the  trap  density  distribution  as  a function  of  energy  and  trap  distance 
from  the  Si-Si02  interface.  In  this  figure,  the  trap  density  inside  the  bandgap  increases  as 
the  energy  becomes  closer  to  the  conduction  and  valence  band  edges  due  to  tailing 
effects.  The  trap  density  versus  trap  distance  is  an  exponentially  decreasing  function 
away  from  the  interface. 


-'8V 


Energy,  E ( eV)  0 0 Trap  Distance,  x ( nm  ) 


Figure  6.2  Trap  density  distribution  as  a function  of  energy  and  trap  distance  from  the 
Si-Si02  interface.  ( Vgs  = 0.7  V,  Vro  = 0.45  V,  = 6.25xl018  cm' 13eV\  Tc 
= 300  K,  Tv  = 600  K,  & = £v=  0.2,  & = 10  V1,  £ = -l.OxlO7  cm1 ) 
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6.3.3  Low  Frequency  Noise  Model 

The  tunneling  assisted-thermally  activated  (TATA)  time  constant  model  may  be 
used  to  explain  1 If  noise  in  MOSFETs,  where  the  carriers  are  assumed  to  be  tunneling  to 
traps  located  in  the  oxide  layer,  and  to  be  followed  by  a thermal  activation  process  in  the 
oxide  trap  states.  For  traps  far  from  the  Fermi  level,  the  traps  contribute  an  insignificant 
amount  of  noise  because  most  of  the  time  the  traps  are  either  empty  or  filled.  Integrating 
over  the  noise  components  of  the  occupancy  of  traps  within  the  elemental  volume 
AxAyAz,  the  noise  power  spectral  density  due  to  all  the  traps  in  the  gate  oxide  is  obtained 
[Hun90a,  Hun90b] 


Teff(E,x ) 


W C f f ■ ^E,wyM-mr-+  Wx) 


dxdzdE  (6.15) 


where  Nj{Ej;,y,z)  is  the  charge  density  of  traps  and  f,  is  the  probability  that  a trap  is 
filled  given  by  the  Fermi  factor 


/,=- 


1 


1 + exp 


f E-E^ 


fit 


kT 


(6.16) 


with  E the  trap  energy  level  and  Efr  the  trap  quasi-Fermi  level.  Note  that  the  term  /,( 1 -/,) 
is  a peaked  function  about  the  electron  quasi-Fermi  level  and  reflects  the  fact  that  only 
the  traps  near  the  electron  quasi-Fermi  level  contribute  to  fluctuations  or  noise.  Assuming 
that  the  energy  and  spatial  distribution  of  traps  is  given  by  Equation  (6.12)  and/,(l-/,) 
behaves  like  a delta  function,  the  integral  in  Equation  (6.15)  can  be  approximated  as  an 
analytical  expression 
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s„T(f)  = - 


Ay-W-kT ■ Nto  exp[^,|E^  - E,  ] 

q^-tun.eff 

0 

1 

1 

■ [iKTeffO  exp U,H,effEfn  ‘ ‘ COS 

, 40b/(]  , 

- (6.17a) 


with 


Y = 1 + - 


tiK-vj 


(l^wn,eJftox 
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V„-Vto 

WJq 


'tun,eff 


1 — 


V ' y y 


V^L 

Wbiq 


(6.17b) 


V ' y y 


where  the  slope  of  the  spectrum  is  directly  related  to  the  gate  bias. 

If  we  consider  a section  of  channel  with  width  W and  length  Ay,  fluctuations  in 
the  trapped  oxide  charge  will  induce  correlated  fluctuations  in  the  channel  carrier  number 
and  mobility.  The  resulting  fractional  change  in  the  local  drain  current  can  be  represented 
as  [Hun90a,  Hun90b] 


yAQn  SAQt  [leg 


8AQTj 


■8AQt 


(6.18) 


where  A Qn  = Q„W  Ay,  A<2r=  QtW  Ay,  and  Qn  and  QT  are  total  charge  of  channel  carriers 
and  occupied  traps  per  unit  area.  The  sign  of  the  mobility  term  is  chosen  either  positive 
for  acceptor-like  traps  or  negative  for  donor-like  traps.  The  ratio  R0  = SAQn  / 8AQT  in 
Equation  (6.18)  describes  the  efficiency  of  the  drain  current  modulation  corresponding  to 
a variation  of  the  trapped  charge.  The  coefficient  Rq  depends  on  the  bias  and  on  the 
spatial  location  of  the  trapping  center  within  the  oxide  layer  and  is  given  by  [Pac99] 


R 


o 


c. 


Kx  ~ x, 


+ (Cdep  + Ci,  + Cinv  ) 


'l  + ^-' 

v Kx  -x,j 


*ax  ~ X'  Qn 
Kx  Qn  + Qn 


(6.19a) 
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with 


(6.19b) 


(6.19c) 


where  Cinv  is  the  capacitance  of  the  inversion  layer,  Ninv  is  the  areal  density  of  electrons  in 
the  inversion  layer,  is  the  depletion-layer  capacitance,  C„  is  the  interface  trap 
capacitance,  and  the  distance  of  charged  trap  from  the  Si-SiC>2  interface,  xt,  is  obtained  by 
rewriting  Equation  (6.11a)  in  terms  of  the  characteristic  frequency.  On  the  basis  of 
Matthiessen’s  rule  [Jay89],  we  can  readily  find  that 


where  S(x,,Niriv)  is  a scattering  rate  which  is  a function  of  the  distance  x,  of  the  trapping 
state  [Jay 89]  and  of  the  carrier  density  Ninv  [Hun90c],  and  is  strongly  affected  by  channel 
quantization  [Pac99].  In  order  to  properly  account  for  these  effects  in  ultrathin  oxide 
MOSFETs,  we  employ  an  empirical  model  for  the  scattering  coefficient  S(xt,Ninv)  in  a 
linear  approximation  as  [Jay89,  VanOO] 


(6.20) 


soWrm«  * V^T 


1 


(6.21) 
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where  rTmx  ~ (£ox  + £si  + C*p  + Q ) *s  the  maximum  distance  from  the  scattering 

charge  for  which  the  effect  of  image  charges  can  be  neglected  and  So  is  the  scattering 
constant  in  cm/C1/2Vs.  Substituting  Equations  (6.19)  and  (6.20)  into  Equations  (6.18), 
and  using  Equations  (6.2b)  and  (6.3),  the  power  spectral  density  of  the  local  current 
fluctuations  is  obtained  as 


S„(f)  = (AId)2  = 


QnWAy 


«o  ± S/t* 


.~\2 


<1  J 


S A2r  (/) 


with 


Qn  = <lNinV  = Cox [Vg5  - VT0  - (1  - kbulk  )V(y)] 
S,QT(f)  = q2-S^T(f) 


(6.22a) 


(6.22b) 

(6.22c) 


where  S4Gr(/)  is  the  power  spectral  density  of  fluctuations  in  the  total  charge  of 

occupied  traps  over  the  area  W Ay.  Performing  the  integration  over  the  channel  length, 
the  total  noise  power  spectral  density  of  the  drain  current  is  obtained  as 


StJ(y,f)Ay  dy 


qkT-1] 


(6.23) 
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Combining  kbulk  = 0.5  - kbulk  and  Cox  = £ox/tox,  and  replacing  the  integration  over  y by 


integration  over  Qn  yield 
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Substituting  Equations  (6.19a)  and  (6.21)  into  Equation  (6.24a)  and  integrating, 
get  an  analytical  expression  for  the  drain  current  noise  power  in  the  linear  region 
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In  the  saturation  region,  the  total  noise  power  can  be  expressed  as  the  sum  of  the  noise 
components  due  to  the  triode  region,  L-Lsat,  and  the  pinch-off  region,  Lsat.  Using  Equation 
(6.23)  and  integrating  over  each  region,  we  can  obtain 
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In  order  to  incorporate  a more  realistic  model  for  the  trap  distribution,  we  may  use 
Equation  (6.14)  describing  the  conduction  and  valence  band  tails.  In  the  linear  region,  it 
is  straightforward  to  show  that  a noise  contribution  due  to  traps  near  the  conduction  band, 
Sfd  (/),  a noise  contribution  due  to  traps  near  the  valence  band,  Sj  (/),  and  a noise 

contribution  due  to  traps  near  the  midgap,  S"  (/) , are  given  by 
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where  7/c  and  iVv  are  the  effective  densities  of  states  in  the  conduction  and  valence  bands, 
respectively. 


6.4  Experiments  and  Parameter  Extraction 

The  devices  used  are  n+  poly-Si  gate  nMOSFETs  which  were  fabricated  using  a 
remote  plasma  nitrided  oxide  (RPNO)  process  [Lee02,  NicOO].  About  1. 7-4.3  nm  of  gate 
oxide  was  grown  with  both  furnace  and  rapid-thermal  based  processes.  Nitridation  was 
performed  at  room-temperature  by  exposing  the  gate  oxide  to  a short,  high-density, 
remote  helicon-based  nitrogen  discharge.  This  was  then  followed  by  post-nitridation 
annealing  in  an  inert  or  low  partial-pressure  oxygen  ambient.  Four  terminal  devices  were 
fabricated  with  standard  processing  techniques  on  p-Si  substrate  with  n+  poly-Si  gates, 
using  LOCOS  isolation  to  define  the  active  regions. 

Parameter  extraction  is  very  important  for  the  characterization,  modeling,  and 
simulation  of  MOS  devices  and  circuits.  The  device  parameters  should  be  extracted  and 
determined  so  that  the  model  accurately  predicts  experimental  data.  It  is  necessary  to  be 
careful  in  ultrathin  oxide  MOSFETs  because  the  mobility  degrades  at  high  gate  voltages 
resulting  in  a negative  transconductance,  as  shown  in  Figure  6.3  (a)  for  the  equivalent 
oxide  thickness,  tox,eq  = 2.2  nm.  The  equivalent  oxide  thickness  tox,eq  can  be  calculated  as 
tox,eq  = (ZoJzsiON)  tsio\,  where  the  nitrided  oxide  and  oxide  dielectric  constant,  EsiON  and 
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£r«,  are  assumed  to  be  5.7  £o  and  3.9  Eo,  respectively.  The  quadratic  dependence  of  gate 
voltage  through  the  coefficient  a2  in  Equation  (6.5b)  accounts  for  the  mobility 

degradation  and  the  negative  transconductance.  Figure  6.3  (b)  illustrates  the  typical  plots 

1/2 

of  Ids  / gm  as  a function  of  the  gate  voltage  V The  straight  line  extrapolations  from 
high  gate  voltages  intercept  the  x-axis,  allowing  the  determination  of  the  threshold 
voltage  Vro  for  different  drain  voltages  Vds.  It  is  possible  to  calculate  the  bulk  charge 
parameter  £ in  Equation  (6.3)  from  the  extracted  values  of  Vm  = 0.45  V.  In  addition, 
obtaining  £ enables  the  determination  of  khulk  and  kbulk  expressed  in  Equations  (6.3)  and 
(6.5c),  respectively.  It  was  found  that  £ and  kbuik  are  1.13  and  0.625,  respectively.  The 
low  field  mobility  parameter  /i„o  and  the  degradation  parameter  a2  are  extracted  from 
Figure  6.3  (c).  In  high  gate  voltage  regions,  the  curves  show  a linear  relationship, 
allowing  the  determination  of  /uno  from  slope  and  a2  from  extrapolation  to  the  y-axis.  The 
extracted  values  of  Uno  and  a2  were  550  cm2V'1s'1  and  0.33  V'2.  It  is  worth  noting  that  the 
presence  of  a2  for  ultrathin  oxide  leads  to  a great  impact  of  mobility  degradation  on  the 
current-voltage  and  noise  characteristics  [Mas98].  Finally,  the  parasitic  resistance  R and 
the  mobility  degradation  parameter  or0  can  be  obtained  from  plots  of  a0  vs  Vgs-Vjx)  for 
several  values  of  Vds ■ Figure  6.3  (d)  is  plotted  substituting  the  extracted  parameters  into 
Equations  (6.5a)  - (6.5c).  By  plotting  the  constant  values  in  high  gate  voltage  regions  as 
a function  of  V*  and  applying  a least  squares  fit,  R = 6 H and  a0  = 0.01  V"1  are  obtained. 

In  ultrathin  oxide  MOSFETs,  the  two  key  noise  parameters  Ro  and  S are  very 
sensitive  to  the  trap  location  and  the  inversion  charge  density.  Thus,  the  accurate 
prediction  of  noise  behavior  depends  on  a satisfactory  estimate  of  these  parameters. 
Rearranging  Equations  (6.22a)  and  (6.23)  for  y = 1,  Efi,  =£„  and  Vgs-Vn  < 4 (pbJq,  we  get 
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(a) 


(b) 


(c) 


(d) 


Figure  6.3  Static  characteristics  and  parameter  extraction  for  ultrathin  oxide  nMOSFET 
with  tox,eq  = 2.2  nm  ( tSioN  = 3.2  nm  ),  W = 10  fim,  and  L = 0.14  yum  at  T = 
300  K.  (a)  transconductance  gm  and  drain  current  /*  vs  Vgs  for  different  V*, 
(b)  plots  of  Ids  / gm1'2  vs  Vgs  determining  Vto  and  £ (c)  plots  of 
d(l/Ids)/dVgs  vs  l/(VgS-VTf  determining  fi„o  with  slope  and  a2  with  intercept, 
and  (d)  plots  of  a i vs  V„  determining  R and  ao. 


Ratio,  R0  Noise  Power,  S J°/Id  (Hr'°) 
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(a) 


(b) 


Figure  6.4  Noise  parameter  extraction  (a)  plots  of  Sid'2  / Id  vs  1 /Ninv  for  the 
determination  of  Ro  with  slope  and  5 with  intercept,  (b)  plots  of  Su'2  / Id  vs 
1 /Ni„v  for  the  determination  of  Ro  and  S as  a function  of  Ninv  from  the  least 
squares  fit,  (c)  plots  of  Ro  vs  Ni/lv  for  extracted  data  and  model  ( x,  = 0.5  tox, 
Nu  = 2.0xl012  cm'2eV"'  ) , and  (d)  plots  of  S vs  2Vinv  for  extracted  data  and 
model  ( x,  = 0.5  tox,  So  = 8.0xl09  cm/C1/2Vs  ). 
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As  illustrated  in  Figure  6.4  (a),  Ro  and  S are  extracted  from  plots  of  the  variation 
of  the  normalized  drain  current  Sid/Idm  with  1/Ninv.  R0  is  determined  from  the  slope  and 
S is  obtained  by  extrapolating  the  measured  data  in  low  density  regions  to  the  y-axis.  R0 
and  S were  found  to  be  0.61  [Pac99]  and  l.lxlO'16  Vs  [VanOO],  respectively.  In  order  to 
extract  Ro  and  S as  a function  of  Ninv,  a least  squares  fit  is  applied  to  the  measured  data. 
The  trap  distance  x,  and  the  scattering  parameter  So  determined  from  Figure  6.4  (b)  were 
found  to  be  0.50  [Lee02]  and  8.0xl09  cm/C1/2  Vs  [VanOO],  respectively.  Using  the 
extracted  values  and  the  theoretical  models,  Equations  (6.19a)  and  (6.21),  R0  and  S are 
plotted  as  a function  of  Ninv,  as  shown  in  Figure  6.4  (c)  and  Figure  6.4  (d). 


6.5.1  Static  Characteristics 

Prior  to  the  noise  model  verification,  the  static  characteristics  are  analyzed  and 
validated  by  using  the  extracted  values  of  mobility  and  I-V  parameters.  In  ultrathin  oxide 
MOSFETs,  the  mobility  dependence  on  the  gate  voltage  should  be  accurately  modeled, 
since  the  strong  normal  field  produces  abnormal  I-V  and  noise  characteristics.  Figure  6.5 
(a)  shows  the  comparison  between  the  measured  data  and  the  theoretical  model  equation 
(6.5b)  of  effective  inversion  carrier  mobility  with  respect  to  the  gate  voltage  at  low  drain 
voltage.  The  effective  mobility  was  calculated  by  using  the  set  of  parameters  obtained 
from  the  extraction  procedure  ( Vm  = 0.45  V,  kbuik  = 0.625,  fin0  = 550  cm2  VV1,  R = 6VL, 
ao  = 0.01  V'1,  and  ct2  = 0.33  V'2  for  nMOSFETs  with  t0Xteq  = 2.2  nm,  W = 10  //m,  L = 


6.5  Simulation  Results  and  Discussion 


0.14  /urn ). 
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(a) 


(b) 

Figure  6.5  Comparison  of  the  mobility  and  I-V  models  to  measured  data  for  nMOSFETs 
with  tax, eg  = 2.2  nm,  W = 10  /urn,  L = 0.14  //m,  Vto  = 0.45  V,  kbuik  = 0.625,  /i„0 
= 550  cm2  VV\  R = 6 fl,  a0  = 0.01  V1,  or2  = 0.33  V'2 , Nsub  = 6xl017  cm'3, 
vSat  = 8xl06  cm/s,  and  Xj  = 0.1  ftm  at  T = 300  K.  (a)  effective  mobility  finiejf 
versus  gate  voltage  Vgs  at  Vds  = 50  mV  and  (b)  Id-Vds  characteristics  for  Vgs  = 
0.6,  0.8,  1.0  V 
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It  is  observed  that  the  quadratic  mobility  parameter  «2  enables  the  reproduction  of 
the  strong  mobility  degradation  at  high  gate  voltage.  The  I-V  characteristics  of  both  the 
measured  data  and  the  simulation  results  are  plotted  in  Figure  6.5  (b).  The  simulation  was 
carried  out  with  the  set  of  mobility  parameters  and  the  process  parameters  Nsub  = 6xl017 
cm"3  and  Xj  = 0.1  fim  for  vsat  = 8xl06  cm/s.  The  results  clearly  show  an  excellent 
agreement  between  the  theoretical  prediction  of  the  model  and  the  measured  data  in  both 
the  linear  and  saturation  regions.  This  is  due  to  the  fact  that  the  model  incorporates 
several  effects  of  mobility  degradation,  parasitic  resistance,  carrier  velocity  saturation, 
and  channel  length  modulation. 

6.5.2  Ilf7  Noise  Characteristics 

Noise  measurements  were  performed  using  a HP3582A  spectrum  analyzer  in  the 
frequency  range  of  1 Hz  to  25  kHz.  The  Brookdeal  5004  low  noise  voltage  amplifier  was 
used  for  the  drain  current  noise.  Figure  6.6  shows  a typical  drain  current  noise  spectrum 
for  nMOSFETs  biased  at  Vds  = 0.2  V and  = 1.0  V with  tox,eq  = 2.2  nm,  W = 10  fim, 
and  L = 0.14  /rm.  As  shown  in  Figure  6.6  (a)  and  (b),  the  noise  model  was  verified  by  the 
simulation  results  for  the  trap  distribution  Ntcm(Ejc)  and  Nt'(Ejc)  with  a fitting  parameter 
of  Tc  and  £i,  respectively.  It  can  be  seen  that  the  frequency  dependence  of  1 If  noise 
spectrum  is  described  accurately  by  the  trap  distribution  model  as  a function  of  energy 
and  trap  distance.  This  means  that  the  difference  in  the  noise  behavior  observed  by 
various  process  technologies  can  be  well  explained  by  adjusting  the  trap-related  model 
parameters  of  Equations  (6.12)  and  (6.14). 

The  bias  dependence  of  the  drain  current  noise  power  of  the  measured  data  and 
simulation  results  is  plotted  in  Figure  6.7.  As  illustrated  in  Figure  6.7  (a),  the  noise  power 
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in  the  linear  region  shows  a more  gradual  dependence  on  the  gate  voltage  than  in 
saturation.  This  behavior  can  be  understood  by  different  degradation  of  mobility  with  the 
gate  voltage.  Figure  6.7  (b)  shows  the  variations  of  the  noise  power  with  the  drain 
voltage  for  different  gate  voltage  values.  A good  agreement  is  observed  in  both  the  linear 
and  saturation  regions. 

In  order  to  investigate  the  effect  of  the  gate  oxide  thickness  on  the  noise  behavior 
in  MOSFETs  with  direct  tunneling  gate  leakage  current,  the  variations  of  the  normalized 
drain  current  noise  with  the  oxide  thickness  are  plotted  in  Figure  6.8  (a).  It  was  confirmed 
that  the  normalized  noise  power  decreases  as  the  oxide  thickness  becomes  thinner 
[KolOl,  Mom98].  This  is  attributed  to  the  fact  that  the  noise  model  is  based  on  the  carrier 
capture-emission  process  via  the  oxide  trap.  In  addition,  the  influence  of  the  doping 
profile  on  the  \/f  noise  behavior  is  investigated.  The  bulk  charge  parameter  £ depends  on 
the  doping  profile  in  the  near-surface  region.  For  pMOS  devices  fabricated  with 
uniformly  doping,  the  twin-tub,  and  the  n-well  technologies,  the  values  of  £are  1.0,  1.4, 
and  1.7,  respectively  [Kru88].  For  nMOS  devices  from  the  twin-tub  and  the  n-well 
technologies,  the  values  of  £ are  1.0  and  0.6,  respectively  [Kru88].  Figure  6.8  (b)  shows 
the  normalized  noise  power  as  a function  of  the  bulk  charge  parameter  £ for  different 
drain  voltages  in  nMOSFETs.  The  noise  power  decreases  with  decreasing  £,  implying 
that  a technology  process  with  lower  value  of  £ is  optimum  for  noise  reduction. 

The  impact  of  the  quadratic  mobility  parameter  a2  becomes  important  in  thinner 
oxide  thickness.  The  variations  of  the  drain  current  noise  power  with  the  drain  current  as 
a parameter  of  or2  are  illustrated  in  Figure  6.9.  The  noise  power  decreases  with  increasing 


Drain  Current  Noise  Specral  Density,  Sid  ( A2/Hz  ) 
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Figure  6.6  Frequency  dependence  of  noise  power  Sid  of  measured  data  and  simulation 
results  for  (a)  the  trap  distribution  NtC’  M(Ejc)  with  Nm  = 6.25xl018  cm"3eV'\ 
£•  = £v  = 0.2,  £2  = 10  V1,  and  £3  = -lxlO7  cm'1  and  (b)  the  trap  distribution 
Nt‘  {Ej c)  with  Ntd  = 2.5xl017  cm'Vv1,  £>  = 10  V1,  and  & = -lxlO7  cm1  for 
nMOSFETs  with  tOXieq=  2.2  nm,  W=  10  /im,  L = 0.14  //m,  Vn  = 0.45  V,  kbuik 
= 0.625,  So  = 8.0xl09  cmJCm  Vs,  fa  = 8xl0'7  cm/V,  and#  = 3xl0'6  cm/V 
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Figure  6.7  Bias  dependence  of  noise  power  Sid  of  measured  data  and  simulation  results  at 
/=  10  Hz  as  a function  of  (a)  the  gate  voltage  V*.  with  the  drain  voltage  V* 
as  a parameter  and  (b)  the  drain  voltage  V*  with  the  gate  voltage  Vgs  as  a 
parameter  for  nMOSFETs  with  tox,eq  = 2.2  nm,  IV  = 10  pm,  L = 0.14  pm, 

= 0.45  V,  = 0.625,  and  So  = 8.0xl09  cm/C1/2  Vs.  The  trap  distribution  is 
assumed  to  be  NtC’m(Ejc). 
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(a) 


(b) 

Figure  6.8  Characteristics  of  the  normalized  drain  current  noise,  (a)  Variations  of  the 
normalized  drain  current  noise  SiJll  with  the  equivalent  oxide  thickness  tox,eq 
and  (b)  variations  of  the  normalized  drain  current  noise  SiJId  with  the  bulk 
charge  parameter  £ in  the  linear  and  saturation,  for  nMOSFETs  with  W = 10 
Hm,  L = 0.14  //m,  Vto  = 0.45  V,  kbuik  = 0.625,  and  So  = 8.0xl09  cm/C1/2  Vs. 
The  trap  distribution  is  assumed  to  be  NtCM{Ejc). 
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Figure  6.9  Variations  of  the  drain  current  noise  power  Su  with  the  drain  current  Id  as  a 
parameter  of  the  quadratic  mobility  degradation  factor  a2  at/=  10  Hz  and 
V*  = 0.2  V. 


Drain  Current  Noise  Specral  Density,  Su  ( A2/Hz  ) 
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Figure  6.10  Variations  of  the  drain  current  noise  power  Sid  with  the  normalized  trap 
distance  x,ltOXteq  as  a parameter  of  the  scattering  coefficient  So  for  the 
number  fluctuation  ( Nta  ' A term  ) and  the  correlated  fluctuation  ( 
NtbCM  B*+  NtcCM  C*  term  ) at  Vds  = 0.2  V,  Vgi  = 1.0  V,  and /=  10  Hz, 
assuming  that  the  trap  distance  xt  is  independent  of  frequency/. 
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<y 2 at  higher  drain  current.  Therefore,  the  accurate  extraction  of  aj  is  essential  to  the  noise 
modeling  in  ultrathin  oxide  MOS  devices  [Mas98]. 

In  order  to  explain  a broader  set  of  noise  data,  a unified  1//  model  has  been 
developed  [Hun90a,  Hun90b,  Jay89,  VanOO].  This  model  is  based  on  the  fact  that 
fluctuations  in  the  occupancy  of  the  oxide  traps  induces  correlated  fluctuations  in  the 
carrier  number  and  surface  mobility.  As  shown  in  Figure  6.10,  the  comparison  of  the 
number  fluctuation,  NtaC,m  A*  term  in  Equations  (6.27),  and  the  correlated  mobility 
fluctuation,  Ntb  ' B +Njc  ’ C term  in  Equation  (6.27),  is  demonstrated  as  a function  of 
the  trap  distance  in  the  linear  region.  The  number  fluctuation  contribution  gets  smaller  for 
the  scattering  coefficient  So  = 8.0x10s  - 8.0xl09  cm/C1/2Vs,  and  gets  larger  for  So  = 
8.0xl07  cm/C1/2Vs,  exhibiting  downward  variations  with  the  trap  distance.  This  is  mainly 
due  to  the  fact  that  the  mobility  fluctuation  depends  on  the  scattering  rate  S(xt,Ni„v ) which 
is  larger  for  charge  scattering  sites  nearer  to  the  interface. 

6.6  Conclusion 

A 1 If  drain  current  noise  model  for  advanced  MOSFETs  operating  in  both  the 
linear  and  saturation  regions  was  presented.  Physically,  the  model  is  based  on  the 
tunneling  assisted-thermally  activated  process  of  carriers,  incorporating  a more  realistic 
description  of  trap  distribution.  The  model  includes  the  effects  of  the  two-dimensional 
electric  fields,  the  parasitic  resistance,  the  doping  profile,  and  the  mobility  degradation  on 
the  noise  behavior.  Especially,  the  impacts  of  the  trap  location  and  the  inversion  carrier 
density  on  both  the  trapping  efficiency  ratio  and  the  scattering  rate  have  been 
quantitatively  estimated  and  taken  into  account  for  ultrathin  oxide  MOSFETs.  It  was 
shown  that  the  noise  characteristics  are  accurately  predicted  by  the  proposed  model 
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imposing  no  physical  limitation  on  the  trap  distribution  and  trapping-detrapping  process. 
By  using  the  extracted  model  parameters,  the  bias  dependence  of  the  noise  behavior  was 
well  described  in  both  the  linear  and  saturation  regions.  Furthermore,  it  was  recognized 
from  the  simulation  results  that  the  effects  of  process  technologies  on  the  noise 
characteristics  can  be  analyzed  by  the  bulk  charge  parameter,  the  quadratic  mobility 
degradation  factor,  parasitic  resistance,  and  trap-related  parameters.  The  results 
confirmed  that  the  noise  model  is  applicable  to  deep-submicron  MOSFETs  with  ultrathin 
gate  oxide.  The  model  is  analytically  represented  and  therefore  suitable  for 
implementation  in  circuit  simulators. 


CHAPTER  7 

CHARGE-BASED  QUANTUM  CORRECTION  MODEL 

7.1  Introduction 

With  the  introduction  of  MOS  devices  aggressively  scaled  to  deep  submicron 
channel  length,  a ultrathin  gate  oxide  and  high  doping  concentration  are  required  to 
increase  the  drive  current  and  to  minimize  undesirable  short  channel  effects.  These 
design  rules  result  in  strong  transverse  electric  fields  at  the  Si-Si02  interface  even  near 
the  threshold  of  inversion.  This  gives  rise  to  significant  band  bending  of  the  energy 
bands,  leading  to  the  confinement  of  electron  motion  in  the  direction  perpendicular  to  the 
interface.  In  these  conditions,  electrons  in  the  inversion  layer  form  actually  a quasi  two- 
dimensional  electron  gas  (2-DEG)  with  quantized  energy  levels  within  the  channel 
[And82,  Ste72].  These  quantum  mechanical  (QM)  effects  have  significant  impact  on  the 
performance  of  nanoscale  MOSFETs  and  thus  the  classical  models  are  inadequate  for 
state-of-the-art  MOS  structures  [Bac83,  Dor92,  Ohk90]. 

The  quantum  nature  of  electrons  in  the  inversion  layer  is  typically  studied  by 
solving  Schrodinger’s  and  Poisson’s  equations  self-consistently  [Mog86,  Ste72],  and  also 
considering  the  simplified  models  used  in  circuit  simulators  [Dor94],  A computationally 
efficient  model,  which  utilizes  the  first  three  subbands  in  a quantized  inversion  layer,  has 
been  described  for  predicting  the  QM  effects  on  the  inversion  layer  charge  density  and 
distribution  [Har92].  The  enhancement  in  the  centroid  and  the  finite  thickness  of  the 
inversion  layer  are  one  of  the  main  QM  effects  on  MOSFET  behavior  [Liu99,  Lop97]. 
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This  reduces  the  inversion  layer  charge  density  and  the  total  gate  oxide  capacitance 
[LarOl,  Nak98,  OmuOO,  Tak95],  compared  to  classical  models  without  QM  effects,  and 
causes  the  threshold  voltage  shift  effect  [Jan94,  MaOO,  Vas99]. 

For  a compact  circuit  simulation  model  of  direct  tunneling  of  electrons  form 
nMOS  inversion  layer,  a variational  approximation  is  employed  to  calculate  the  surface 
potential  and  the  energy  of  the  first  quantized  level  in  the  inversion  layer  without  the 
requirement  for  iterative  solution  [CleOl],  and  an  efficient  model  is  demonstrated  by 
using  both  QM  calculations  for  the  substrate  and  a modified  WKB  approximation  for  the 
transmission  probability  [Yan99].  It  is  shown  from  the  quantum  charge  sheet  model  that 
QM  effects  on  the  drain  current  are  expected  to  be  due  to  QM  effects  on  gradients  in  both 
the  inversion  charge  and  in  the  surface  potential  [Ip98].  Based  on  an  analytical  charge- 
control  model,  the  drain  current  model  fully  comprising  QM  effects  as  well  as  short 
channel  effects  is  developed  [MaOl]. 

Besides  QM  effects  on  the  static  characteristics,  QM  effects  have  been  observed 
in  the  trap  and  noise  behavior  [CelOO,  LukOO,  Pac99,  Pal97,  Sie92].  For  heavily  doped 
MOSFETs,  extracted  interface  trap  capture  cross  sections  are  found  to  be  an  order  of 
magnitude  larger  than  values  from  the  classical  theory  [Sie92].  In  order  to  investigate 
QM  effects  on  the  capture  and  emission  times  of  random  telegraph  signal  (RTS)  noise, 
the  interaction  between  an  oxide  trap  and  an  electron  in  the  subbands  of  the  inversion 
layer  is  theoretically  modeled  and  analyzed  in  a Si-SiC>2  interface  [CelOO,  LukOO,  Pal97]. 
Regarding  the  1 //  noise  models  of  drain  current,  numerical  models  are  employed  to 
account  for  QM  impact  on  the  extracted  value  of  the  oxide  trap  density  [Pac99]. 
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However,  there  are  no  analytical  1 If  and  thermal  noise  models,  which  incorporate  QM 
effects  and  are  applicable  to  circuit  simulations. 

In  this  Chapter,  we  will  present  a charge-based  quantum  correction  model  for 
current  and  noise  in  nanoscale  MOSFETs.  Section  7.2  describes  the  quantization  effects 
on  the  inversion  charge  parameters,  such  as  the  surface  potential,  the  threshold  voltage, 
the  inversion  charge  density,  and  gate  capacitance  model.  In  Section  7.3,  on  the  basis  of 
the  quantum  correction  model  of  the  inversion  layer  charge,  a gate  leakage  current  and  1 If 
noise  model  are  formulated.  Section  7.4  gives  the  analytical  expressions  for  the  drain 
current  and  noise  model  accounting  for  quantum  and  short  channel  effects.  The  quantized 
charge  parameters  are  calculated  from  the  lowest  energy  level,  and  the  comparison 
between  the  classical  and  quantum  models  is  shown  in  Section  7.5.  Finally,  Section  7.6 
presents  some  conclusions. 

7.2  Quantum  Correction  for  Compact  Modeling 

7.2.1  Schrodinger-Poisson’s  Equations 

The  high  channel  doping  levels  and  large  transverse  electric  fields  result  in  the 
confinement  of  electron  motion  in  the  x direction  perpendicular  to  the  interface,  so  that 
the  conduction  band  within  the  channel  is  split  into  discrete  subbands.  Based  on  the 
effective  mass  approximation  and  a parabolic  band  structure,  the  envelope  wave  function 
£/  for  the  ith  valley  and y'th  subband  satisfies  1-D  Schrodinger’s  equation  [Ans82,  Ste72] 
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+ ql//(x) 


Cij(x)  = EijCij(x ) 
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where  q is  the  electron  charge,  h is  the  reduced  Planck  constant,  mXi  is  the  effective  mass 

perpendicular  to  the  surface  ( mx\  =0.916  mo,  mX2  = 0.19  mo ),  Eg  is  the  energy  state  of  the 
zth  valley  and  y'th  subband,  and  is  the  electrostatic  potential  which  is  the  solution  of 
1-D  Poisson’s  equation 

d [p(x)  - n(x)  + Nj  ( x ) - N~  (*)]  (7.2) 

ax  e 

where  esi  is  the  silicon  dielectric  constant,  n(x)  is  the  electron  density,  p(x)  is  the  hole 
density,  and  N/(x)  and  Na(x)  are  the  donor  and  acceptor  concentrations,  respectively.  By 
solving  Equations  (7.1)  and  (7.2)  self-consistently  under  appropriate  boundary 
conditions,  the  2-D  electron  density  of  the  inversion  layer  is  obtained  by  summing  over 
all  subbands  as  follows  [Har92,  Ste72] 
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where  k is  the  Boltzmann’s  constant,  T is  the  temperature,  g,  is  the  degeneracy  of  the 
energy  subband  ( g\  = 2,  g2  = 4 ),  m^  is  the  density-of-states  effective  mass  ( m^i  = mX2, 
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nidi  = {mx\mx2)m  ),  E/ is  the  Fermi  energy  level,  <pf-  v,h  ln(  NSUb/rii ) is  the  Fermi  potential, 
if/s  is  the  surface  potential,  Eg  is  the  energy  bandgap,  Fs  is  the  surface  electric  field,  v,h  = 
kT/q  is  the  thermal  voltage,  Nsub  is  the  substrate  doping  density,  n,  is  the  intrinsic  carrier 
density,  and  Nij  is  the  sheet  charge  density  for  the y'th  subband  from  the  ith  valley. 

7.2.2  Surface  Potential  and  Threshold  Voltage  Model 

Due  to  the  nonzero  value  of  the  quantized  charge  centroid,  the  surface  potential 
determined  by  quantum  effects  is  greater  than  predicted  by  the  classical  model  [Dor94]. 
Multiplying  the  Poisson  equation  (7.2)  by  x and  integrating  it,  we  get 
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where  Xdep  is  the  depletion  layer  width,  Qn  is  the  total  electron  charge  density,  and  xnQM  is 
the  average  centroid  position  from  the  interface.  Therefore,  it  is  observed  from  Equation 
(7.5)  that  the  enhancement  of  the  surface  potential  is  given  by  [Lop97] 
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Due  to  the  quantum-mechanical  behavior  of  inversion  layer  electrons,  more  band  bending 
is  required  to  achieve  the  same  inversion  charge  as  the  classical  value,  resulting  in  the 
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increase  in  the  threshold  voltage.  The  threshold  voltage  VY  can  be  easily  calculated  by 
including  an  additional  quantum  shift  term  as  follows  [MaOO,  MaOl] 


Vtqm  =VT  + AVj 
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For  the  polysilicon-Si02-p-Si  MOS  structure,  the  relationship  between  the  gate-to-source 
voltage  Vgs  and  the  surface  potential  tf/s  is  classically  described  by  [CleOl,  LarOl,  Mas02] 


V„=VFB+Vr,+V^+^l+^-  (7.10) 

^ OX  '-'ox 

with 

where  Vfb  is  the  flat-band  voltage,  Vpoiy  is  the  polysilicon  voltage  drop,  Npoiy  is  the 
polysilicon  doping  concentration,  Cox  is  the  oxide  capacitance  per  unit  area,  eox  is  the 
oxide  dielectric  constant,  tox  is  the  oxide  thickness,  Qs  is  the  total  charge  in  the 
semiconductor  per  unit  area,  and  the  trapped  charge  in  the  oxide,  Qit,  is  included  to 
account  for  the  trapping-detrapping  effects  at  the  Si-Si02  interface.  In  order  to  obtain  the 
quantum  model,  it  is  necessary  to  replace  the  classical  Qs  with  the  quantized  charge 
density  as  follows 
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where  Qd  is  the  depletion  charge  density  per  unit  area. 

7.2.3  Inversion  Charge  and  Capacitance  Model 

As  the  channel  length  and  oxide  thickness  continue  to  be  scaled  to  nanometer 
range  in  current  CMOS  devices,  the  quantization  effects  should  be  taken  into  account  for 
an  accurate  inversion  charge  model.  By  using  Equations  (7.5),  (7.7),  and  (7.12),  an 
expression  for  the  inversion  layer  charge  is  obtained  as 
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where  the  denominator  toxQM  of  Equation  (7.14)  is  termed  the  electrical  oxide  thickness. 
Including  the  parasitic  source  and  drain  resistance  Rs,d  and  incorporating  the  drain-to- 
source  voltage  V*,  the  channel  charge  density  is  expressed  by 


er  (y)  = cr  ■(v„-l„R,J,-Vrm -2 kukVy(y)) 


(7.15) 


where  Id  is  the  drain  current,  Vy  (y)  is  the  channel  potential  at  any  given  point  y along  the 
channel  and  kbuik  is  the  bulk  charge  parameter. 
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Since  the  quantum  effects  lead  to  an  increase  in  the  inversion  layer  thickness,  the 
inversion  capacitance  decreases  and  thus,  the  total  gate  capacitance  is  reduced.  The 
inversion  charge  capacitance  is  given  by  [Liu99,  MaOl] 


with 


1 _J_+XCM 


C?M  C.  £, 


n si 


xQM  = 


^x-n(x)dx  _ 5.047 + 0.77 IE 
J n(x)dx 


0.113  + E 


Eeff=- 


2 1 


f f 

W lh  In 

1 + exp 

l V 

y -y  -V 

Vgs  v poly  VT 


w 


JJ 


+ 2(vr -V„-2^,) 


\£si<lNsub 


4 0fC, 


/ ^ ox 


(7.16) 


(7.17a) 


(7.17b) 

(7.17c) 


where  C„  is  the  classical  inversion  capacitance,  X<2M  is  the  inversion  layer  thickness,  Eeg 
is  the  effective  electric  field  in  inversion  layer  in  MV/cm,  and  nsw  is  the  subthreshold 
swing  parameter.  For  very  high  polysilicon  doping  levels  ( Npoiy  > lxlO19  cm'3  ),  no 
quantum  behavior  has  to  be  considered  to  determine  the  surface  potential  at  the 
polysilicon-Si02  interface,  because  the  all  the  subbands  are  squeezed  out  of  the  potential 
well  [LarOl].  According  to  the  classical  theory,  the  poly  silicon  capacitance  is  calculated 
by  differentiating  the  polysilicon  charge  density  Qpoiy  as  follows  [LarOl] 
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7.3  Quantum  Effects  on  Gate  Leakage  Current  and  Noise  Model 
7.3.1  Gate  Leakage  Current  Model 

In  the  ultrathin  gate  oxide  regime  ( < 3nm  ),  the  classical  models  without 
including  quantum  effects  of  the  inversion  layer  overestimate  the  value  of  the  gate 
leakage  current  and  noise,  because  the  peak  charge  density  is  situated  at  a distance  ( « 1 
nm  ) away  from  the  Si-Si02  interface.  This  leads  to  an  enhancement  in  the  oxide 
thickness  and  a reduction  in  the  oxide  voltage.  Therefore,  the  quantum  effects  need  to  be 
taken  into  account  by  using  the  correction  for  the  oxide  thickness  and  the  oxide  voltage, 
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The  BSIM4  gate  tunneling  current  model  can  predict  all  the  significant  direct  tunneling 
current  components  which  are  electron  tunneling  from  the  conduction  band  in  the 
substrate  to  the  gate,  electron  tunneling  from  the  valence  band  in  the  substrate  to  the  gate, 
and  hole  tunneling  from  the  valence  band  in  the  gate  to  the  substrate.  Taking  into  account 
quantum  effects  of  the  oxide  thickness  and  the  oxide  voltage,  the  BSIM4  tunneling 
current  model  becomes  [Lee02] 
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where  h is  the  Planck  constant,  mox  is  the  electron  effective  mass  of  the  conduction  band 
in  the  oxide,  A is  a fitting  parameter  for  three  current  components,  (pbo  is  the  Si-SiC>2 
barrier  height  ( 3.1  eV  for  electrons  and  4.5  eV  for  holes  ),  and  (f>b  is  the  tunneling  barrier 
height. 

7.3.2  Gate  Current  Noise  Model 

One  possible  mechanism  for  1 If  noise  of  the  gate  leakage  current  is  for  traps  to 
affect  locally  the  barrier  height  or  shape  due  to  thermal  noise  of  the  frequency-dependent 
conductance  of  the  oxide  slow  traps  [Lee03].  This  fluctuation  in  turn  modulates  the 
tunneling  transmission,  and  causes  fluctuations  in  the  tunneling  current.  The  quantum 
effects  for  the  spectral  noise  density  are  accounted  for  by  the  correction  of  the  oxide 
voltage  Vox  and  the  electron  subband  energy  Ey , and  are  described  by  [Lee03] 
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where  m±  is  the  transversal  electron  effective  mass,  E±  is  the  transversal  energy 
component,  Ag  is  the  gate  tunneling  area,  GPQM  ( / ) and  CPQM  ( f ) are  the  frequency- 
dependent  parallel  conductance  and  capacitance  of  oxide  traps,  respectively.  Assuming 
that  the  oxide  trap  density  is  spatially  distributed  over  a physical  distance  d into  the 
oxide,  the  electrical  trap  distance  <fiM  becomes  d + {e ox  / e si  in  the  quantum 

approach,  leading  to 
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7.4  Quantum  Effects  on  Drain  Current  and  Noise  Model 
7.4.1  Mobility  and  Drain  Current  Model 

As  the  gate  oxide  thickness  and  the  substrate  doping  level  are  thinner  and  higher, 
the  normal  electric  field  increases  significantly  and  the  influence  of  the  quantized 
inversion  layer  become  important.  A correct  and  simple  model  to  account  for  these 
effects  is  required  to  predict  the  mobility  and  I-V  characteristics  of  nanoscale  MOSFETs. 
The  quadratic  mobility  degradation  factor  is  necessary  for  an  effective  mobility  model  of 
ultrathin  gate  oxide  nMOSFETs,  and  an  analytical  mobility  model  is  given  by 

[Mas98] 
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where  Hn o is  the  low-field  electron  mobility,  and  afM  and  are  the  mobility 
degradation  factor.  Introducing  an  effective  drain-to- source  voltage  V£Meff , and 

considering  the  channel  length  shortening  effects  and  parasitic  resistance  Rs,d,  the  drain 
current  for  short-channel  MOSFETs  biased  in  the  linear  and  saturation  region  is  then 
obtained  as  [Cho92,  MaOl] 
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where  W is  the  channel  width,  L is  the  channel  length,  Ec  is  the  lateral  critical  electric 
field  at  which  carriers  travel  with  the  saturated  velocity  vsat  near  the  drain  section  of  the 
length  Lsal,  and  A is  a small  quantity  to  adjust  the  smoothness  at  the  transition  region  near 
Vdsar  [Cho92]. 

7.4.2  RTS  ( Random  Telegraph  Signal ) Noise  Model 

If  we  consider  a section  of  channel  with  width  W and  length  Ay,  fluctuations  in 
the  trapped  oxide  charge  will  induce  correlated  fluctuations  in  the  channel  carrier  number 
and  mobility.  The  resulting  fractional  change  in  the  local  drain  current  can  be  represented 


as  [Hun90b] 
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(7.28) 


where  A<2n  = Q„  W A y,  &Qt=  Qt  W Ay,  and  Qn  and  Qt  are  total  charge  of  channel  carriers 
and  occupied  traps  per  unit  area.  The  sign  of  the  mobility  term  is  chosen  either  positive 
for  acceptor-like  traps  or  negative  for  donor-like  traps.  The  ratio  R0  = ~SAQn  / SkQT  in 

Equation  (7.28)  describes  the  efficiency  of  the  drain  current  modulation  corresponding  to 
a variation  of  the  trapped  charge.  Considering  the  effects  of  the  quantization  and 
polysilicon  depletion,  the  coefficient  R%M  depends  on  the  surface  potential  through 
charge  capacitance  and  is  given  by  [Hun90b] 


where  C„  is  the  interface  trap  capacitance.  On  the  basis  of  Matthiessen’s  rule  [Jay89],  we 
can  readily  find  that 


(7.29) 
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where  S QM  {xfM , Q®M ) is  a scattering  rate  which  is  a function  of  the  distance 
xfM  (=  xt  +(£ox/  £st  )xfM  ) of  the  trapping  state  [Jay89]  and  of  the  carrier  charge  density 


qQm  [vanQ0],  ancj  js  strongly  affected  by  channel  quantization  [Pac99].  In  order  to 
properly  account  for  these  effects  in  ultrathin  oxide  MOSFETs,  we  employ  an  empirical 
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model  for  the  scattering  coefficient  SQM  [xfM  ,Q®M ) in  linear  approximation  as  [Jay89, 
VanOO] 
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where  rf"  ~(£0JC  + £„-)/(c®M  + C,.,)  is  the  maximum  distance  from  the  scattering 

charge  for  which  the  effect  of  image  charges  can  be  neglected  and  So  is  the  scattering 
constant  in  cm/C  Vs.  Assuming  that  the  carrier  distribution  is  uniform  along  the 
channel,  the  RTS  amplitude  noise  becomes  [CelOO] 
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7.4.3  Low  and  High  Frequency  Noise  Model 

According  to  the  number  fluctuation  theory,  the  total  power  spectral  density  of  the 
drain  current  noise  at  low  frequencies  is  obtained  by  calculating  the  fluctuations  in  the 
number  of  occupied  traps  and  performing  the  integration  over  the  channel  length. 
Assuming  that  the  oxide  trap  density  N,(E)  is  uniformly  distributed  over  the  space  into 
the  oxide  and  the  trapping  time  constant  is  an  exponentially  increasing  function  of  the 
trap  distance  xt,  the  total  1// noise  power  of  the  drain  current  in  the  linear  region  ( Vd  < 
Vdsat ) is  given  by  [Hun90b] 
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with 
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where  ytun  « 108  cm'1  is  the  exponential  factor  of  tunneling  process,  N,(Efn)  is  the  oxide 
trap  density  around  the  electron  quasi-Fermi  level  E^,  and  and  are  carrier 

charge  densities  for  RSfd  = 0 at  the  source  and  drain  ends  of  the  channel,  respectively.  In 
the  saturation  region,  the  total  1 If  noise  power  is  obtained  by  including  the  noise 
contribution  arising  from  the  pinch-off  region  Lsat,  and  is  represented  as  [Hun90b] 
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It  has  been  shown  that  the  high  frequency  noise  originates  from  thermal  fluctuation  in  the 
channel  and  is  expresses  as  [Wan94] 


C”'=4*7'-f-0>,,,vJ  (7.38) 

where  Qn  is  the  total  inversion  layer  charge  density.  For  the  thermal  noise  modeling  of 
short  channel  devices,  the  total  inversion  charge  density  is  derived  by  considering  the 
velocity  saturation  of  carriers  due  to  the  high  lateral  electric  field,  and  is  turned  out  to  be 


[Wan94] 
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Substituting  Equation  (7.39)  into  Equation  (7.38),  the  thermal  noise  spectral  densities  of 
drain  current  in  the  linear  and  saturation  region  are  given  by 


For  the  purpose  of  developing  a model  applicable  to  circuit  simulation,  it  is  very 
important  to  extend  classical  simulators  to  include  QM  effects  with  computational 
efficiency.  For  the  quantized  energy  levels  of  inversion  layers,  a simple  and  approximate 
model  can  be  obtained  by  considering  the  dominant  contribution  of  carriers  in  the  lowest 
energy  subbands  E\o  [CleOl,  Dor94].  Figure  7.1  shows  the  quantized  charge  parameters 
calculated  from  the  lowest  energy  level  E\q  for  the  substrate  doping  concentrations,  Nsub  = 
3xl017  cm'3,  6xl017  cm'3,  and  9xl017  cm"3.  As  the  gate-source  voltage  increases,  the 
sheet  carrier  density  NnQM  in  the  energy  level  E\q  increases  and  the  inversion  layer 
centroid  xnQM  moves  more  closely  to  the  Si-Si02  interface  for  each  doping  concentrations, 
as  shown  in  Figure  7.1  (a)  and  (b),  respectively.  These  results  are  confirmed  by  numerical 
results  and  other  models  [CleOl,  Lop97,  Ohk90],  The  QM  effects  on  the  surface  potential 
and  the  threshold  voltage  are  depicted  in  Figure  7.2.  Figure  7.2  (a)  clearly  illustrates  that, 


(7.40a) 


(7.40b) 
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at  higher  gate-source  voltage  and  thinner  oxide  thickness,  the  transverse  electric  field 
increases  appreciably  and  thus  the  energy  gap  between  the  bottom  of  the  conduction  band 
and  the  lowest  energy  level  is  enhanced,  leading  to  noticeable  QM  effects  [CleOl,  Dor94, 
MaOl].  It  is  obvious  from  Equation  (7.9)  that  the  threshold  voltage  shift  due  to  QM 
effects  becomes  more  pronounced  as  the  oxide  thickness  and  the  substrate  doping 
concentration  increase,  as  shown  in  Figure  7.2  (b).  Figure  7.3  (a)  and  (b)  show  the 
inversion  charge  density  and  capacitance  as  a function  of  the  gate-source  voltage 
according  to  the  classical  and  QM  models,  respectively.  The  deviation  between  the 
classical  and  QM  models  for  the  inversion  charge  density  increases  significantly  for 
higher  doping  concentration,  and  for  the  inversion  layer  capacitance,  the  results  coincide 
in  weak  inversion  and  deviate  in  strong  inversion  [Bac83]. 

By  considering  the  QM  effects  on  the  oxide  voltage  and  the  oxide  thickness,  the 
gate  leakage  tunneling  current  and  noise  model  can  be  formulated.  The  oxide  voltage 
due  to  QM  effects  is  calculated  using  the  additional  surface  potential  to  achieve  the  same 
inversion  charge  as  the  classical  value,  as  shown  in  Figure  7.4  (a).  For  ultrathin  oxide 
MOS  structure,  the  oxide  thickness  shift  due  to  QM  effects  has  a greater  impact  on  the 
gate  tunneling  current  than  the  oxide  voltage  shift  due  to  QM  effects.  Therefore,  the 
extraction  procedure  of  the  electrical  oxide  thickness  accounting  for  QM  effects  is 
required  to  accurately  predict  the  tunneling  current  through  the  gate  oxide.  The  use  of  the 
physical  oxide  thickness  ignoring  QM  effects  leads  to  erroneous  results  in  the  calculation 
of  the  gate  tunneling  current,  as  shown  in  Figure  7.4  (b).  Quantitatively,  it  is  found  from 
Equation  (7.22)  that  the  exact  knowledge  of  the  oxide  thickness  is  also  essential  to  the 
evaluation  of  the  gate  current  noise  model.  As  illustrated  in  Figure  7.5  (a),  the  deviation 
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(a) 


(b) 


Figure  7.1 


Quantized  charge  parameters  calculated  from  the  lowest  energy  level  £10  for 
three  substrate  doping  concentrations,  (a)  sheet  carrier  density  NnQM  as  a 


function  of  the  gate-source  voltage  V„s 
position  from  the  Si-SiCh  interface,  xnQ  , 


and  (b)  inversion  layer  centroid 
as  a function  of  NnQM  ( Vfb  = -0.9 


V,  F=300K  ) 
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(a) 


Figure  7.2  Quantum  effects  on  (a)  the  surface  potential  i jjs  versus  the  gate-source  voltage 
Vv  and  (b)  the  additional  threshold  voltage  shift  A VtQM  versus  the  substrate 
doping  density  Nsub  for  tox  = 2.5  and  4.5  nm.  ( Vfb  = -0.9  V,  Eij  = £10,  Nit  = 
2xlOn  cm'2,  Npoiy  = 2xl020  cm'3,  T=  300  K) 
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Gate-Source  Voltage,  Vgs  (V) 


(a) 


(b) 

Figure  7.3  Quantum  effects  on  (a)  the  inversion  charge  density  Qno  versus  the  gate- 
source  voltage  Vgs  and  (b)  the  inversion  layer  capacitance  Cn  versus  Vgs  for 
Nsub  = 3xl017  cm'3  and  3xl018  cm'3.  ( tox  = 2.5  nm,  VFB  = -0.9  V,  Eg  = Eio, 
Nit  = 2xlOu  cm'2,  Npoiy  = 2xl020  cm'3,  T = 300  K ) 
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(a) 


(b) 

Figure  7.4  Quantum  effects  on  (a)  the  oxide  voltage  Vox  versus  the  gate-source  voltage 
Vgs  and  (b)  the  gate  leakage  current  density  Jg  versus  V„s.  ( tox  = 2.5  nm,  Nsub 
= 3xl017  cm'3,  VFb  = -0.9  V,  Eq  = El0t  Nit  = 2xlOn  cm  , Npoty  = 2xl020  cm'3, 
7=300  K) 
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(a) 


(b) 

Figure  7.5  Quantum  effects  on  (a)  the  normalized  1// noise  of  gate  leakage  current  Sig/Ig 
versus  the  gate-source  voltage  for  two  oxide  thickness  values  and  (b)  the 
normalized  1 If  noise  of  gate  leakage  current  Sig/Ig  versus  the  frequency  /.  ( 
Nsub  = 3xl017  cm'3,  VFB  = -0.9  V,  Eg  = El0,  Nit  = 2xlOu  cm'2,  Npoly  = 2xl020 
cm'3,  T=  300  K ) 
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(a) 


(b) 

Figure  7.6  Quantum  effects  on  (a)  the  electron  effective  mobility  versus  the  gate- 
source  voltage  Vgs  for  two  oxide  thickness  values  and  (b)  the  total  drain 
current  Id, total  versus  the  drain-source  voltage  Vds  for  two  gate-source  voltages 
( Nsub  = 3xl017  cm'3,  VFB  = -0.9  V,  Ey  = Eio , vsal  = 8xl06  cm/s,  Rs,d  = 10  H, 
hulk  = 0.625,  T = 300  K ) 
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(b) 

Figure  7.7  Quantum  effects  on  (a)  the  number  fluctuation  factor  Ro  versus  the  gate- 
source  voltage  with  the  poly-silicon  depletion  effects  and  (b)  the 
scattering  coefficient  S versus  the  gate-source  voltage  Vgs  for  two  substrate 
doping  densities  ( tox  = 2.5  nm,  Vfb  = -0.9  V,  = £10 , T = 300  K ) 
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Figure  7.8  Quantum  effects  on  the  fractional  RTS  amplitude  A ljld  versus  the  surface 
potential  ips  ( tox  = 2.5  nm,  NSUb  = 3xl017  cm'3,  Vfb  = -0.9  V,  Etj  = E\o , So 
=8xl09  cm/C1/2Vs,  kbuik  = 0.625,  T = 300  K ) 
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(a) 


(b) 

Figure  7.9  Quantum  effects  on  (a)  the  1 If  noise  power  of  drain  current  versus  the  drain- 
source  voltage  Vds  and  (b)  the  thermal  noise  power  of  drain  current  versus  the 
drain-source  voltage  Vds  for  two  gate-source  voltages  ( Nsub  = 3xl017  cm'3, 
Vfb  = -0.9  V,  Eij  = Eio , vsat  = 8xl06  cm/s,  Rs,d  = 10  fl,  kbulk  = 0.625,  T=  300 
K ) 
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Figure  7.10  Quantum  effects  on  the  total  noise  power  of  drain  current  versus  the 
frequency  / for  two  substrate  doping  concentrations  ( VFB  = -0.9  V,  Eij  = 
£io,  vsal  = 8xl06  cm/s,  Rsj  = 10  fl,  kbuik  = 0.625,  T = 300  K ) 
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between  the  classical  and  QM  model  increases  with  decreasing  oxide  thickness.  Figure 
7.5  (b)  clearly  shows  a Ilf  frequency  dependence  with  the  power  factor  y close  to  unity. 

On  the  basis  of  the  inversion  charge  model  including  QM  effects,  the  drain 
current  and  noise  model  is  evaluated  for  ultrathin  oxide  MOS  device  ( tox  = 1.5~2.5  nm  ). 
As  shown  in  Figure  7.6  (a),  the  QM  effects  on  the  effective  mobility  become  noticeable 
with  thinner  oxide  thickness  and  higher  gate-source  voltage.  Figure  7.6  (b)  illustrates  that 
QM  models  give  lower  drain  current  values  with  respect  to  the  classical  model  and  QM 
effects  become  more  significant  in  the  saturation  region.  Prior  to  the  evaluation  of  1 If 
noise  model,  it  is  necessary  to  investigate  QM  effects  on  the  RTS  fluctuations  which  are 
the  origin  of  1 If  noise.  The  QM  and  poly  silicon  depletion  effects  on  the  number 
fluctuation  factor  R0  versus  the  gate-source  voltage  are  plotted  in  Figure  7.7  (a).  The  QM 
and  polysilicon  depletion  effects  are  appreciably  enhanced  in  the  strong  inversion.  Figure 
7.7  (b)  shows  QM  effects  on  the  scattering  coefficient  5 as  a function  of  the  gate-source 
voltage.  For  higher  substrate  doping  concentration,  QM  effects  become  more  noticeable. 
Based  on  number  and  mobility  fluctuations,  the  RTS  amplitude  versus  the  surface 
potential  is  evaluated  using  the  QM  correction  charge  model,  as  shown  in  Figure  7.8. 
Note  that  in  the  case  of  number  fluctuations,  the  classical  and  QM  models  exactly 
coincide  in  the  strong  inversion  and  somewhat  deviate  in  the  weak  inversion.  For  the 
evaluation  of  1// drain  current  noise,  it  is  assumed  that  the  number  fluctuation  appreciably 
contributes  to  the  total  current  noise.  As  illustrated  in  Figure  7.9  (a),  the  classical  and  QM 
models  of  1 If  drain  current  noise  behave  in  a qualitatively  similar  fashion  to  the  RTS 
noise  in  the  weak  and  strong  inversion,  showing  the  deviation  in  the  saturation  region. 
Figure  7.9  (b)  illustrates  the  thermal  noise  spectral  power  density  of  drain  current  as  a 
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function  of  the  drain-source  voltage  for  different  gate-source  voltage.  It  is  seen  that  QM 
effects  are  also  important  to  the  thermal  noise  modeling.  As  shown  in  Figure  7.10,  it  is 
confirmed  that  for  higher  doping  concentration  and  thinner  oxides,  QM  effects  in  the 
inversion  layer  should  be  considered  for  accurate  noise  modeling  of  nanoscale 
MOSFETs. 

7.6  Conclusion 

The  current  and  noise  models  accounting  for  QM  effects  in  nanoscale  MOSFETs 
were  presented.  Based  on  the  BSIM  framework,  simple  and  analytical  models  were 
formulated  by  considering  the  QM  effects  due  to  the  lowest  energy  level  in  the  inversion 
layer.  The  inversion  charge  parameters  were  studied  and  treated  quantum-mechanically, 
because  good  modeling  of  current  and  noise  depends  strongly  on  the  accurate  estimate  of 
the  inversion  charge  due  to  QM  effects.  The  gate  leakage  current  and  noise  models  were 
evaluated  by  using  the  quantum  correction  model  for  the  oxide  voltage  and  oxide 
thickness.  By  accounting  for  the  QM  effects  through  the  inversion  charge  density  model, 
the  quantum  models  of  the  drain  current  and  noise  were  obtained  and  compared  with  the 
classical  models.  The  simulation  results  showed  that  1)  the  quantum  correction  of  the 
oxide  thickness  is  essential  to  the  accurate  prediction  of  the  gate  leakage  current  and 
noise  model,  2)  QM  effects  on  the  1//  noise  and  the  thermal  noise  of  drain  current  is 
significant  in  the  saturation  region,  similar  to  the  drain  current  characteristics,  and  3)  QM 
effects  on  the  noise  model  become  more  pronounced  for  higher  doing  concentration  and 
thinner  oxide  thickness.  Therefore,  as  the  MOS  devices  are  scaled  down  to  nanometer 
regime,  existing  current  and  noise  models  should  incorporate  more  accurate,  physics- 
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based  descriptions  which  account  for  quantization  effects  observable  in  the  inversion 
layer. 


CHAPTER  8 

COMPACT  NOISE  MODELING  AND  IMPLEMENTATION 

8.1  Introduction 

Aggressive  downscaling  of  MOSFETs  has  been  driven  by  significant  advances  in 
CMOS  technologies,  leading  to  increasing  demands  for  compact  model  (CM) 
development.  The  compact  modeling  of  MOS  devices  is  critical  for  the  computer-aided 
design  (CAD)  of  digital  and  analog  VLSI  circuits  [Joa98,  Tsi94],  There  are  several 
recently  published  compact  models  of  MOSFETs,  such  as  the  BSIM4  model  from 
Berkeley  [BerOl],  the  MM11  model  from  Philips  [PhiOl],  the  EKV  model  [Enz96],  and 
the  MISNAN  model  [Boo91].  Among  them,  BSIM3v3  was  selected  as  the  first  and 
public  domain  MOSFET  model  for  standardization  by  the  compact  model  council 
(CMC).  Although  BSIM4  has  demonstrated  the  capability  to  predict  sub-micron  device 
characteristics,  it  is  still  necessary  to  develop  a generic  and  practical  model  flexible  to  the 
significant  changes  in  modem  CMOS  devices  with  nanometer  dimensions  [BerOl]. 
Especially,  an  accurate  and  physical  noise  modeling  is  essential  for  the  use  of  modem 
CMOS  devices,  because  the  noise  becomes  increasingly  larger  with  shrinking  devices 
[Cel99,  Sim99].  With  the  downscaling  of  gate  oxide  thickness,  the  effects  of  the  gate 
tunneling  current  noise  become  significantly  critical,  and  should  thus  be  incorporated  to 
accurately  predict  the  overall  noise  performance  [Mom98].  The  noise  models  should 
describe  all  kinds  of  the  noise  sources,  including  the  important  physical  effects  of 
nanoscale  MOSFET,  such  as  quantum  mechanical  effects  in  the  inversion  layer  [And82]. 
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Table  8. 1 Comparison  of  BSIM4  and  UF_Noise  for  noise  models 


MOS  Noise  Model 

BSIM4 

UF_Noise 

Gate  leakage  current 
noise 

1// noise  model 

Not  included 

Implemented 

Shot  noise  model 

Not  included 

Implemented 

1// noise  of 
drain  current 

Number  of  trap-related 
parameters 

3 

2 

Trapping-detrapping 

process 

Tunneling 

T unneling-thermal 

Trap  distribution 

Uniform 

Over  space  and  energy 

Frequency  exponent 

Default  (1.0) 

Bias-dependent 

Parasitic  resistance 

Not  included 

Included 

Correlation  noise 

Not  included 

Developed 

QM  correction  noise 

Not  included 

Developed 

Applicable  oxide  thickness 

> 3.0  nm 

< 3.0  nm 

As  shown  in  Table  8.1,  the  BSIM  does  not  contain  the  noise  model  of  gate 
leakage  tunneling  current  as  well  as  the  correlation  noise  model  between  the  gate  and 
drain,  which  are  no  longer  negligible  in  untrathin  oxide  MOSFETs.  In  addition,  BSIM 
noise  model  of  drain  current  has  several  drawbacks  which  are  physical  limitations  for 
trapping-detrapping  process,  unrealistic  trap  distribution  over  the  energy  and  space,  no 
capability  of  explaining  the  gate  bias  dependence  of  the  frequency  exponent,  and  no 
consideration  of  parasitic  series  resistance  [Hun90a,  Hun90b].  In  order  to  overcome  the 
problems  mentioned  above,  we  have  developed  the  physics-based  compact  model  of  the 
gate  leakage  current  and  noise  in  ultrathin  gate  oxide  MOSFETs,  in  Chapter  2 and  3.  It 
was  found  for  nitrided  oxide  that  the  results  simulated  by  the  proposed  current  model 
have  better  agreement  with  the  experimental  results  than  those  under  BSIM  model.  Also, 
the  newly  proposed  noise  model  can  accurately  predict  the  excess  noise  behavior  of  MOS 
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transistors  with  ultrathin  gate  oxide.  In  addition,  for  the  overall  noise  performance  of 
ultrathin  oxide  MOSFETs  at  low  frequencies,  the  correlation  noise  between  the  gate  and 
the  drain  should  be  modeled,  and  the  noise  model  has  been  developed  in  Chapter  5. 
Furthermore,  we  have  proposed  the  newly  developed  drain  current  noise  model  and  QM 
correction  model,  based  on  sound  physical  reasoning,  as  shown  in  Chapter  6 and  7. 

In  this  Chapter,  the  physics-based  compact  noise  model  will  be  expressed  in  such 
a way  that  the  models  incorporate  important  physical  effects  and  are  suitable  for  circuit 
simulators,  including  QM  effects  on  the  gate  leakage  current  noise,  the  drain  current 
noise,  and  the  correlation  noise.  An  emphasis  on  the  trap-related  parameters  of  noise 
models  is  laid  to  make  the  models  adaptable  to  the  variations  in  different  process 
technologies  and  to  make  its  parameters  easily  extractable  from  measured  data.  With  the 
help  of  an  accurate  and  generally  applicable  compact  noise  model,  we  can  successfully 
implement  the  compact  noise  models  into  BSIM  format. 

8.2  UF_Noise  : Gate  Leakage  Current  and  Noise  Models 

8.2.1  Model  Implementation  in  BSIM4 

The  gate  leakage  current  model  is  based  on  the  inelastic  trap-assisted  tunneling 
process  (ITAT).  The  model  describes  electrons  tunneling  to  deep-lying  trap  states  and, 
immediately,  released  to  positions  deeper  than  the  original  position,  and  subsequently 
tunneling  to  the  gate  [Lee02].  In  terms  of  trap-related  parameters,  X,,  4>,,  CCS,  and  NOID, 
defined  in  Table  8.2,  the  expressions  for  the  gate  leakage  current  density  with  the  drain- 
to-source  voltage  V*=  0 are  rewritten  as 
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where  9 is  the  electron  charge,  Pi  is  the  reduced  Planck  constant,  eox  is  the  oxide  dielectric 
constant,  <ph  is  the  Si-Si02  barrier  height,  mox  is  the  effective  mass  of  tunneling  carriers, 
t0- 1 is  the  oxide  thickness,  Eioss  is  the  energy  loss,  Fox  1,2  is  the  local  electric  field,  and  C is 
the  correction  function.  The  equations  includes  the  effects  of  charge  stored  in  the  trap 
states  and  the  dependence  of  the  energy  loss  on  the  oxide  electric  field  [Lee02]. 


Table  8.2  Trap-related  parameters  of  noise  model  of  gate  leakage  current 


Symbols 

Description 

Default 

Unit 

X, 

Trap  distance  from  Si-Si02  interface 

0.5  to. 

nm 

Barrier  height  of  the  oxide  trap  states 

2.5 

eV 

CCS 

Capture  cross  section 

5.0xl0'14 

cm2 

NOID 

Oxide  slow  trap  density 

2.0xl012 

cm'2  eV'1 

NOIE 

Oxide  fast  trap  density 

2.0xl012 

cm'2  eV'1 
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The  dependence  of  gate  leakage  current  on  the  drain-source  voltage  is  described 
by  the  source-drain  partition  model  which  has  been  already  implemented  in  BSIM4 
[CaoOO].  The  gate  tunneling  current  is  given  by 
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where  h is  the  Planck  constant,  W is  the  channel  width,  L is  the  channel  length,  Pigcd  is  a 
fitting  parameter  with  a default  value  of  1,  VY  is  the  threshold  voltage,  and  V»  is  the  gate- 
to- source  voltage. 

One  possible  mechanism  for  1 If  noise  contribution  is  for  traps  to  affect  locally 
the  barrier  height  or  shape  due  to  thermal  noise  of  the  frequency-dependent  conductance 
of  the  oxide  slow  traps  [Ghi87,  Kle78].  This  fluctuation  in  turn  modulates  the  tunneling 
transmission,  and  causes  the  fluctuation  in  the  tunneling  current.  Three  assumptions  are 
made  in  implementing  the  1 If  noise  model  of  gate  leakage  current  into  BSIM4.  1)  The 
oxide  slow  traps  are  uniformly  distributed  over  the  energy  and  space,  2)  The  parallel  trap 
conductance  and  capacitance  is  constant  over  low  frequencies,  and  3)  The  carrier 
tunneling  is  effective  for  trap  levels  around  the  quasi-Fermi  level.  Thus,  the  1 If  noise 
model  of  gate  leakage  current  can  be  simplified  as 
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with 

Vox=Vgs-VFB-il/s-Vpoly  (8.6) 

where  A,  = 2^2 mox<j)b  Ih  ,/  is  the  frequency,  C„x  is  the  oxide  capacitance,  Vox  is  the  oxide 

voltage,  if/s  is  the  surface  potential,  Vfb  is  the  flat-band  voltage,  and  Vpoiy  is  the  polysilicon 
voltage  drop. 

The  mechanism  for  shot  noise  associated  with  generation-recombination  process 
can  be  explained  by  carrier  transport  based  on  trap-assisted  tunneling  [IanOO,  Lee02]. 
The  generation  rate,  i.e.,  the  transition  rate  from  the  substrate  and  the  gate  to  the 
unoccupied  trap  in  the  oxide,  and  recombination  rate,  i.e.,  the  transition  rate  from  the 
occupied  trap  to  the  substrate  and  the  gate,  are  caused  by  the  statistical  occupancy 
fluctuation  of  fast  oxide  traps  in  the  oxide.  The  shot  noise  model  is  given  by 


Ss,ho,(f)  = 2qI  F 


Fano 


(8.7) 


with 


pi  Fano  y 


1 fast 


T fasti  1 + 4 7t2f2T 


2 *1,2 
fast 


+- 


* fast 
2 r2J2 


tl  1 4-  Ajr1  fLT 

L fasti  J i 


(8.8a) 


fast 


1 1 

— + — 


^ fast  ^ fast  1 ^ fast  2 


(8.8b) 


173 


(8.8c) 


= ±0.5  q-CCS  NOIE 

Tfas,  2 dN 


(8.8d) 


where  FFano  is  the  Fano  factor,  Tfasn,2  is  the  local  time  constants,  N is  the  total  number  of 


carriers,  and  g and  r are  the  uniform  generation  and  recombination  rate,  respectively.  The 
signs  of  Equations  (8.8c)  and  (8.8d)  depend  on  the  positive  or  negative  correlation 
between  tunneling  current  components.  As  illustrated  in  Figure  8.1,  for  example, 
different  signs  of  the  time  constants  lead  to  enhanced  shot  noise. 

Due  to  the  quantized  energy  levels  in  the  inversion  layer,  the  peak  charge  density 
is  situated  at  a distance  away  from  the  Si-SiC>2  interface.  This  results  in  the  enhancement 
in  the  oxide  thickness  and  the  reduction  in  the  oxide  voltage.  As  shown  in  Chapter  7, 
therefore,  the  quantum  effects  are  taken  into  account  by  using  the  correction  for  the  oxide 
thickness,  t , and  the  oxide  voltage,  V°M  , as  follows 


where  xf*  is  the  average  centroid  position  from  the  interface  and  t®xM  is  the  QM 
corrected  surface  potential,  as  defined  in  Chapter  7. 

8.2.2  Benchmark  Tests 

For  the  practical  application  of  the  compact  noise  model  of  gate  leakage  current, 


fQM  | £ox  rQM 

lox  ~lox^  „ An 
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(8.9a) 


(8.9b) 


the  models  are  implemented  into  BSIM4  in  public  domain,  as  shown  in  Figure  8.2. 
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x 10'5 


Figure  8.1  Enhancement  and  suppression  in  the  Fano  factor,  depending  on  the  signs  of 
local  time  constants  (/=  10  kHz ) 
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Experimental  Verification  for  Current  and  Noise 


Parameter  Extraction 


Trap-related 

Parameters 


Basic  MOSFET  Parameters 


Trap-related 

Parameters 


BSIM4  Standard  Format 


Figure  8.2  Implementation  of  the  gate  leakage  current  and  noise  models  into  the  BSIM4 
public  domain  platform 
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fdd 


Vds 


(a) 


(b) 


Figure  8.3  Benchmark  test  of  the  noise  model  of  gate  leakage  current,  (a)  a test  circuit 
of  a single  MOSFET  and  (b)  an  input  deck  to  simulate  the  gate  leakage 
current  noise 
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(a) 


(b) 


Figure  8.4  Comparison  of  the  benchmark  test  results  and  measured  data  of  the  gate 
leakage  current,  (a)  the  current  density  versus  the  gate-source  voltage,  and 
(b)  the  simulated  current  versus  the  drain-source  voltage  for  different  gate- 
source  voltages  ( W = 10  /urn,  L = 10//m,  toxeq  = 2.2  nm,  NOID  = 2xl0'‘  cm" 
2tWl ). 
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(a) 


(b) 

Figure  8.5  Comparison  of  the  benchmark  test  results  and  measured  data  of  the  gate 
leakage  current  noise,  (a)  the  noise  characteristics  as  a function  of  frequency 
for  different  gate-source  voltages  and  (b)  the  simulated  noise  characteristics 
as  a function  of  drain-source  voltage  for  two  channel  lengths  ( tox,eq  = 2.2  nm, 
NOID=NOIE  = 2xl012  cm'2eV‘ ) 
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The  basic  device  parameters  are  also  used  for  the  current  and  noise  models,  but 
the  trap-related  parameters  are  newly  introduced  into  simulators.  The  validation  for  the 
compact  model  is  derived  from  the  comparisons  of  simulation  results  and  measurements. 
All  tests  are  to  be  done  using  a test  circuit  of  a single  MOSFET,  as  depicted  in  Figure  8.3 

(a) .  An  input  deck  in  SPICE  format  is  provided  in  Figure  8.3  (b).  Figure  8.4  (a)  and  (b) 
show  the  comparison  of  the  test  results  and  measured  data  of  gate  leakage  current  as  a 
function  of  gate-source  voltage  and  drain-source  voltage,  respectively.  It  is  observed  that 
the  simulation  results  have  good  agreement  with  measurements.  Also,  Figure  8.5  (a)  and 

(b)  illustrate  the  noise  characteristics  as  a function  of  frequency  for  different  gate-source 
voltages  and  as  a function  of  drain-source  voltage  for  two  channel  lengths,  respectively. 
It  is  also  verified  that  the  compact  noise  model  can  accurately  predict  the  noise  behavior 
over  the  frequency  and  bias  ranges.  As  shown  in  Figure  8.5  (b),  the  downscaling  of  the 
channel  length  reduces  the  gate  leakage  current  noise  at  a given  oxide  thickness. 
However,  since  the  drain  current  noise  increases  with  downscaling  of  channel  length,  a 
trade-off  is  needed  to  optimize  the  overall  noise  performance. 

8.3  UF_Noise  : Drain  Current  and  Correlation  Noise  Models 
8.3.1  Drain  Current  Noise  Model 

In  order  to  explain  the  1 If  noise  of  drain  current  in  ultrathin  oxide  MOSFETs,  the 
tunneling  assisted-thermally  activated  (TATA)  time  constant  model  is  employed,  which 
is  physically  consistent  with  the  gate  leakage  tunneling  process  [Lee02,  Tew94].  Since 
the  noise  is  strongly  affected  by  trap-related  properties,  the  noise  model  accounts  for  a 
realistic  description  of  trap  distribution.  In  addition,  the  model  should  include  the  effects 
of  the  two-dimensional  electric  fields,  the  parasitic  resistance,  the  doping  profile,  and  the 
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mobility  degradation  on  the  noise  behavior,  as  described  in  Chapter  6.  The  QM 
correction  noise  model  can  be  obtained  by  accounting  for  the  QM  effects  through  the  in 
the  oxide  capacitance,  the  threshold  voltage,  and  the  drain  current  in  the  charge  density 
parameters  as  well  as  the  gate  oxide  thickness  [MaOl].  The  QM  correction  model  for  1 If 
drain  current  noise  in  the  linear  region,  accounting  for  all  effects  mentioned  above,  is 
given  by 
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where  Vto  is  the  threshold  voltage  with  VdS  = 0 V,  fauik  is  the  bulk  charge  parameter,  R is 
the  parasitic  resistance,  fii  and  /?2  are  the  mobility  parameters,  /r„o  is  the  maximum 
mobility  in  the  inversion  layer,  Teffo  is  an  effective  attempt  time  constant,  NTa,  Ntb,  and 
Ntc  depend  on  the  trap  distribution  in  the  oxide,  and  Xth,eff  and  A tUn,eff  are  an  exponential 
factor  of  thermal  activation  process  and  tunneling  process,  respectively.  Unlike  BSIM 
noise  model,  A*,  B*,  and  C*  are  not  empirical  constants  but  are  analytical  expressions  as 
a function  of  Vv  and  Vds,  as  described  in  Chapter  6. 


181 


In  the  saturation  region,  the  noise  power  spectral  density  is  expressed  as  the  sum 
of  the  noise  components  due  to  the  triode  region,  L-Lsat,  and  the  pinch-off  region,  Lsat, 
and  is  given  by 


C(/)=c»(/)  L 
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where  is  the  drain  saturation  voltage,  and  A^a,,  B*Ual,  and  C*Uat  are  analytically 
expressed  as  a function  of  Vgs  and  V*ar. 

8.3.2  Comprehensive  noise  model 

Due  to  the  presence  of  the  gate  leakage  current  in  ultrathin  oxide  MOSFETs,  the 
comprehensive  noise  performance  at  low  frequencies  is  evaluated  by  taking  into  account 
the  cross  correlation  coefficient,  &F(f),  between  the  drain  and  the  gate  current  noise,  and 
is  expressed  as  a function  of  the  current  ratio,  Ap=Ig/Id,  and  the  noise  ratio, 

SFF(f)/SFF(f),  as  follows 
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In  addition  to  the  influence  of  the  correlation  on  the  low  frequency  noise 

performance,  the  shot  noise  due  to  the  gate  leakage  current,  Ss.ho'(f ) = 2 ql  F Fano , has  an 

1 8 8 

impact  on  the  noise  behavior,  and  thus  should  be  incorporated  into  the  minimum  noise 
figure  F^fn  (/)  as  follows  [Cap88,  Puc74] 
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where  8 mg  =dId  !dvgs  is  the  transconductance,  gg  = dIg/dVgs  is  the  gate  conductance, 
and  Cgs  is  the  gate-to-source  capacitance. 

8.3.3  Simulation  results 

From  the  viewpoint  of  compact  noise  modeling,  the  model  should  provide  an 
accurate  prediction  of  noise  behavior,  describing  all  the  important  physical  effects 
without  an  increase  of  parameters.  For  nMOSFETs  with  ultrathin  oxide  ( = 2.2nm  ),  the 
1 //  noise  of  drain  current  can  be  accurately  predicted  by  the  expressions  of  (8.10)  - 
(8.13).  Figure  8.6  shows  the  comparison  of  measured  data  and  simulation  results  carried 
out  with  single  trap-related  parameter  under  the  assumption  of  uniform  trap  distribution 
in  the  oxide.  Excellent  agreements  between  measurement  and  simulation  are  observed 
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(a) 


(b) 


Figure  8.6  Comparison  of  the  simulation  results  and  measured  data  of  the  drain  current 
noise,  (a)  the  noise  characteristics  as  a function  of  frequency  for  two 
channel  lengths,  (b)  the  noise  characteristics  as  a function  of  drain-source 
voltage  for  two  channel  lengths,  assuming  that  the  trap  is  uniformly 
distributed.  ( tox,eq  = 2.2  nm,  Nta  = lxlO18  cm^eV'1 ). 
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(a) 


(b) 

Figure  8.7  Simulation  of  the  minimum  noise  figure,  (a)  the  characteristics  of  minimum 
noise  figure  as  a function  of  frequency  for  different  gate-source  voltagse,  (b) 
the  characteristics  of  minimum  noise  figure  as  a function  of  gate-source 
voltage  for  different  drain-source  voltages  ( tox,eq  =2.2  nm,  NTa=  lxlO18  cm' 
3eV\  NOID  = 2xl012  cm'VV'1 ). 
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over  the  frequency  and  in  both  the  linear  and  saturation  regions  for  two  channel  lengths, 
as  shown  in  Figure  8.6  (a)  and  (b),  respectively. 

For  the  overall  noise  evaluation  of  advanced  CMOS  devices,  the  compact  noise 
models  at  low  frequencies  should  include  the  correlation  effects  between  the  gate  and  the 
drain,  together  with  the  gate  and  drain  current  noise  models.  Unlike  the  capacitive 
coupling  at  high  frequencies,  the  low  frequency  correlation  is  mainly  attributed  to  the 
trap-assisted  tunneling  process,  which  is  a common  origin  of  the  gate  and  drain  current 
noise.  With  the  help  of  an  analytical  expression  of  correlation  coefficient,  it  is  possible  to 
evaluate  the  minimum  noise  figure,  as  illustrated  in  Figure  8.7.  Figure  8.7  (a)  and  (b) 
show  the  characteristics  of  minimum  noise  figure  as  a function  of  frequency  for  different 
gate-source  voltages  Vgs  and  the  gate-source  voltage  for  different  drain-source  voltages 
Vds,  respectively.  The  increasing  minimum  noise  figure  with  increasing  V*  and 
decreasing  Vd  can  be  interpreted  as  being  due  to  the  gate  leakage  current  and  its  noise 
components,  greatly  affecting  the  overall  noise  performance  at  low  frequencies. 

8.4  Conclusion 

The  comprehensive  and  physics-based  compact  noise  models  for  advanced 
CMOS  devices  were  presented.  The  models  incorporate  important  physical  effects  in 
nanoscale  MOSFETs,  such  as  the  gate  leakage  current  and  noise,  the  low  frequency 
correlation  effect  between  the  drain  and  the  gate,  the  trap-related  phenomena,  and  QM 
effects  in  the  inversion  layer.  Especially,  the  gate  leakage  current  and  noise  models  were 
implemented  in  BSIM4  simulator  and  were  verified  with  the  comparison  of  simulation 
results  and  measurements.  For  ultrathin  oxide  MOSFETs,  the  drain  current  noise  model 
was  improved  by  including  the  tunneling  assisted-thermally  activated  process,  the 
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realistic  trap  distribution,  the  parasitic  resistance,  and  mobility  degradation.  The 
expression  of  correlation  coefficient  was  analytically  described,  enabling  the  overall 
noise  performance  to  be  evaluated.  With  the  consideration  of  QM  effects,  the 
comprehensive  low  frequency  noise  performance  was  simulated  over  the  entire  bias 
range.  Based  on  the  simulation  results  and  measurements,  it  is  found  that  as  modem 
CMOS  technology  scales  down  to  nanometer  dimension,  the  noise  model  should  be 
replaced  by  more  accurate  and  generally  applicable  physics-based  model. 


CHAPTER  9 

SUMMARY  AND  RECOMMENDATIONS  FOR  RELATED  WORK 

9.1  Summary 

The  research  presented  in  this  dissertation  explored  the  development  of  physics- 
based  compact  noise  models  of  nanoscale  MOSFETs,  which  are  applicable  to  circuit 
simulation  and  can  be  used  to  upgrade  the  existing  noise  models  in  the  BSIM  simulator. 
The  models  include  the  noise  behavior  due  to  the  gate  leakage  current,  the  drain  current, 
the  cross  correlation  between  the  gate  and  the  drain,  and  the  quantization  in  the  inversion 
layer  of  ultrathin  oxide  MOS  devices.  It  was  also  demonstrated  in  Chapter  4 that  1 If 
noise  of  the  gate  leakage  current  is  directly  applicable  to  extract  the  trap  properties  of  an 
ultrathin  oxide  MOS  structure.  Especially,  the  compact  noise  model  of  the  gate  leakage 
current  was  newly  implemented  into  BSIM4,  and  the  utility  of  the  modeling  was  verified 
by  a comparison  of  the  simulation  results  with  measurements,  as  shown  in  Chapter  8. 

A short  description  of  the  noise  models  developed  in  this  work  will  be  outlined 
and  improvements  are  highlighted  in  terms  of  noise  modeling.  The  models  offer  the 
following  unique  features : 

1)  For  the  characterization  of  the  enhanced  conduction  observed  in  the  thin 
nitrided  oxide  devices,  the  gate  leakage  current  model  is  based  on  inelastic  trap-assisted 
tunneling  mechanism  combined  with  a semi-empirical  direct  tunneling  current  model, 
including  the  effects  of  charge  stored  in  the  trap  states  and  the  dependence  of  the  energy 
loss  on  the  oxide  field,  as  shown  in  Chapter  2. 
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2)  A physics-based  noise  model  of  gate  leakage  current  is  developed,  which  is 
physically  based  on  the  trap-related  tunneling  processes,  as  described  in  Chapter  3.  In 
order  to  explain  the  excess  noise  behavior,  low  and  high  frequency  noise  models  are 
formulated  using  a barrier  height  fluctuation  model  and  Lorentzian-modulated  shot  noise 
model  for  the  gate  leakage  current,  respectively. 

3)  In  Chapter  5,  for  a comprehensive  evaluation  of  low  frequency  noise 
performance  in  ultrathin  gate  oxide  MOSFETs,  the  model  includes  the  gate  leakage 
current  and  the  cross  correlation  coefficient  between  the  drain  and  the  gate.  By  using  the 
partition  noise  theory  and  the  BSIM4  partition  current  model,  the  correlation  coefficient 
at  low  frequencies  is  modeled.  The  noise  performance  characteristics  at  low  frequencies, 
such  as  minimum  noise  figure,  optimum  source  resistance  and  reactance,  are  evaluated  by 
introducing  the  frequency-dependent  noise  parameters,  P^F , RLF,  and  CfF. 

4)  A Ilf  drain  current  noise  model  of  ultrathin  oxide  MOSFETs  operating  in  both 
the  linear  and  saturation  regions  is  formulated  on  the  basis  of  the  tunneling  assisted- 
thermally  activated  process  of  carriers,  incorporating  a realistic  description  of  trap 
distributions,  as  shown  in  Chapter  6.  The  model  incorporates  the  effects  of  the  two- 
dimensional  electric  field,  the  parasitic  resistance,  the  doping  profile,  and  the  quadratic 
mobility  degradation  factor  on  the  noise  behavior. 

5)  In  Chapter  7,  quantum-mechanical  (QM)  effects  on  the  noise  characteristics 
have  been  taken  into  account  for  the  first  time,  extending  the  applicability  of  the  model  to 
ultrathin  («1.5nm)  MOS  devices.  By  considering  the  QM  effects  due  to  the  lowest 
energy  level  in  the  inversion  layer,  quantum  noise  models  are  obtained  and  compared 
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with  classical  models.  It  is  observed  that  QM  effects  on  the  noise  model  become  more 
pronounced  for  higher  doing  concentration  and  thinner  oxide  thickness. 

6)  With  the  help  of  analytical  and  physical  noise  models  applicable  to  nanoscale 
MOSFETs,  the  noise  models  can  be  implemented  into  general-purpose  circuit  simulators, 
providing  a comprehensive  evaluation  of  noise  characteristics  in  advanced  CMOS 
devices. 

9.2  Recommendations  for  Future  Work 

Although  the  noise  models  in  this  dissertation  are  focusing  on  the  noise  behaviors 
of  n-channel  MOSFETs,  the  models  are  equally  applicable  to  p-channel  MOSFETs  if  the 
appropriate  changes  for  voltage  polarity,  barrier  height,  and  some  physical  parameters  are 
made.  Nevertheless,  more  work  has  to  be  done  on  the  model  verification  by  noise 
measurements  over  the  entire  bias  range. 

As  described  in  Chapter  3,  the  shot  noise  due  to  the  gate  leakage  current  is 
suppressed  or  enhanced,  depending  on  the  sign  of  the  characteristic  time  constants.  While 
the  Lorentzian-modulated  shot  noise  model  can  accurately  predict  the  noise  behavior,  it  is 
still  necessary  to  identify  the  physical  origin  of  transport  process  responsible  for  the  shot 
noise.  A microscopic  analysis  of  shot  noise  can  be  done  by  combining  Coulomb 
interaction  with  trapping-detrapping  process  in  a MOS  structure,  and  simulations  can  be 
performed  by  means  of  an  ensemble  Monte  Carlo  simulator  self-consistently  coupled 
with  a Poisson  solver  [Vas02]. 

Besides  this  work  on  compact  noise  modeling,  future  research  must  explore  the 
RTS  noise  modeling,  the  RF  compact  noise  modeling,  and  low  noise  design.  Due  to  the 
continuous  downscaling  of  CMOS  devices,  it  may  occur  that  only  one  oxide  trap  with  an 


190 


energy  within  a few  kT  of  the  surface  Fermi  level  exists.  Under  these  conditions,  the 
current  fluctuations  behave  discretely  and  are  known  as  RTS  noise.  The  analysis  of  RTS 
noise  provides  the  opportunity  of  studying  the  interaction  of  inversion  carriers  with 
interface  states,  and  the  RTS  noise  is  also  strongly  correlated  with  the  presence  of  an 
excess  leakage  current.  Therefore,  an  accurate  and  general  RTS  noise  model  is  needed  to 
predict  the  noise  behavior  of  future  devices.  In  addition,  continuous  downscaling  of 
CMOS  devices  makes  CMOS  technology  a viable  choice  for  RF  applications. 
Nevertheless,  the  design  of  RF  ICs  for  real  products  remain  a challenge  due  to  the  strong 
constraints  on  noise.  Particularly,  for  analog  applications,  such  as  front-end  circuits  for 
wireless  communication,  noise  issues  are  of  primary  importance.  This  requires  the 
development  of  RF  compact  noise  models  for  advanced  CMOS  devices.  Furthermore, 
from  the  practical  viewpoint  of  optimum  noise  performance,  the  noise  should  be  reduced 
as  a function  of  the  process  technology  [CheOlb,  Jay90,  Lai99,  ToiOl],  channel 
engineering  [Che02,  Chr02,  Mat97,  Okh99,  Yos98],  and  circuit  techniques  [Blo91, 


Gie99,  Hsu99,  KluOO,  WelOO]. 
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